Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,920,794 G→T intergenic (‑423/+132) RPPX_08560 ← / ← RPPX_08565 16S ribosomal RNA/CoA‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,920,7940GT98.6% 160.6 / ‑3.6 72intergenic (‑423/+132)RPPX_08560/RPPX_0856516S ribosomal RNA/CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (31/40);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTTGATCGTAGCCCTTACAGC‑‑A‑‑CCTTCCGAATCGTCTAACTCATTGAATCTCAAGGAGTTTGTCGTTCCGTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCA  >  CP009974/1920691‑1920941
                                                                                                           |                                                                                                                                                   
cTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTACCTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgag                                                                                                      <  1:2243879/151‑1 (MQ=11)
 ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga                                                                                                     >  2:4716883/1‑151 (MQ=2)
 ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga                                                                                                     <  1:3606830/151‑1 (MQ=2)
   ggAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGACGTGGGGCGCATTATACGGGGATTCGAGAgg                                                                                                   >  2:4451016/1‑151 (MQ=2)
   ggAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAgg                                                                                                   <  2:1026015/151‑1 (MQ=2)
   ggAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAgg                                                                                                   <  2:3060508/151‑1 (MQ=2)
    gAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTTGAGAggg                                                                                                  >  1:2175135/1‑151 (MQ=2)
    gAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATGATAAGGGGATTGGAGCggg                                                                                                  >  1:2175116/1‑151 (MQ=2)
      aGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCg                                                                                                >  1:974192/1‑151 (MQ=2)
      aGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTGTCGTTCCGTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCg                                                                                                <  1:1857849/151‑1 (MQ=11)
         gcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCa                                                                                             >  2:4412377/1‑151 (MQ=2)
         gcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCa                                                                                             <  2:1156168/151‑1 (MQ=2)
           tGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCGCTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAAc                                                                                           <  2:923169/151‑1 (MQ=2)
            gTCCACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAAcc                                                                                          <  2:1789953/151‑1 (MQ=2)
              cAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCtt                                                                                        <  2:1009439/151‑1 (MQ=2)
              cAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCtt                                                                                        >  1:4597702/1‑151 (MQ=2)
              cAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCtt                                                                                        >  2:4650512/1‑151 (MQ=2)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTATCTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                       <  2:3484699/151‑1 (MQ=2)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                       <  2:461927/151‑1 (MQ=11)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                       <  2:5227266/151‑1 (MQ=2)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                       <  1:2547308/151‑1 (MQ=2)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                       <  1:714560/151‑1 (MQ=2)
                 ccATCTTCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTaa                                                                                     <  1:4357038/151‑1 (MQ=2)
                 ccATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTaa                                                                                     >  1:972799/1‑151 (MQ=2)
                    tcatcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAttt                                                                                  <  2:2474804/151‑1 (MQ=2)
                    tcatcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAttt                                                                                  >  1:4022223/1‑151 (MQ=2)
                    tcatcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAttt                                                                                  <  2:3429108/151‑1 (MQ=2)
                     catcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTc                                                                                 >  1:5039975/1‑151 (MQ=2)
                     catcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTc                                                                                 <  1:415163/151‑1 (MQ=2)
                      atcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGCTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATCATAACGCGATTCGCTATGGCGTCAACCTTTAATTTCa                                                                                >  2:3751597/1‑151 (MQ=2)
                                                                           cactacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCt                               <  2:4660602/144‑1 (MQ=255)
                                                                            actacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATTTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTc                              >  2:4695901/7‑151 (MQ=21)
                                                                            actacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTTATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTc                              <  2:5231299/145‑1 (MQ=21)
                                                                             ctacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcc                             <  2:153164/146‑1 (MQ=255)
                                                                             ctacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcc                             >  1:4518252/6‑151 (MQ=255)
                                                                             ctacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcc                             >  1:5067171/6‑151 (MQ=255)
                                                                              cactttaCTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGAGATTCGAGAGGGCGTCAAGCTTTAATTTCAATAACCGTTGATAGCACAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTccc                            <  2:3512838/145‑1 (MQ=21)
                                                                               ccttAACTCATTGAATCTCAAGGAGTATTTCATACCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcccc                           <  2:2525626/148‑1 (MQ=21)
                                                                               acttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcccc                           <  1:4463541/148‑1 (MQ=255)
                                                                               acttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcccc                           >  2:2726395/4‑151 (MQ=255)
                                                                               acgcaaCTCATTGAATCACAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGCGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcccc                           <  2:4819392/147‑1 (MQ=21)
                                                                               accTAACTCATGGAATCTCAAGGAGTTTTTCGTTCCTATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcccc                           <  2:2829402/149‑1 (MQ=255)
                                                                                 tgaaCTCATTGAATCTCAAGGAGTTTTTCGTTCCGGTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGAGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGc                         <  2:3260942/149‑1 (MQ=255)
                                                                                    aCTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTAt                      >  1:4907401/1‑151 (MQ=255)
                                                                                    aCTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTAt                      >  1:3687562/1‑151 (MQ=255)
                                                                                    aCTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTAt                      <  2:4613580/151‑1 (MQ=255)
                                                                                     cTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTAtt                     >  1:1265222/1‑151 (MQ=255)
                                                                                      tCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTc                    >  1:4615079/1‑151 (MQ=255)
                                                                                         ttGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATc                 >  1:2061870/1‑151 (MQ=255)
                                                                                         ttGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATc                 >  2:588341/1‑151 (MQ=255)
                                                                                          tGAATGTCAAGGAGTTTTTCGTTCCGATGTCGCGGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTGATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCt                <  2:5317872/151‑1 (MQ=255)
                                                                                          tGAATCTCAAGGAGTGTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCt                <  2:2399366/151‑1 (MQ=255)
                                                                                           gAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTAATAAAAAAAAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTa               <  2:3928027/151‑1 (MQ=18)
                                                                                            aaTCTCAAGGAGTTTTGCGTTCCGAGGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTag              <  2:3250990/151‑1 (MQ=255)
                                                                                             aTCTCGAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTaga             <  1:265533/151‑1 (MQ=255)
                                                                                             aTCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTaga             <  1:1090279/151‑1 (MQ=255)
                                                                                             aTCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTaga             <  2:1403186/151‑1 (MQ=255)
                                                                                             aTCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTaga             >  1:2186791/1‑151 (MQ=255)
                                                                                              tctcAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACACAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTagag            >  2:3692263/1‑151 (MQ=255)
                                                                                                tcAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGca          >  1:3132233/1‑151 (MQ=255)
                                                                                                tcAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGca          <  1:3729773/151‑1 (MQ=255)
                                                                                                 cAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGcac         >  2:1051702/1‑151 (MQ=255)
                                                                                                 cAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGcac         >  1:2373197/1‑151 (MQ=255)
                                                                                                  aaGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCt                <  2:4397733/143‑1 (MQ=255)
                                                                                                  aaGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCt                >  1:4397733/1‑143 (MQ=255)
                                                                                                  aaGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGcaca        >  1:3180844/1‑151 (MQ=255)
                                                                                                    ggAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACAcc      <  1:1240480/151‑1 (MQ=255)
                                                                                                     gAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCa     >  1:2711864/1‑151 (MQ=255)
                                                                                                     gAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGAAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCa     <  2:619929/151‑1 (MQ=255)
                                                                                                      aGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAg    <  1:2092386/151‑1 (MQ=255)
                                                                                                      aGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAg    <  2:4487869/151‑1 (MQ=255)
                                                                                                       gTTTTTCGTTCCGATATCGCTGGAAGTGGGGAGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCATTGAGAGCGCAGAAAAAACAAAGCGACGAGGCCTGACGGCCTCCCCGCTATCCATCTAGAGCACAccctc   >  2:3804920/1‑148 (MQ=25)
                                                                                                        ttttTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGca  <  2:3032650/151‑1 (MQ=255)
                                                                                                           |                                                                                                                                                   
CTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTTGATCGTAGCCCTTACAGC‑‑A‑‑CCTTCCGAATCGTCTAACTCATTGAATCTCAAGGAGTTTGTCGTTCCGTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCA  >  CP009974/1920691‑1920941

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