Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 882,672 G→T 100% L42I (CTC→ATC)  ybjJ ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273882,6720GT97.9% 610.2 / 3.5 193L42I (CTC→ATC) ybjJputative transporter
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base T (101/88);  total (101/92)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 198 total aligned reads displayed.

CCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAG  >  CP009273/882527‑882816
                                                                                                                                                 |                                                                                                                                                
ccGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATaa                                                                                                                                                <  1:1878750/148‑1 (MQ=255)
   aTCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAatat                                                                                                                                             >  2:3366794/1‑148 (MQ=255)
   aTCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAatat                                                                                                                                             <  2:3026611/148‑1 (MQ=255)
    tCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATc                                                                                                                                            >  1:2957646/1‑148 (MQ=255)
       aTGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATcgcg                                                                                                                                         >  2:3339848/1‑148 (MQ=255)
        tGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCgg                                                                                                                                        >  2:2120425/1‑148 (MQ=255)
         gcgcAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGa                                                                                                                                       <  2:900991/148‑1 (MQ=255)
          cgcAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGAt                                                                                                                                      >  1:1897332/1‑148 (MQ=255)
           gcAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATa                                                                                                                                     <  2:3247763/148‑1 (MQ=255)
            cAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAg                                                                                                                                    <  2:2676216/148‑1 (MQ=255)
             aGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGc                                                                                                                                   >  1:666931/1‑148 (MQ=255)
               gACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAg                                                                                                                                 >  1:1528229/1‑148 (MQ=255)
               gACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAg                                                                                                                                 >  2:351635/1‑148 (MQ=255)
                aCATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAgg                                                                                                                                >  2:25336/1‑148 (MQ=255)
                aCATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAgg                                                                                                                                >  2:25049/1‑148 (MQ=255)
                    cGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTa                                                                                                                            <  1:1165480/148‑1 (MQ=255)
                     gTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTAc                                                                                                                           <  2:3477011/148‑1 (MQ=255)
                      tGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACg                                                                                                                          >  1:533533/1‑148 (MQ=255)
                         ccAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGt                                                                                                                       <  2:28790/148‑1 (MQ=255)
                           aGGATGACGTTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTCCGGGTTg                                                                                                                     >  2:2618241/1‑148 (MQ=255)
                           aGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTg                                                                                                                     <  1:2233358/148‑1 (MQ=255)
                             gATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGcc                                                                                                                   >  1:169444/1‑148 (MQ=255)
                               tGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCa                                                                                                                 <  1:1214648/148‑1 (MQ=255)
                                gACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCCGGCGTACGGGTTGCCCAg                                                                                                                >  2:2943739/1‑148 (MQ=255)
                                 aCATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAgg                                                                                                               <  1:522045/148‑1 (MQ=255)
                                 aCATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAgg                                                                                                               <  2:20716/148‑1 (MQ=255)
                                   aTTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACTGGTTGCCCAGGAc                                                                                                             >  1:544147/1‑148 (MQ=255)
                                   aTTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGAc                                                                                                             >  2:1039221/1‑148 (MQ=255)
                                     tACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGc                                                                                                           <  2:1228075/148‑1 (MQ=255)
                                     tACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGc                                                                                                           <  1:1785541/148‑1 (MQ=255)
                                      aCGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGcc                                                                                                          >  2:1639801/1‑148 (MQ=255)
                                      aCGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGcc                                                                                                          <  1:3401248/148‑1 (MQ=255)
                                       cGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCa                                                                                                         >  2:2459721/1‑148 (MQ=255)
                                         tgtCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCAtt                                                                                                       <  1:38411/148‑1 (MQ=255)
                                            cccAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAAc                                                                                                    <  1:3334730/148‑1 (MQ=255)
                                            cccAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAAc                                                                                                    <  2:3091270/148‑1 (MQ=255)
                                                aaGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGc                                                                                                >  2:2224151/1‑148 (MQ=255)
                                                 aGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGcc                                                                                               >  1:1911974/1‑148 (MQ=255)
                                                  gCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCt                                                                                              <  1:1015731/148‑1 (MQ=255)
                                                  gCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCt                                                                                              <  2:1739757/148‑1 (MQ=255)
                                                     tttCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCCTGACCAGGCCTGGc                                                                                           >  1:287172/1‑148 (MQ=255)
                                                     tttCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGc                                                                                           >  1:823022/1‑148 (MQ=255)
                                                      ttCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCa                                                                                          >  2:396779/1‑148 (MQ=255)
                                                       tCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCaa                                                                                         >  1:141645/1‑148 (MQ=255)
                                                            aaCCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAaaga                                                                                    <  2:854991/148‑1 (MQ=255)
                                                            aaCCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAaaga                                                                                    <  1:668179/148‑1 (MQ=255)
                                                               cACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAagaaga                                                                                 <  1:155525/148‑1 (MQ=255)
                                                                aCGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAgaagaa                                                                                <  2:3106336/148‑1 (MQ=255)
                                                                aCGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAgaagaa                                                                                <  2:2296418/148‑1 (MQ=255)
                                                                 cGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAAc                                                                               <  1:2876134/148‑1 (MQ=255)
                                                                     gAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATaa                                                                           <  2:1213184/148‑1 (MQ=255)
                                                                      aGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATaaa                                                                          <  1:97574/148‑1 (MQ=255)
                                                                       gCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAAc                                                                         <  1:1311611/148‑1 (MQ=255)
                                                                          gagaATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCCTTAACAGGCCTGGCAAAAAGAAGAACATAAACAgc                                                                      >  2:1894757/1‑148 (MQ=255)
                                                                          gagaATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAgc                                                                      >  2:1394576/1‑148 (MQ=255)
                                                                           agaATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAgcg                                                                     <  1:3339848/148‑1 (MQ=255)
                                                                           agaATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAgcg                                                                     >  2:2922733/1‑148 (MQ=255)
                                                                            taaTACCGATAATCTAACCGATCGACAGACCAAAGAGAACACCGCCAATTTCAGCGATCTAGCCAGAGATAATATCGCGGATAGCTGGCGTACGGGTTGCCCAGGACGCAATTAACAGGCCTGGCAAAAAGAAGAACATAAACAgcgc                                                                    <  1:3280430/147‑1 (MQ=255)
                                                                               tACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCa                                                                 >  1:2456304/1‑148 (MQ=255)
                                                                                        tCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCg                                                        <  2:1433917/148‑1 (MQ=255)
                                                                                           aaCCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGttt                                                     >  1:301271/1‑148 (MQ=255)
                                                                                           aaCCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCAGTACCAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGttt                                                     >  2:1838901/1‑148 (MQ=255)
                                                                                            aCCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTc                                                    >  1:3240118/1‑148 (MQ=255)
                                                                                             ccGATCGACAGACCAAAGAGAACAGCGCCCATTTCAGCGTTCGAGACAGAGATAATAGCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCa                                                   >  1:1349734/1‑148 (MQ=255)
                                                                                             ccGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCa                                                   >  1:1349771/1‑148 (MQ=255)
                                                                                             ccGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCa                                                   <  1:1492002/148‑1 (MQ=255)
                                                                                             ccGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCa                                                   >  2:531273/1‑148 (MQ=255)
                                                                                              cGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCaa                                                  <  2:1288028/148‑1 (MQ=255)
                                                                                              cGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCaa                                                  <  1:798101/148‑1 (MQ=255)
                                                                                                aTCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATg                                                >  2:1614364/1‑148 (MQ=255)
                                                                                                aTCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATg                                                >  1:67394/1‑148 (MQ=255)
                                                                                                     cAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTa                                           >  2:2893623/1‑148 (MQ=255)
                                                                                                       gACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACg                                         >  1:314837/1‑148 (MQ=255)
                                                                                                       gACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACg                                         >  1:3451041/1‑148 (MQ=255)
                                                                                                        aCCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGt                                        >  2:1077142/1‑148 (MQ=255)
                                                                                                        aCCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGt                                        >  2:1077206/1‑148 (MQ=255)
                                                                                                          cAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTga                                      >  2:2744162/1‑148 (MQ=255)
                                                                                                            aaGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTgaag                                    <  2:1528229/148‑1 (MQ=255)
                                                                                                            aaGAGAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTgaag                                    <  2:135999/148‑1 (MQ=255)
                                                                                                             agagAACACCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTgaaga                                   <  1:3366794/148‑1 (MQ=255)
                                                                                                                    acCGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACg                            <  1:2320879/148‑1 (MQ=255)
                                                                                                                     ccGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACgg                           >  1:2813308/1‑148 (MQ=255)
                                                                                                                     ccGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACgg                           <  1:349895/148‑1 (MQ=255)
                                                                                                                     ccGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAAAATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACgg                           <  1:349852/148‑1 (MQ=255)
                                                                                                                      cGCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGt                          >  2:362495/1‑148 (MQ=255)
                                                                                                                       gCCCATTTCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACTTGAAGAATTTACGGTc                         >  2:2278206/1‑148 (MQ=255)
                                                                                                                            tttCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGa                    >  2:3469300/1‑148 (MQ=255)
                                                                                                                              tCAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATc                  <  1:3412169/148‑1 (MQ=255)
                                                                                                                               cAGCGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCa                 >  2:2059542/1‑148 (MQ=255)
                                                                                                                                  cGATCGAGACAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGt              <  1:1077206/148‑1 (MQ=255)
                                                                                                                                       gagaCAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCaaaa         <  1:3123610/148‑1 (MQ=255)
                                                                                                                                       gagaCAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCaaaa         <  1:2522669/148‑1 (MQ=255)
                                                                                                                                          aCAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATaa      >  2:1350556/1‑148 (MQ=255)
                                                                                                                                          aCAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATaa      >  2:2792181/1‑148 (MQ=255)
                                                                                                                                          aCAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATaa      >  2:3290126/1‑148 (MQ=255)
                                                                                                                                          aCAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGGGAAGAATTTACGGTCATAGATCACGTCAAAATaa      >  2:218663/1‑148 (MQ=255)
                                                                                                                                           cAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATaag     <  2:169444/148‑1 (MQ=255)
                                                                                                                                           cAGAGATAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATaag     <  1:1376214/148‑1 (MQ=255)
                                                                                                                                            agagatAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATaaga    >  2:1916792/1‑148 (MQ=255)
                                                                                                                                              agataATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAATAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATaagaag  <  1:1039221/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAG  >  CP009273/882527‑882816

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