Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,216,180 C→G 7.4% A3G (GCC→GGC)  PP_5493 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,216,1800CG 7.4% 212.7 / 7.6 68A3G (GCC→GGC) PP_5493hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/36);  new base G (3/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GGATCATTCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGAGGGCGCCTATGAAGTG  >  NC_002947/2216034‑2216328
                                                                                                                                                   |                                                                                                                                                    
ggATCATTCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGcgcgccgc                                                                                                                                                   <  2:367424/150‑1 (MQ=255)
  aTCATTCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAAATTAAGTGCAGAAGAAATGtgcgccgcgc                                                                                                                                                 <  2:1731603/150‑1 (MQ=255)
      ttCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCtttt                                                                                                                                             >  1:1041513/1‑150 (MQ=255)
      ttCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCtttt                                                                                                                                             >  2:955923/1‑150 (MQ=255)
      ttCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCtttt                                                                                                                                             <  1:2440615/150‑1 (MQ=255)
      ttCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCtttt                                                                                                                                             >  1:2265414/1‑150 (MQ=255)
      ttCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCtttt                                                                                                                                             >  1:1337289/1‑150 (MQ=255)
       tCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTc                                                                                                                                            <  2:1002492/150‑1 (MQ=255)
       tCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTc                                                                                                                                            <  1:1820427/150‑1 (MQ=255)
         gCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCgg                                                                                                                                          >  2:1374135/1‑150 (MQ=255)
            ttAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGt                                                                                                                                       <  2:1220497/150‑1 (MQ=255)
              aGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTgg                                                                                                                                     <  1:1321469/150‑1 (MQ=255)
              aGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTgg                                                                                                                                     <  2:2050161/150‑1 (MQ=255)
                  ggTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTaa                                                                                                                                 <  2:1811250/150‑1 (MQ=255)
                      tctcACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGCGGTTAATGtt                                                                                                                             >  1:1285108/1‑150 (MQ=255)
                        tcACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTAt                                                                                                                           >  1:1650321/1‑150 (MQ=255)
                          acacGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATcc                                                                                                                         <  2:2451482/150‑1 (MQ=255)
                          acacGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATcc                                                                                                                         >  2:1896516/1‑150 (MQ=255)
                            acGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCt                                                                                                                       <  2:2265414/150‑1 (MQ=255)
                             aaataGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGGATGCATCTGAACTGATGTGCG‑TCGCCTGAGCAGAAGACTTCACAAGAGAATAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTc                                                                                                                      >  1:543481/5‑150 (MQ=255)
                                 aGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGa                                                                                                                  <  2:1337289/150‑1 (MQ=255)
                                  gTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGt                                                                                                                                       >  2:116865/1‑128 (MQ=255)
                                  gTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGt                                                                                                                                       <  1:116865/128‑1 (MQ=255)
                                      aTTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACt                                                                                                             <  1:1217627/150‑1 (MQ=255)
                                        tCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTcc                                                                                                           >  2:2013169/1‑150 (MQ=255)
                                               tttCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCt                                                                                                    <  1:955923/150‑1 (MQ=255)
                                                    cAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCAATCAAGAAGATACCCTGCGAACATCt                                                                                               >  2:1134075/1‑150 (MQ=255)
                                                       tGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGc                                                                                            <  2:876696/150‑1 (MQ=255)
                                                                    cgaatgcgaTCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAgg                                                                               <  2:543481/150‑1 (MQ=255)
                                                                    cgaatgcgaTCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAgg                                                                               <  2:1041513/150‑1 (MQ=255)
                                                                      aatgcgaTCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCa                                                                             >  2:2192935/1‑150 (MQ=255)
                                                                         gcgaTCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTa                                                                          <  1:2024065/150‑1 (MQ=255)
                                                                                     gCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCg                                                              <  2:1681564/150‑1 (MQ=255)
                                                                                      cATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTACGTGCCGAAGAAATTCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCgg                                                             <  2:1243962/150‑1 (MQ=255)
                                                                                       atatCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCggg                                                            <  2:228387/150‑1 (MQ=255)
                                                                                          tCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCt                                                                                               <  1:773387/112‑1 (MQ=255)
                                                                                          tCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCt                                                                                               >  2:773387/1‑112 (MQ=255)
                                                                                            aaCTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACaa                                                       >  1:1317014/1‑150 (MQ=255)
                                                                                            aaCTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACaa                                                       >  1:1732604/1‑150 (MQ=255)
                                                                                            aaCTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGCTATCCCTCAAGAAAACTCCCCGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACaa                                                       <  2:1324310/150‑1 (MQ=255)
                                                                                             aCTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAgg                                                                               >  2:2178870/1‑125 (MQ=255)
                                                                                             aCTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAgg                                                                               <  1:2178870/125‑1 (MQ=255)
                                                                                               tGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACaaaaa                                                    <  1:2192935/150‑1 (MQ=255)
                                                                                                    tgAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTaaa                                               <  1:1896516/150‑1 (MQ=255)
                                                                                                     gAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAg                                              <  1:1602897/150‑1 (MQ=255)
                                                                                                       gTTTGATTCATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGtt                                            <  2:1317014/150‑1 (MQ=255)
                                                                                                        tttGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTaaa                                               <  1:207745/146‑1 (MQ=255)
                                                                                                        tttGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTaaa                                               >  2:207745/1‑146 (MQ=255)
                                                                                                        tttGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTaa                                                >  1:1031222/1‑145 (MQ=255)
                                                                                                        tttGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTaa                                                <  2:1031222/145‑1 (MQ=255)
                                                                                                                    gCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCggg                                                            <  1:1266660/121‑1 (MQ=255)
                                                                                                                    gCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCggg                                                            >  2:1266660/1‑121 (MQ=255)
                                                                                                                    gCA‑CCAACTTAAGTGAAGAAGAAAAGCGCGGCGCGCATTTc                                                                                                                                            <  1:232020/41‑1 (MQ=17)
                                                                                                                    gCA‑CCAACTTAAGTGAAGAAGAAAAGCGCGGCGCGCATTTc                                                                                                                                            >  2:232020/1‑41 (MQ=17)
                                                                                                                     cA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAgg                              >  1:1494493/1‑150 (MQ=255)
                                                                                                                        ccAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCa                            >  2:938722/1‑150 (MQ=255)
                                                                                                                                  tGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTATGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGa                  >  1:1455329/1‑150 (MQ=255)
                                                                                                                                   gaagaagaaATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGAg                 <  1:1703852/150‑1 (MQ=255)
                                                                                                                                   gaagaagaaATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTCGAg                 <  2:1494493/150‑1 (MQ=255)
                                                                                                                                    aagaagaaATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGAgg                <  2:1821934/150‑1 (MQ=255)
                                                                                                                                       aagaaATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGAGGgcg             >  1:726420/1‑150 (MQ=255)
                                                                                                                                       aagaaATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGAGGgcg             <  1:829003/150‑1 (MQ=255)
                                                                                                                                        cgaaAAGCGCGGCGCGCATTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCg                                                              <  1:811896/99‑1 (MQ=255)
                                                                                                                                        cgaaAAGCGCGGCGCGCATTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCg                                                              >  2:811896/2‑100 (MQ=255)
                                                                                                                                        cgaaAAGCGCGGCGCGCATTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGAGGgcgc            >  2:12200/2‑150 (MQ=255)
                                                                                                                                         gaaATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCtt                                                                                              <  2:2430273/67‑1 (MQ=255)
                                                                                                                                         gaaATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCtt                                                                                              >  1:2430273/1‑67 (MQ=255)
                                                                                                                                         gaaATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGAGGGCGcc           >  2:14542/1‑150 (MQ=255)
                                                                                                                                          aaaTGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTaaa                                               >  1:927893/1‑113 (MQ=255)
                                                                                                                                          aaaTGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTaaa                                               <  2:927893/113‑1 (MQ=255)
                                                                                                                                                  gccgcgcTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGAGGGCGCCTATGAAGTg  <  1:2013169/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGATCATTCGCCTTAGGCGGTCTCTCACACGGCAGTGCATTCCACTATTTCCCAATGCTACAGTAATGCGAATGCGATCTGGTATGCATATCAACTGATGTGAGTTTGATTAATGAGCA‑CCAACTTAAGTGAAGAAGAAATGCGCGCCGCGCTTTTCGGAGTGGTTAATGTTATCCCTCAAGAAAACTCCCTGCGCTCATCTTGCAGTCCTGCTAAGGCAGTATCAAAAAACTCGGGGACAAAAATTAAAGTTACTTTGCACGTAGGCAACGAGTTTGAGGGCGCCTATGAAGTG  >  NC_002947/2216034‑2216328

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