Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,506,454 C→T 22.7% intergenic (‑269/‑67) mntH ← / → nupC manganese/divalent cation transporter/nucleoside (except guanosine) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,506,4540CT22.7% 108.7 / 30.1 63intergenic (‑269/‑67)mntH/nupCmanganese/divalent cation transporter/nucleoside (except guanosine) transporter
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base T (5/9);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

TTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGACCGCAAAA  >  CP009273/2506310‑2506599
                                                                                                                                                |                                                                                                                                                 
ttGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGttt                                                                                                                                               <  1:763698/149‑1 (MQ=255)
    tgaCTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCt                                                                                                                                           >  1:1027977/1‑149 (MQ=255)
      aCTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTa                                                                                                                                         <  2:1131404/149‑1 (MQ=255)
       cTTTAATGTTATGGTTATGTTTAATGGTAGCACATTTACATAACTTTCAGCGTCCATACACAACAGAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTaa                                                                                                                                        <  1:769153/149‑1 (MQ=255)
             tGTTATGGTTATGTTTAATGTTAGCACATTTACATAACTTTCGGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGt                                                                                                                                  <  2:75950/149‑1 (MQ=255)
                          tttAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAACTGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAAtgtg                                                                                                                     >  2:998704/1‑149 (MQ=255)
                          tttAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAAtgtg                                                                                                                     >  2:1205930/1‑149 (MQ=255)
                          tttAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAAtgtg                                                                                                                     >  1:398494/1‑149 (MQ=255)
                          tttAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGAGCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTACCCATGTTACATAGAAtgtg                                                                                                                     >  1:486418/1‑149 (MQ=255)
                          gttAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATCGTGTGTGTCAGATCTCGTTTTCCTTAACCATCTTCCCTGGAAtgtg                                                                                                                     >  2:943849/2‑149 (MQ=255)
                           ttAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGc                                                                                                                    <  2:714163/149‑1 (MQ=255)
                             aaTGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCAc                                                                                                                  <  2:1023199/149‑1 (MQ=255)
                             aaTGTTAGAACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCAc                                                                                                                  <  1:168960/149‑1 (MQ=255)
                                gtcagcatgtttACTCAACTTTCAGCTTCCATACACAACCTAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTTACCATGTTACATAGAATGTGCACGGa                                                                                                               <  2:238775/140‑1 (MQ=255)
                                gctAGAACATTTAAATAACTTTCAGATTCCATACACAACATAGCAGAATTGTATGACAGATCAGTATTTTTGAAGTCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGa                                                                                                               <  2:584308/147‑1 (MQ=255)
                                gatAGCTCATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGa                                                                                                               <  2:214552/147‑1 (MQ=255)
                                    gCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGGAAtt                                                                                                           >  1:632228/1‑149 (MQ=255)
                                    gCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAAtt                                                                                                           >  1:1186683/1‑149 (MQ=255)
                                     cacaTTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCCGACCACGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAAttt                                                                                                          >  1:231187/1‑149 (MQ=255)
                                        aTTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAAc                                                                                                       <  1:1050387/149‑1 (MQ=255)
                                          tgaCATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCt                                                                                                     <  2:230440/147‑1 (MQ=255)
                                              aTAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCt                                                                                                 >  2:1035090/1‑149 (MQ=255)
                                              aTAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCt                                                                                                 >  1:273246/1‑149 (MQ=255)
                                               tAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGCAGTCATTATAGTGTGTGTCAGATCTCGTTTGCCTTAACCATGTTAAATAGAATGTGCCCGGCAATTTAACCTGCCTc                                                                                                >  2:669624/1‑149 (MQ=255)
                                                   tttCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCatat                                                                                            >  1:640289/1‑149 (MQ=255)
                                                     tCAGCTTCCATACACAACATAGCGGAAATGTATGACAGGTCACTATTTTTGAAGCCTGTCACAGGACGGCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATAttt                                                                                          >  1:668703/1‑149 (MQ=255)
                                                     gcAGCTGCCATACACAACATAGCAGAAATGTATGACAGATCACTGTTTGTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATAttt                                                                                          <  1:1089544/148‑1 (MQ=255)
                                                      cAGGTTCCATACACAACTTAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGGCAGATCTTGTTTTCCTTATCCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTg                                                                                         <  2:1019702/149‑1 (MQ=255)
                                                       aGCTTCCAGAGACAACAGATCGTAACGGTATGACAGATCGCTCTTTTTGAAGCCTGGCACAGGAGGTCATTATACTGTGTGTCAGATCTGGTCTTCCTTAACTATGTTACATAGAATGTGCACGGACATTTAACCTGCCTCATATTTgg                                                                                        <  2:1011014/149‑1 (MQ=255)
                                                            ccATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCCCGGAAATTTAACCTGCCTCATATTTAGAGCaa                                                                                   >  2:657607/1‑149 (MQ=255)
                                                                 tgcaacaTAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTACCCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATg                                                                              <  2:741516/147‑1 (MQ=255)
                                                                 cacaACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTGAACCATGTTACATAGAATGTGCAAGGGATTTTTACCTGCCTCATATTTGGAGCAAATATg                                                                              <  2:881987/149‑1 (MQ=255)
                                                                   caacaTAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGa                                                                            >  1:554137/1‑149 (MQ=255)
                                                                    aacaTAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGAc                                                                           >  2:19530/1‑149 (MQ=255)
                                                                         aGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGt                                                                      <  1:1205930/149‑1 (MQ=255)
                                                                          gtaGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGCGTGTGTCAGATCTCGTTTTCCCTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTc                                                                     <  2:748461/147‑1 (MQ=255)
                                                                           cAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTcc                                                                    >  1:849997/1‑149 (MQ=255)
                                                                               aaTGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCa                                                                <  2:1158736/149‑1 (MQ=255)
                                                                                 tGTATGACAGATCACTATTTTTGAAGCCGGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCAGATTTGGAGCAACTATGGACCACGTCCTTCAtt                                                              >  2:699917/1‑149 (MQ=255)
                                                                                   tctGACAGATCACTATTTTTGAATCCGGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCAtttt                                                            <  2:127763/147‑1 (MQ=255)
                                                                                   tATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCAtttt                                                            >  1:226485/1‑149 (MQ=255)
                                                                                   tATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCAtttt                                                            >  1:1154605/1‑149 (MQ=255)
                                                                                                tATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACCTAGAATGTGCACGGAAATTTAACCTGCCTCATCTTTGCAGCAATTATGGACCGCGTCCTTCCTTCTGTACTGGCACTTg                                               >  2:170395/1‑149 (MQ=255)
                                                                                                tATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTg                                               >  1:567548/1‑149 (MQ=255)
                                                                                                    tttGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCgt                                           >  1:1047145/1‑149 (MQ=255)
                                                                                                         aGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGttgttg                                      <  2:298857/149‑1 (MQ=255)
                                                                                                                  caGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCg                             <  1:637496/149‑1 (MQ=255)
                                                                                                                   aGGACGTCTTTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGc                            <  2:51892/149‑1 (MQ=255)
                                                                                                                   aGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGc                            >  2:1008670/1‑149 (MQ=255)
                                                                                                                         tCATTATAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTCCACGGAAATTTACCCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCActgct                      <  1:699917/149‑1 (MQ=255)
                                                                                                                          cATTTTAGTGTGTGTCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAACATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGTTAATCGCActgcgg                     <  2:937337/149‑3 (MQ=255)
                                                                                                                          cATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCActgctg                     >  2:1125622/1‑149 (MQ=255)
                                                                                                                          cATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTCGTTGCGATTCTCGCActgctg                     >  2:282593/1‑149 (MQ=255)
                                                                                                                            gtatggcgtgtgtCAGATCTCGTTTGCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGt                   <  2:779545/142‑1 (MQ=255)
                                                                                                                              ataGTGTGTGTCAGAACTCGTTTTCCTTAACCATGTTACACAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTaa                 <  1:981555/149‑1 (MQ=255)
                                                                                                                                aGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAagc               <  1:501576/149‑1 (MQ=255)
                                                                                                                                    tgtgtCAGATCTTGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCg           >  1:403541/1‑149 (MQ=255)
                                                                                                                                     gtgtCAGACCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGa          <  1:485376/149‑1 (MQ=255)
                                                                                                                                       ttCAGACCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCTTCATATTTGGAGCAAATATGGACCGCGTCCTTCTTTTTGTACTGGCCCTTGCCGTTGTTGCGATTTTCGCATTGCTGGTAAGCAGCGAcc        <  2:605681/148‑1 (MQ=255)
                                                                                                                                       gtCAGATTGTGTTTGCCTTCACCATGTCACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTAGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGAcc        <  2:824372/149‑1 (MQ=255)
                                                                                                                                        tCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGACCg       >  2:637634/1‑149 (MQ=255)
                                                                                                                                            aTCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGACCGCaaa   >  1:1121968/1‑149 (MQ=255)
                                                                                                                                            aTCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGACCGCaaa   >  2:1125797/1‑149 (MQ=255)
                                                                                                                                             tctcGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGNACTNCNGGTAAGCACCGACCGCaaaa  >  1:675644/1‑149 (MQ=255)
                                                                                                                                             tctcGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGACCGCacaa  >  2:215606/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTGATGACTTGAATGTTATGGTTATGTTTAATGTTAGCACATTTACATAACTTTCAGCTTCCATACACAACATAGCAGAAATGTATGACAGATCACTATTTTTGAAGCCTGTCACAGGACGTCATTATAGTGTGTGTCAGATCTCGTTTTCCTTAACCATGTTACATAGAATGTGCACGGAAATTTAACCTGCCTCATATTTGGAGCAAATATGGACCGCGTCCTTCATTTTGTACTGGCACTTGCCGTTGTTGCGATTCTCGCACTGCTGGTAAGCAGCGACCGCAAAA  >  CP009273/2506310‑2506599

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: