Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,732,025 A→T 5.8% L196* (TTG→TAG) ‡ hutI ← imidazolonepropionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,732,0250AT 5.8% 314.5 / 9.1 105L196* (TTG→TAG) ‡hutIimidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base A (47/50);  new base T (3/3);  total (51/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 106 total aligned reads displayed.

CGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTTG  >  NC_002947/5731879‑5732174
                                                                                                                                                   |                                                                                                                                                      
cGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGAGCCGGCGAGAACGTGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCGTCGCCGGCCAAgg                                                                                                                                                     >  2:1030507/1‑150 (MQ=255)
    cAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCgg                                                                                                                                                 >  2:1398102/1‑150 (MQ=255)
     aaCGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTTGCCGGCCTTGGCCGGc                                                                                                                                                >  1:1703833/1‑150 (MQ=255)
     aaCGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGc                                                                                                                                                >  2:1370290/1‑150 (MQ=255)
     aaCGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGc                                                                                                                                                <  1:1538867/150‑1 (MQ=255)
       cGATTACAGTTGTTCGGCAGGCAGATTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGGACAGCCGCTCGACCTGGGCCGGCGATAACGCGAGGTGTCCGCAGAAGCCAGCCACCGCCCCCAGCAGGCCTTCGCAGGCCAAGGCCGGCaa                                                                                                                                              <  2:340697/150‑1 (MQ=255)
        gaggtgaGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAAc                                                                                                                                             <  2:400135/144‑1 (MQ=255)
        gATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAAc                                                                                                                                             >  2:164866/1‑150 (MQ=255)
         aTGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTTGCCGGCCTTGGCCGGcgaag                                                                                                                                            <  2:1703833/150‑5 (MQ=255)
         aTGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGATCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCTTGGCCGGcgaag                                                                                                                                            <  2:1022636/150‑5 (MQ=255)
             cAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGGCCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAtctc                                                                                                                                        <  1:505800/150‑1 (MQ=255)
              aGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCg                                                                                                                                       >  2:121772/1‑150 (MQ=255)
               gttgttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGt                                                                                                                                      <  1:189352/150‑1 (MQ=255)
                ttgttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTc                                                                                                                                     >  1:558476/1‑150 (MQ=255)
                 tgttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTca                                                                                                                                    >  1:767717/1‑150 (MQ=255)
                   ttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTcaca                                                                                                                                  <  2:827037/150‑1 (MQ=255)
                   ttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGGcaca                                                                                                                                  >  1:585796/1‑150 (MQ=255)
                    tCGGCATGCAGCTTGCCCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAg                                                                                                                                 <  1:189971/150‑1 (MQ=255)
                      ggCATGCAGCTTGACGGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACATAGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAt                                                                                                                               <  2:357176/150‑1 (MQ=255)
                      ggCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAt                                                                                                                               <  1:392234/150‑1 (MQ=255)
                        cATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAtgt                                                                                                                             <  2:657935/150‑1 (MQ=255)
                         aTGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCg                                                                                                                                                  <  2:1049831/128‑1 (MQ=255)
                         aTGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCg                                                                                                                                                  >  1:1049831/1‑128 (MQ=255)
                          tGCAGCTTGACCGGCAGGCCCAGTGCGCGCGGCTTGAGGGACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGTCAAGGCCGGCAACATCTCGTCACAGTTGTgc                                                                                                                           <  2:391701/150‑1 (MQ=255)
                            cAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTgcgc                                                                                                                         <  1:197076/150‑1 (MQ=255)
                              gCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGACCAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGa                                                                                                                       <  2:1399998/150‑1 (MQ=255)
                              gCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGa                                                                                                                       <  1:598397/150‑1 (MQ=255)
                               cTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCTGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGAt                                                                                                                      >  2:1314790/1‑150 (MQ=255)
                                  gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTa                                                                                                                   <  2:1331468/150‑1 (MQ=255)
                                  gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTa                                                                                                                   <  2:558476/150‑1 (MQ=255)
                                   aCCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAg                                                                                                                  >  2:1634707/1‑150 (MQ=255)
                                         aGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGtcgtcg                                                                                                            <  2:939733/150‑1 (MQ=255)
                                           gCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCTGTCCTTCGCCGGCCAAGGCCg                                                                                                                                                  >  1:27021/1‑110 (MQ=255)
                                           gCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCg                                                                                                                                                  <  2:27021/110‑1 (MQ=255)
                                            cccAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGcc                                                                                                        >  1:905517/1‑150 (MQ=255)
                                             ccAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGccc                                                                                                       >  1:214674/1‑150 (MQ=255)
                                             ccAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGccc                                                                                                       >  2:299404/1‑150 (MQ=255)
                                              cAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCGACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCg                                                                                                      >  1:632916/1‑150 (MQ=255)
                                                 ttCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGCTGTGCGCGACGTAGGGGGCTGCCCGAcc                                                                                                   >  2:901318/1‑150 (MQ=255)
                                                  tCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGAccc                                                                                                  <  1:1244286/150‑1 (MQ=255)
                                                   cgcgcgCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCcg                                                                                                 <  1:299404/150‑1 (MQ=255)
                                                     cgcgcCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCcgcg                                                                                               <  2:1215488/150‑1 (MQ=255)
                                                      gcgcCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATACACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGt                                                                                              >  2:65046/1‑150 (MQ=255)
                                                          cTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTc                                                                                          >  2:955463/1‑150 (MQ=255)
                                                            tggtgaACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCg                                                                                        >  1:1221825/4‑150 (MQ=255)
                                                             gatgaACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCgg                                                                                       >  2:1761762/1‑150 (MQ=255)
                                                                gaACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGTTCCGCGTACTCCggggg                                                                                    <  1:1499218/150‑1 (MQ=255)
                                                                  aCAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCa                                                                                  >  1:439266/1‑150 (MQ=255)
                                                                    aGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCaaa                                                                                <  1:1350468/150‑1 (MQ=255)
                                                                    aGCCGCTCGACCTGGGCCGGCGAGAAAACTAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCaaa                                                                                <  1:1435894/150‑1 (MQ=255)
                                                                     gCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAg                                                                               >  1:436146/1‑150 (MQ=255)
                                                                     gCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGGTCGCAGAAGGCATCCACCGCATCCACCAGGCCGTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAg                                                                               <  2:130823/150‑1 (MQ=255)
                                                                        gCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAtctc                                                                                                                                        <  2:45894/91‑1 (MQ=255)
                                                                        gCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAtctc                                                                                                                                        >  1:45894/1‑91 (MQ=255)
                                                                        gCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGccc                                                                                                       <  1:748712/123‑1 (MQ=255)
                                                                        gCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGccc                                                                                                       >  2:748712/1‑123 (MQ=255)
                                                                           cGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTAATCCGGGGGCAAAGCGTgcg                                                                         <  1:65046/150‑1 (MQ=255)
                                                                                  ggCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCaaa                                                                                >  1:401750/1‑136 (MQ=255)
                                                                                      ggCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTc                                                                                                                                     >  2:1624502/1‑80 (MQ=255)
                                                                                         gagaACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTg                                                           <  1:121772/150‑1 (MQ=255)
                                                                                         gagaACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGAGGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTg                                                           >  2:800460/1‑150 (MQ=255)
                                                                                               gcgAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCggg                                                                                      <  1:1107879/117‑1 (MQ=255)
                                                                                               gcgAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCggg                                                                                      >  2:1107879/1‑117 (MQ=255)
                                                                                               gcgAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgcacgc                                                     >  2:804668/1‑150 (MQ=255)
                                                                                               gcgAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgcacgc                                                     >  2:192720/1‑150 (MQ=255)
                                                                                                 gAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgcacgcac                                                   >  1:1285830/1‑150 (MQ=255)
                                                                                                    gtgtTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGt                                                <  2:767717/150‑1 (MQ=255)
                                                                                                       ttCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCa                                                                   <  1:1693551/128‑1 (MQ=255)
                                                                                                       ttCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCa                                                                   >  2:1693551/1‑128 (MQ=255)
                                                                                                       ttCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCAGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAg                                             >  1:1404496/1‑150 (MQ=255)
                                                                                                          gCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCgg                                          <  1:804668/150‑1 (MQ=255)
                                                                                                           cAGAAGGCATCCACCGCATCCACCAGGCCTTCGTCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGc                                         <  2:984160/150‑1 (MQ=255)
                                                                                                           cAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGc                                         <  1:1761762/150‑1 (MQ=255)
                                                                                                             gAAGGCATCCACCGCATCCACCAGGCCTTCGCACGAACCAGCAAGAAACCCCGTGGCTCAGATGGGCGCGGTGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAg                                       >  1:976465/1‑150 (MQ=255)
                                                                                                             gAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAg                                       <  2:1221825/150‑1 (MQ=255)
                                                                                                                  catccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgt                                  >  1:500348/1‑150 (MQ=255)
                                                                                                                   atccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCCCCGTCGGCGGCAGCtcgtc                                 >  2:110341/1‑150 (MQ=255)
                                                                                                                    tccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgtcg                                >  2:1137013/1‑150 (MQ=255)
                                                                                                                    tccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgtcg                                >  1:489070/1‑150 (MQ=255)
                                                                                                                                accaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGCTGTAGTCGGCGGCCCGTCCCGCGTACTCCGGGGTCAAAGCGTGCGCCGCCATGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGcggcgg                    <  2:214674/150‑1 (MQ=255)
                                                                                                                                accaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGcggcgg                    <  1:1314790/150‑1 (MQ=255)
                                                                                                                                accaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGcggcgg                    >  1:500991/1‑150 (MQ=255)
                                                                                                                                accaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGcggcgg                    <  2:1913860/150‑1 (MQ=255)
                                                                                                                                 ccaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCggg                   <  1:1030507/150‑1 (MQ=255)
                                                                                                                                 ccaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCggg                   >  1:1853739/1‑150 (MQ=255)
                                                                                                                                 ccaGGCCTTCGCCGGCCAAGGCCGATAGGCGCGCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCggg                   >  1:528180/1‑150 (MQ=255)
                                                                                                                                    ggCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGccc                >  2:1297580/1‑149 (MQ=255)
                                                                                                                                     gCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCAc               <  2:489070/150‑1 (MQ=255)
                                                                                                                                       cTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCa                                        >  1:74865/1‑123 (MQ=255)
                                                                                                                                       cTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCa                                        <  2:74865/123‑1 (MQ=255)
                                                                                                                                       cTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTc             >  1:945680/1‑150 (MQ=255)
                                                                                                                                       cTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTc             >  1:343659/1‑150 (MQ=255)
                                                                                                                                       cTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTc             <  2:343659/150‑1 (MQ=255)
                                                                                                                                        ttCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCg            <  1:1137013/150‑1 (MQ=255)
                                                                                                                                         ttgCCGGCCTTGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTgc                                                                          <  1:1597319/85‑1 (MQ=255)
                                                                                                                                         ttgCCGGCCTTGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTgc                                                                          >  2:1597319/3‑87 (MQ=255)
                                                                                                                                         ttgCCGGCCTTGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTGGACGGCGAGGCGGGGGGCCACTCgc           >  1:1662767/3‑150 (MQ=255)
                                                                                                                                            ccGGCCCAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCgcagc        <  1:1022066/150‑1 (MQ=255)
                                                                                                                                                gCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATc     >  1:1616064/1‑150 (MQ=255)
                                                                                                                                                   aaGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGCCGTAGGTAGCGATCGTTTTCTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTTg  <  2:528180/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTTG  >  NC_002947/5731879‑5732174

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