Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,325,559 T→A 7.7% noncoding (133/149 nt) PP_mr52 → PrrF ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,325,5590TA 7.7% 273.1 / 9.7 88noncoding (133/149 nt)PP_mr52PrrF ncRNA
Reads supporting (aligned to +/- strand):  ref base T (42/39);  new base A (5/2);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACCTT  >  NC_002947/5325416‑5325708
                                                                                                                                               |                                                                                                                                                     
cGGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATCGACCCGGCAATTTGCCg                                                                                                                                                 <  2:2360419/150‑1 (MQ=255)
  gggCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCggg                                                                                                                                               >  1:1240653/1‑150 (MQ=255)
      aaaGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCtt                                                                                                                                           >  2:2149198/1‑150 (MQ=255)
       aaGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCttt                                                                                                                                          <  2:928384/150‑1 (MQ=255)
         gggATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCttttt                                                                                                                                        >  1:768628/1‑150 (MQ=255)
         gggATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCttttt                                                                                                                                        <  1:1995534/150‑1 (MQ=255)
                 cgtcTAAATGAGAATTGTTATAAGTATCGCAACTGGTCGCGAGACTGCTGGGATAAGCCAAGAGCCCCTGCTTCGGACTCTCAGATTAGCTCCTCATCAGGCTAATCACGCCTATTGACCCGCCAATTTCCCGGGCCTATTTTTGCCtgt                                                                                                                                <  2:121189/147‑1 (MQ=255)
                        aTGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGAGTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTCGCCTGTGCCtttt                                                                                                                         >  1:1140444/1‑150 (MQ=255)
                         tGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCTGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTc                                                                                                                        >  2:211995/1‑150 (MQ=255)
                         tGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTc                                                                                                                        <  1:212637/150‑1 (MQ=255)
                          gagaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                       <  2:1841140/150‑1 (MQ=255)
                             aaTTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCttttttt                                                                                                                                      <  2:1914570/132‑1 (MQ=255)
                             aaTTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCttttttt                                                                                                                                      >  1:1914570/1‑132 (MQ=255)
                              aTTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGAt                                                                                                                   <  1:211995/150‑1 (MQ=255)
                                      tattatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGcgccg                                                                                                           >  2:1219270/1‑150 (MQ=255)
                                       attatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGcgccgc                                                                                                          <  1:2003364/150‑1 (MQ=255)
                                        ttatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGGCCCGGCAATTTGCCGGGTCTTTTTTTGCATGTGCCTTTTCAGGATCAGCGCCGCa                                                                                                         >  1:1322925/1‑150 (MQ=255)
                                                 cTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGc                                                                                                                                    >  1:1670768/1‑114 (MQ=255)
                                                 cTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGc                                                                                                                                    <  2:1670768/114‑1 (MQ=255)
                                                       gcgAGACTGGTTGGAACAGCTGAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACaataat                                                                                          >  2:2020701/1‑150 (MQ=255)
                                                        cgAGACTGGTTGAATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATc                                                                                         <  2:336963/150‑1 (MQ=255)
                                                         gagaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCg                                                                                                                                                 <  2:659684/93‑1 (MQ=255)
                                                         gagaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCg                                                                                                                                                 >  1:1856156/1‑93 (MQ=255)
                                                         gagaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCg                                                                                                                                                 <  2:1856156/93‑1 (MQ=255)
                                                         gagaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCg                                                                                                                                                 >  1:659684/1‑93 (MQ=255)
                                                         gagaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCACGATCAGCGCCGCAGCTGCGCACAATAATcc                                                                                        >  2:825701/1‑150 (MQ=255)
                                                              tGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCttttttt                                                                                                                                      >  2:1560937/1‑99 (MQ=255)
                                                              tGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCttttttt                                                                                                                                      <  1:1560937/99‑1 (MQ=255)
                                                               ggTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTcc                                                                                  <  1:1219270/150‑1 (MQ=255)
                                                                gctGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCg                                                                                 >  1:567298/3‑150 (MQ=255)
                                                                     aTAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                       >  2:896348/1‑107 (MQ=255)
                                                                     aTAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                       <  1:896348/107‑1 (MQ=255)
                                                                          cTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCt                                                                                    >  1:1996592/1‑137 (MQ=255)
                                                                          cTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCt                                                                                    <  2:1996592/137‑1 (MQ=255)
                                                                             agagCCCCTTGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCgt                                                                    <  1:911846/150‑1 (MQ=255)
                                                                             agagCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCACGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCgt                                                                    <  2:342195/150‑1 (MQ=255)
                                                                                   ccTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACca                                                              <  1:711511/150‑1 (MQ=255)
                                                                                    cTGGTTCGGACTCTCAGATTATCTCCTCATCAAGCTAATCACGGTTATTGACTCGGCAATTTGCCGGGTCGTTTTTTGCCTGTGAGTTTTCGGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCCGCATCGCCGGGGTGTGCcac                                                             >  1:1748227/1‑150 (MQ=37)
                                                                                      ggTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGCATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCAGGCATCGCCGCCGTGTACcacac                                                           <  2:294408/150‑1 (MQ=255)
                                                                                       gTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCtgt                                                                                                                                >  2:2396256/1‑80 (MQ=255)
                                                                                       gTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCtgt                                                                                                                                <  1:2396256/80‑1 (MQ=255)
                                                                                       gTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTgcg                                                                                                   <  2:1604064/109‑1 (MQ=255)
                                                                                       gTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTgcg                                                                                                   >  1:1604064/1‑109 (MQ=255)
                                                                                        ttCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACcacacaa                                                         >  2:1740072/1‑149 (MQ=255)
                                                                                        ttCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACAt                                                         >  1:2246453/1‑150 (MQ=255)
                                                                                          cGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATaa                                                       <  1:1348624/150‑1 (MQ=255)
                                                                                           ggACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCGGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATGAGCGCCGCATCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAAt                                                      >  2:53360/1‑150 (MQ=255)
                                                                                             aCTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTCGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCa                                                    <  2:1140444/150‑1 (MQ=255)
                                                                                                 tcAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCCCCGTGTACCACACATAATCAGCCt                                                >  2:420445/1‑150 (MQ=255)
                                                                                                  cAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTc                                               <  2:768628/150‑1 (MQ=255)
                                                                                                      ttATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCCAGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCa                                           <  2:2373148/150‑1 (MQ=255)
                                                                                                       tATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAg                                          <  1:2149198/150‑1 (MQ=255)
                                                                                                         tctcCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGc                                        <  2:1226653/150‑1 (MQ=255)
                                                                                                           tcctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCa                                                                                                                 >  1:360605/1‑75 (MQ=255)
                                                                                                           tcctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCa                                                                                                                 <  2:360605/75‑1 (MQ=255)
                                                                                                            cctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCg                                     <  2:1972784/150‑1 (MQ=255)
                                                                                                             ctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACca                                                              >  1:363658/1‑124 (MQ=255)
                                                                                                             ctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACca                                                              <  2:363658/124‑1 (MQ=255)
                                                                                                                atcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtc                                 >  1:1383352/1‑150 (MQ=255)
                                                                                                                   aGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgct                              >  1:1027697/1‑150 (MQ=255)
                                                                                                                   aGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgct                              >  1:1310865/1‑150 (MQ=255)
                                                                                                                     gCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctcg                            >  1:626122/1‑150 (MQ=255)
                                                                                                                       tAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTc                                               <  2:694012/129‑1 (MQ=255)
                                                                                                                       tAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTc                                               >  1:694012/1‑129 (MQ=255)
                                                                                                                       tAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcg                                <  1:1598304/144‑1 (MQ=255)
                                                                                                                       tAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcg                                >  2:1598304/1‑144 (MQ=255)
                                                                                                                            aCGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTAc                                                                <  2:1962932/107‑1 (MQ=255)
                                                                                                                            aCGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTAc                                                                >  1:1962932/1‑107 (MQ=255)
                                                                                                                              ggTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTca                   >  1:494868/1‑150 (MQ=255)
                                                                                                                              ggTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTca                   >  2:2176742/1‑150 (MQ=255)
                                                                                                                                    agACCCGGCAAATTGCAGGGTCTTTTTTTGACTGTGGCTCTTCGCGGTCAGCGCCGCGGCTTCGCACAGCAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTcagccagc             >  2:1480775/2‑150 (MQ=255)
                                                                                                                                     gACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCt                                       >  1:1831538/1‑123 (MQ=255)
                                                                                                                                     gACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCt                                       <  2:1831538/123‑1 (MQ=255)
                                                                                                                                      aCCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCcagcag           >  2:710367/1‑150 (MQ=255)
                                                                                                                                      aCCCGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGTGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAgg                                         <  1:1648059/120‑1 (MQ=255)
                                                                                                                                      aCCCGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGTGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAgg                                         >  2:1648059/1‑120 (MQ=255)
                                                                                                                                      aCCCGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCcagcag           >  1:1495804/1‑150 (MQ=255)
                                                                                                                                        ccGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACca                                                              >  2:1428259/1‑97 (MQ=255)
                                                                                                                                        ccGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACca                                                              <  1:1428259/97‑1 (MQ=255)
                                                                                                                                        ccGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAgcagca         >  2:1647734/1‑150 (MQ=255)
                                                                                                                                         cGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCt                                                                                    >  1:2075814/1‑74 (MQ=255)
                                                                                                                                         cGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCt                                                                                    <  2:2075814/74‑1 (MQ=255)
                                                                                                                                         cGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCaa        <  1:2176742/150‑1 (MQ=255)
                                                                                                                                         cGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGACTCCCCAAGAGCGGTTCTCTCCCCCACCGCATCCAGCAGCaa        >  1:1736851/1‑150 (MQ=255)
                                                                                                                                          ggCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctc                             >  2:489656/1‑128 (MQ=255)
                                                                                                                                          ggCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctc                             <  1:489656/128‑1 (MQ=255)
                                                                                                                                             aaTTTGCCGGGTCTTTTTTTTTTTTTTTTTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTCCCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATAcc    >  2:2350447/1‑150 (MQ=255)
                                                                                                                                             aaTTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGCATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGACATCAGCAATAcc    >  1:1062689/1‑150 (MQ=255)
                                                                                                                                               tttGCCGGGTCTTTTTTTGCCTGAGCTTATCCAGGTGGAGGTCCTGAGCTGCGCACCATAACCCTAGAGCCGCTTCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACCtt  <  1:420445/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CGGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACCTT  >  NC_002947/5325416‑5325708

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: