Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,901,415 C→A 100% P104P (CCG→CCT pstB ← phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,901,4150CA96.9% 607.0 / 9.0 195P104P (CCG→CCTpstBphosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base A (107/82);  total (111/84)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 195 total aligned reads displayed.

AGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTG  >  CP009273/3901270‑3901559
                                                                                                                                                 |                                                                                                                                                
 gTAACCGCTCTGGTGCAATTTCTCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGa                                                                                                                                               <  2:1296734/148‑1 (MQ=255)
  tAACCGCTCTGGTGCAATTTCTCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGaa                                                                                                                                              <  2:3225837/148‑1 (MQ=255)
  tAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGaa                                                                                                                                              >  2:3421701/1‑148 (MQ=255)
   aaCCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGaaa                                                                                                                                             >  2:3556184/1‑148 (MQ=255)
   aaCCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGaaa                                                                                                                                             >  1:2681520/1‑148 (MQ=255)
   aaCCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGaaa                                                                                                                                             >  1:1496862/1‑148 (MQ=255)
    aCCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAAc                                                                                                                                            >  1:1390312/1‑148 (MQ=255)
    aCCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAAc                                                                                                                                            <  1:1682402/148‑1 (MQ=255)
     ccGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACg                                                                                                                                           >  1:603867/1‑148 (MQ=255)
         tctGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGt                                                                                                                                       <  2:2860984/148‑1 (MQ=255)
         tctGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGt                                                                                                                                       <  1:566670/148‑1 (MQ=255)
         tctGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGt                                                                                                                                       <  2:1085619/148‑1 (MQ=255)
         tctGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGt                                                                                                                                       <  2:2664273/148‑1 (MQ=255)
            ggTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCgg                                                                                                                                    <  1:2795517/148‑1 (MQ=255)
                cAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTc                                                                                                                                >  2:212839/1‑148 (MQ=255)
                cAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTc                                                                                                                                <  1:967042/148‑1 (MQ=255)
                     tATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGaa                                                                                                                           >  1:2561379/1‑148 (MQ=255)
                       tCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAaca                                                                                                                         >  2:2582277/1‑148 (MQ=255)
                             gTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGc                                                                                                                   <  2:1751213/148‑1 (MQ=255)
                              tttCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGcc                                                                                                                  <  1:1246011/148‑1 (MQ=255)
                                tCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCa                                                                                                                >  1:870123/1‑148 (MQ=255)
                                   ttCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACtt                                                                                                             <  1:241866/148‑1 (MQ=255)
                                   ttCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACtt                                                                                                             <  2:1246621/148‑1 (MQ=255)
                                    tCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACttt                                                                                                            >  1:2082430/1‑148 (MQ=255)
                                        caATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTcgca                                                                                                        >  1:1329486/1‑148 (MQ=255)
                                         aaTGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTcgcac                                                                                                       >  1:1636041/1‑148 (MQ=255)
                                           tGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTcgcacgc                                                                                                     >  1:1917018/1‑148 (MQ=255)
                                           tGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTcgcacgc                                                                                                     >  1:952246/1‑148 (MQ=255)
                                           tGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTcgcacgc                                                                                                     <  2:3218819/148‑1 (MQ=255)
                                             cGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACgcag                                                                                                   >  2:3456609/1‑148 (MQ=255)
                                                 tttGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGt                                                                                               <  2:1151805/148‑1 (MQ=255)
                                                 tttGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGt                                                                                               <  2:2925866/148‑1 (MQ=255)
                                                  ttGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTg                                                                                              >  1:2860014/1‑148 (MQ=255)
                                                    ggTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCg                                                                                            >  1:3087986/1‑148 (MQ=255)
                                                     gTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGa                                                                                           >  2:477178/1‑148 (MQ=255)
                                                     gTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGa                                                                                           >  2:2734087/1‑148 (MQ=255)
                                                       cAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGata                                                                                         >  1:1131616/1‑148 (MQ=255)
                                                       cAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGata                                                                                         <  2:1812256/148‑1 (MQ=255)
                                                         aTGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATc                                                                                       <  2:2235912/148‑1 (MQ=255)
                                                          tGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATcc                                                                                      <  2:1579968/148‑1 (MQ=255)
                                                           gCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCt                                                                                     >  2:65121/1‑148 (MQ=255)
                                                            cccACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTg                                                                                    <  2:1558232/148‑1 (MQ=255)
                                                            cccACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTg                                                                                    >  2:3290303/1‑148 (MQ=255)
                                                            cccACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTg                                                                                    >  1:1566489/1‑148 (MQ=255)
                                                              cACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTgag                                                                                  >  1:3311334/1‑148 (MQ=255)
                                                                cTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGagag                                                                                <  2:587022/148‑1 (MQ=255)
                                                                cTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGagag                                                                                >  2:3018460/1‑148 (MQ=255)
                                                                    aCGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTgg                                                                            <  1:2964986/148‑1 (MQ=255)
                                                                     cgcgCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGt                                                                           >  1:3134767/1‑148 (MQ=255)
                                                                     cgcgCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGt                                                                           <  1:1550384/148‑1 (MQ=255)
                                                                         cTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAgg                                                                       >  1:2397587/1‑148 (MQ=255)
                                                                            gTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGAtg                                                                    <  2:1308439/148‑1 (MQ=255)
                                                                               cATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGAtgttg                                                                 <  1:1602734/148‑1 (MQ=255)
                                                                                 tGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTc                                                               >  1:1001261/1‑148 (MQ=255)
                                                                                    cGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGcc                                                            >  1:665309/1‑148 (MQ=255)
                                                                                        aCGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCa                                                        <  1:1491258/148‑1 (MQ=255)
                                                                                              gagCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGa                                                  <  2:3312022/148‑1 (MQ=255)
                                                                                               agCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGaa                                                 <  2:2195480/148‑1 (MQ=255)
                                                                                               agCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGaa                                                 >  1:1751470/1‑148 (MQ=255)
                                                                                                gCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAt                                                <  2:869298/148‑1 (MQ=255)
                                                                                                 cTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAtt                                               >  2:1923058/1‑148 (MQ=255)
                                                                                                      cAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAcc                                          <  2:1580922/148‑1 (MQ=255)
                                                                                                      cAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAcc                                          <  2:3343039/148‑1 (MQ=255)
                                                                                                      cAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAcc                                          <  1:3199212/148‑1 (MQ=255)
                                                                                                      cAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAcc                                          <  2:1983799/148‑1 (MQ=255)
                                                                                                       aaaCAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                         <  1:1794758/148‑1 (MQ=255)
                                                                                                          cAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTcc                                      >  1:2046865/1‑148 (MQ=255)
                                                                                                           aGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCg                                     >  1:2292496/1‑148 (MQ=255)
                                                                                                            gACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGc                                    >  2:1806720/1‑148 (MQ=255)
                                                                                                             aCGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCa                                   >  1:3529743/1‑148 (MQ=255)
                                                                                                              cGAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCAc                                  >  1:1417032/1‑148 (MQ=255)
                                                                                                               gAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACg                                 <  2:1360525/148‑1 (MQ=255)
                                                                                                               gAACGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACg                                 <  2:2790466/148‑1 (MQ=255)
                                                                                                                aaCGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgc                                >  2:2977254/1‑148 (MQ=255)
                                                                                                                 aCGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgct                               >  2:2224055/1‑148 (MQ=255)
                                                                                                                  cGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctg                              >  1:2334702/1‑148 (MQ=255)
                                                                                                                  cGCCAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctg                              >  2:2374517/1‑148 (MQ=255)
                                                                                                                    ccAAAAGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctgct                            <  1:2377975/148‑1 (MQ=255)
                                                                                                                     cAAAATCGATGATGTCGTAGATGGACATAGGAAACTGCGTCGGTTTCTGGAAAAACATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCAACAAGCACAATATCACCATCAGCACgcagatc                           >  2:2930945/1‑143 (MQ=255)
                                                                                                                       aaaGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCg                         >  2:2197377/1‑148 (MQ=255)
                                                                                                                        aaGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCgg                        <  1:45328/148‑1 (MQ=255)
                                                                                                                        aaGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCgg                        >  2:3498444/1‑148 (MQ=255)
                                                                                                                         aGCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCggg                       <  2:1448740/148‑1 (MQ=255)
                                                                                                                          gCGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGt                      >  2:202909/1‑148 (MQ=255)
                                                                                                                           cGATGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTa                     <  1:2224708/148‑1 (MQ=255)
                                                                                                                             aTGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                   >  1:987769/1‑148 (MQ=255)
                                                                                                                             aTGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                   <  1:2013052/148‑1 (MQ=255)
                                                                                                                             aTGTTGTCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                   >  2:670160/1‑148 (MQ=255)
                                                                                                                              tgttgtCGTAGATGGACATAGGAAACGGCGTCTGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAg                  <  1:3529862/148‑1 (MQ=255)
                                                                                                                              tgttgtCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAg                  >  1:1232975/1‑148 (MQ=255)
                                                                                                                              tgttgtCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAg                  <  2:2187188/148‑1 (MQ=255)
                                                                                                                                   tCGTAGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCaa             <  2:442966/148‑1 (MQ=255)
                                                                                                                                       aGATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAt         <  1:2828503/148‑1 (MQ=255)
                                                                                                                                        gATGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAtt        >  2:110313/1‑148 (MQ=255)
                                                                                                                                         aTGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt       >  1:2358063/1‑148 (MQ=255)
                                                                                                                                         aTGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt       <  1:2831835/148‑1 (MQ=255)
                                                                                                                                         aTGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt       >  2:107020/1‑148 (MQ=255)
                                                                                                                                          tGGACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAtttt      >  1:2206058/1‑148 (MQ=255)
                                                                                                                                           ggACATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttg     <  2:3064094/148‑1 (MQ=255)
                                                                                                                                              cATAGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgttg  <  1:3357011/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTG  >  CP009273/3901270‑3901559

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: