Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,920,803 T→A intergenic (‑432/+123) 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,8030TA98.9% 210.7 / NA 89intergenic (‑432/+123)RPPX_08560/RPPX_0856516S ribosomal RNA/CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (39/47);  minor base G (0/1);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTTGATCGTAGCCCTTACAGC‑‑A‑‑CCTTCCGAATCGTCTAACTCATTGAATCTCAAGGAGTTTGTCGTTCCGTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAGG  >  CP009974/1920691‑1920950
                                                                                                                    |                                                                                                                                                   
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:5227266/151‑1 (MQ=2)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                                <  1:2547308/151‑1 (MQ=2)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                                <  1:714560/151‑1 (MQ=2)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                                <  2:461927/151‑1 (MQ=11)
                 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:415163/151‑1 (MQ=2)
                     catcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTc                                                                                          >  1:5039975/1‑151 (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                                      >  1:4518252/6‑151 (MQ=255)
                                                                             ctacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcc                                      >  1:5067171/6‑151 (MQ=255)
                                                                             ctacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcc                                      <  2:153164/146‑1 (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                               <  2:4613580/151‑1 (MQ=255)
                                                                                    aCTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTAt                               >  1:3687562/1‑151 (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                      <  2:1403186/151‑1 (MQ=255)
                                                                                             aTCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTaga                      >  1:2186791/1‑151 (MQ=255)
                                                                                             aTCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTaga                      <  1:1090279/151‑1 (MQ=255)
                                                                                              tctcAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACACAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTagag                     >  2:3692263/1‑151 (MQ=255)
                                                                                                tcAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGca                   <  1:3729773/151‑1 (MQ=255)
                                                                                                tcAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGca                   >  1:3132233/1‑151 (MQ=255)
                                                                                                 cAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGcac                  >  2:1051702/1‑151 (MQ=255)
                                                                                                 cAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGcac                  >  1:2373197/1‑151 (MQ=255)
                                                                                                  aaGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCt                         >  1:4397733/1‑143 (MQ=255)
                                                                                                  aaGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCt                         <  2:4397733/143‑1 (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           <  1:2474201/147‑1 (MQ=255)
                                                                                                        ttttTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGca           <  2:3032650/151‑1 (MQ=255)
                                                                                                         ttttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcac          >  1:428419/4‑151 (MQ=255)
                                                                                                           ttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcacca        >  1:5007053/2‑151 (MQ=255)
                                                                                                           ttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcacca        <  2:5120121/150‑1 (MQ=255)
                                                                                                            tCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGACCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcaccac       <  1:4077741/151‑1 (MQ=255)
                                                                                                            gcGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcaccac       <  2:4237004/150‑1 (MQ=255)
                                                                                                             cGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACg      >  1:2545249/1‑151 (MQ=255)
                                                                                                              gTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGa     >  1:2590376/1‑151 (MQ=255)
                                                                                                              gTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCCCCCCGCCATTCATCTACAGCACACCACCACCACGa     >  2:5091406/1‑151 (MQ=255)
                                                                                                               ttCCGCTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGaa    <  1:1741399/151‑1 (MQ=255)
                                                                                                               ttCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGaa    <  1:4567049/151‑1 (MQ=255)
                                                                                                                tCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAg   <  1:4300495/151‑1 (MQ=255)
                                                                                                                tCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAg   >  2:436139/1‑151 (MQ=255)
                                                                                                                tCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAg   >  1:781316/1‑151 (MQ=255)
                                                                                                                 ccGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAgg  <  2:3006769/151‑1 (MQ=255)
                                                                                                                 ccGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATACACCTAGAGCACACCAGCACCACGAAgg  >  1:1258328/1‑151 (MQ=255)
                                                                                                                    |                                                                                                                                                   
CTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTTGATCGTAGCCCTTACAGC‑‑A‑‑CCTTCCGAATCGTCTAACTCATTGAATCTCAAGGAGTTTGTCGTTCCGTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAGG  >  CP009974/1920691‑1920950

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