Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,404,772 C→A 88.6% A20A (GCG→GCT yfbS ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,404,7720CA88.6% 428.6 / 53.3 176A20A (GCG→GCTyfbSputative transporter
Reads supporting (aligned to +/- strand):  ref base C (7/13);  new base A (86/70);  total (93/83)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 176 total aligned reads displayed.

CATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTT  >  CP009273/2404627‑2404918
                                                                                                                                                 |                                                                                                                                                  
cATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGc                                                                                                                                                  >  1:2825547/1‑148 (MQ=255)
cATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGc                                                                                                                                                  >  2:769364/1‑148 (MQ=255)
  tCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCaa                                                                                                                                                >  2:3204226/1‑148 (MQ=255)
   cGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCaaa                                                                                                                                               >  1:1887212/1‑148 (MQ=255)
     ccAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACa                                                                                                                                             >  1:2691584/1‑148 (MQ=255)
      cAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACaa                                                                                                                                            >  1:2883953/1‑148 (MQ=255)
      cAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACaa                                                                                                                                            >  1:3205116/1‑148 (MQ=255)
       aataatAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAg                                                                                                                                           >  1:1811469/1‑148 (MQ=255)
       aataaaAAACACGGGGTAAATCATCAAAACTTTTGTATAATAAAACCCTTAAAATACTTCTGTGACTGTCAGCTTTCCGCTTAATCAAAACGCTAAAATAAAAAAAAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAg                                                                                                                                           <  2:2406314/148‑1 (MQ=255)
             aaaCAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATg                                                                                                                                     >  2:2475893/1‑148 (MQ=255)
               aCATGGCGTCAATCTGGTCAACTTTCGTAACAAAAACGCCTTAAAAAAGTTGTGGCACTGTCAGGGTTCCGATTAATCCAACCGCTAAAATAACAAACACCGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGc                                                                                                                                   <  2:1479139/148‑1 (MQ=255)
                catggcgtcAATCTGGACAACGTTAGTATCAGAAAAGCGGGAAAATAGTTCTGGGAGTGTCAGGGTTCCGCTTAATGCAAACGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCa                                                                                                                                  <  2:2729967/145‑1 (MQ=255)
                cacggcggcTATCATTACAACTTTTGGATCATAAATGCCGGAAAATAATTCTGGGAGTGTCTGCTTTCCGCTTATTTCAAACGTTAAAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCa                                                                                                                                  <  2:2730135/145‑1 (MQ=255)
                     cgccAATCAGGACAACTTTAGTATCATAAAAGCCGCAAAATAGTTCTGGGAGTGTCAGCGTTCCGATTAATCAAAACGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGc                                                                                                                             <  2:2108704/145‑1 (MQ=255)
                       gcAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCa                                                                                                                           >  1:2745347/1‑148 (MQ=255)
                               gACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGaa                                                                                                                   >  1:1393309/1‑148 (MQ=255)
                               gACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGaa                                                                                                                   >  2:2051858/1‑148 (MQ=255)
                                 cAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAg                                                                                                                 >  1:3260129/1‑148 (MQ=255)
                                    ctttATGATAAGAAAACCCGTAAAATACTTCTGGGACTTTCAGCTTTCCGATTAATGCAAACGCTAAAATAAAAAACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAAc                                                                                                              <  2:1647081/146‑1 (MQ=255)
                                    cGTTTGTATCAGAAAAGCCGTAAAATACTTCTGGTACTGTCAGCGTTCCGATTAATGCAAACGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAAc                                                                                                              <  2:1647085/148‑1 (MQ=255)
                                      ttAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAAccc                                                                                                            >  1:1984208/1‑148 (MQ=255)
                                       tAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCa                                                                                                           >  2:900231/1‑148 (MQ=255)
                                       tAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCa                                                                                                           >  1:1028275/1‑148 (MQ=255)
                                       tAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCa                                                                                                           >  2:156470/1‑148 (MQ=255)
                                          tataaGAAAAGCCGTAAAATACTTGTGGAAGTTTCAGGTTTCCGATTAATGCAAACGCTAAAATAAAAAAAACCGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAAt                                                                                                        <  2:1998536/144‑1 (MQ=255)
                                          gATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAAt                                                                                                        >  1:1747886/1‑148 (MQ=255)
                                          gATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAAt                                                                                                        >  1:1747631/1‑148 (MQ=255)
                                           aTCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATc                                                                                                       >  1:3113334/1‑148 (MQ=255)
                                           aTCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATc                                                                                                       >  2:1713817/1‑148 (MQ=255)
                                             cAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCaa                                                                                                     >  1:3061541/1‑148 (MQ=255)
                                                  aaGCCGTAAAATACTTCTGGGACTTTCAGCGTTCCGCTTAATGAAAACGCGACAATAAAAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCAc                                                                                                <  2:2690414/148‑1 (MQ=255)
                                                   aGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCAcc                                                                                               >  1:1144363/1‑148 (MQ=255)
                                                     ccGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGt                                                                                             >  2:1955237/1‑148 (MQ=255)
                                                         aaaaTACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTcac                                                                                         >  1:821767/1‑148 (MQ=255)
                                                          aaaTACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACACGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACAGATcaca                                                                                        >  1:3062733/1‑148 (MQ=255)
                                                           aaTACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACaa                                                                                       >  1:1732786/1‑148 (MQ=255)
                                                           aaTACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAAGTCACCGTTCACaa                                                                                       >  1:1884537/1‑148 (MQ=255)
                                                               cTTCTGGGACTGACAGCGTACCGATTAATCCAACCGCGACAAAAAAAAACACCGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTa                                                                                   <  2:206137/148‑1 (MQ=255)
                                                                ttCTGGTATTTTCTGCGTTCCGATTAATGAAAACGCTAAAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTAt                                                                                  <  2:2189823/148‑1 (MQ=255)
                                                                     ggACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTg                                                                             >  2:2650087/1‑148 (MQ=255)
                                                                            cAGCGTTCCGCTTAATCCAACCGCTAAAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTg                                                                      <  2:2880015/148‑1 (MQ=255)
                                                                            cAGCGTTCCGCTTAATCAAAACGCCAAAATAAAAAACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTg                                                                      <  2:274101/148‑1 (MQ=255)
                                                                              gCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAg                                                                    >  1:1936465/1‑148 (MQ=255)
                                                                              gCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAg                                                                    >  2:1256218/1‑148 (MQ=255)
                                                                                 ttCCGATAAATCAAAACGAGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAgggg                                                                 <  2:2359831/148‑1 (MQ=255)
                                                                                  tccgaaaaaTCCAAACCCAAAAAAAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAggggg                                                                <  1:2572966/141‑1 (MQ=255)
                                                                                  tCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAggggg                                                                >  1:747996/1‑148 (MQ=255)
                                                                                  tCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTAACCGTTCACAACTTATCATTGTTAATTGAggggg                                                                >  1:748001/1‑148 (MQ=255)
                                                                                    ctattAATGAAAACGCGAAAATAAAAAACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGAt                                                              <  2:1897415/145‑1 (MQ=255)
                                                                                       ttAATGAAAACGCCAAAATAAAAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGAc                                                           <  2:443178/148‑1 (MQ=255)
                                                                                       ttAATCCAAACGCCACAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGAc                                                           <  2:442640/148‑1 (MQ=255)
                                                                                       ttAATCCAAACCCTAAAAAAACAAACCCCGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGAc                                                           <  1:451569/148‑1 (MQ=255)
                                                                                            ccacaCGCCACAAAAAAACACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGAt                                                      <  1:1647529/144‑1 (MQ=255)
                                                                                             cAACCGCTAAAATAACAACCAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGAtt                                                     <  2:2706235/148‑1 (MQ=255)
                                                                                             cAAACGCGACAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGAtt                                                     <  2:2706231/148‑1 (MQ=255)
                                                                                              aaaCGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTc                                                    >  2:1096083/1‑148 (MQ=255)
                                                                                                ccgcaaaaaTAAAAAACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTg                                                  <  1:2206744/141‑1 (MQ=255)
                                                                                                 cgcgACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGc                                                 >  1:1627542/1‑148 (MQ=255)
                                                                                                 cgccACAAAAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATAAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGc                                                 <  1:2463905/148‑1 (MQ=255)
                                                                                                  gctaaaaTACCAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGcc                                                <  2:2745347/143‑1 (MQ=255)
                                                                                                    caaaaTAACAACCAACGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCAt                                              <  2:2621760/145‑1 (MQ=255)
                                                                                                      cAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATaa                                            >  2:1544321/1‑148 (MQ=255)
                                                                                                      cAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATaa                                            >  2:2601027/1‑148 (MQ=255)
                                                                                                       aaTAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATaaa                                           >  2:2923073/1‑148 (MQ=255)
                                                                                                          aacaaacaAAGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATaaaaaa                                        <  2:1393309/148‑1 (MQ=255)
                                                                                                           aaaaacaacGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAg                                       <  1:2869366/146‑1 (MQ=255)
                                                                                                            aaaacaaaGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGc                                      <  1:2392651/147‑1 (MQ=255)
                                                                                                              accaccgCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGccc                                    <  1:347246/142‑1 (MQ=255)
                                                                                                              acaaacgCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGccc                                    <  2:2810149/142‑1 (MQ=255)
                                                                                                              aaaaacgCGACCGCATCCATACGCTCTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCATTTCACCGTTCACAACTTATCCTTGTTAATTGTGGGGGATGACTTGATTCTGCCATAAAAAAGccc                                    <  2:1126142/142‑1 (MQ=255)
                                                                                                              aaaaaaGCTACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAATAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGccc                                    <  2:2809053/145‑1 (MQ=255)
                                                                                                               acaaaGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGcccc                                   >  1:2053679/1‑148 (MQ=255)
                                                                                                                  aaGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCcgac                                >  1:1535451/1‑148 (MQ=255)
                                                                                                                   cgCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCcgacg                               <  1:545954/147‑1 (MQ=255)
                                                                                                                   aGCGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCcgacg                               >  1:1399521/1‑148 (MQ=255)
                                                                                                                     cGACCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACgag                             >  1:1049193/1‑148 (MQ=255)
                                                                                                                       aCCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGAc                           >  1:3268770/1‑148 (MQ=255)
                                                                                                                       aCCGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGAc                           >  1:1631160/1‑148 (MQ=255)
                                                                                                                        ccGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACg                          >  1:1213753/1‑148 (MQ=255)
                                                                                                                        ccGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACg                          >  1:928204/1‑148 (MQ=255)
                                                                                                                         cGCATCCATACGCACTCTGCCCGTCGCAAACAAGCCGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgg                         <  1:2126913/148‑1 (MQ=255)
                                                                                                                         cGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgg                         <  1:2946519/148‑1 (MQ=255)
                                                                                                                         cGCATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgg                         <  1:2946010/148‑1 (MQ=255)
                                                                                                                           cATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgggg                       <  2:1936465/148‑1 (MQ=255)
                                                                                                                           cATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgggg                       >  2:1468243/1‑148 (MQ=255)
                                                                                                                           cATCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACgggg                       <  1:3194605/148‑1 (MQ=255)
                                                                                                                            aTCCATACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCACCTCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGc                      <  2:3244893/148‑1 (MQ=255)
                                                                                                                               cATACGCACTCTGCCCGTAGCAAACAAGACGCTGGAAACGGCCAGTAATGAAAGAACCCACATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTaa                   <  1:2591215/148‑1 (MQ=255)
                                                                                                                               cATACGCACTCTGCCCGTACCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTaa                   <  1:317251/148‑1 (MQ=255)
                                                                                                                                 tACGCACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAAt                 <  2:3205116/148‑1 (MQ=255)
                                                                                                                                     cACTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATg             <  2:299056/148‑1 (MQ=255)
                                                                                                                                      aCTCTGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGa            <  2:397064/148‑1 (MQ=255)
                                                                                                                                        tctGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATc          <  2:550022/148‑1 (MQ=255)
                                                                                                                                         ctGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa         <  1:156470/148‑1 (MQ=255)
                                                                                                                                         ctGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa         >  2:26478/1‑148 (MQ=255)
                                                                                                                                         ctGCCCGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCa         <  2:2657060/148‑1 (MQ=255)
                                                                                                                                            cccGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGgt      <  2:3091573/148‑1 (MQ=255)
                                                                                                                                             ccGTAGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGgtg     <  2:1728187/148‑1 (MQ=255)
                                                                                                                                               gTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGtt   <  1:900231/148‑1 (MQ=255)
                                                                                                                                                tCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGttt  <  2:1146494/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCTGCCATAAAAAAGCCCCGACGAGACGGGGCTAAATCATGATCAGGTGTTT  >  CP009273/2404627‑2404918

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