Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT50.0% 42.5 / 42.6 64S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (4/28);  new base T (32/0);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 6.43e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604195‑4604490
                                                                                                                                                     |                                                                                                                                                   
aCTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTtttttct                                                                                                                                                    >  1:708628/1‑151 (MQ=255)
  tCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttt                                                                                                                                                  >  1:101748/1‑151 (MQ=255)
   cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttt                                                                                                                                                 <  2:504069/151‑1 (MQ=255)
   cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttt                                                                                                                                                  <  1:713679/150‑1 (MQ=255)
       gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACTCttttttttttttt                                                                                                                                              >  1:717486/1‑148 (MQ=255)
        aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttt                                                                                                                                              >  2:697769/1‑147 (MQ=255)
           aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGccct                                                                                                                                         >  1:1224155/1‑149 (MQ=255)
           aaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCAt                                                                                                                                         <  2:83425/151‑1 (MQ=255)
              gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag                                                                                                                                      <  1:439645/151‑1 (MQ=255)
               cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCtttttttttttt                                                                                                                                                >  1:556118/1‑140 (MQ=255)
                gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGaa                                                                                                                                    <  2:965068/151‑1 (MQ=255)
                  cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGaa                                                                                                                                    >  1:242010/1‑149 (MQ=255)
                      agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGa                                                                                                                                >  1:1223833/1‑149 (MQ=255)
                      agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                              <  2:896872/151‑1 (MQ=255)
                       gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAgcg                                                                                                                             >  2:1248343/1‑151 (MQ=255)
                        tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGAGAgcg                                                                                                                             >  1:665972/1‑150 (MQ=255)
                               ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAgg                                                                                                                     <  2:116910/151‑1 (MQ=255)
                                 tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCCGAGAggg                                                                                                                    >  2:770627/1‑150 (MQ=255)
                                  aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                                                                                   <  2:907869/150‑1 (MQ=255)
                                     gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTa                                                                                                                <  1:814170/150‑1 (MQ=255)
                                                  gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactga                                                                                                  <  2:381090/151‑1 (MQ=255)
                                                   aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactga                                                                                                  <  1:1156560/150‑1 (MQ=255)
                                                    aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGGGGGATTTTTTTTCCTGACTGAAc                                                                                                >  2:576027/1‑151 (MQ=255)
                                                      gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑Tg                                                                                             <  2:550019/151‑1 (MQ=255)
                                                       tttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGAGGGGTTATCTTTCCTGGCTGAAA‑Tgg                                                                                            >  1:234705/1‑151 (MQ=255)
                                                        ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAGGGATTATCTTGCCTGACTGaa                                                                                                 >  2:355114/1‑146 (MQ=255)
                                                          gCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCg                                                                                          <  2:746909/150‑1 (MQ=255)
                                                          gCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGa                                                                                         <  2:353074/151‑1 (MQ=255)
                                                            gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGAtt                                                                                       <  1:669622/151‑1 (MQ=255)
                                                              gTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGGGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGAtt                                                                                       >  2:1202696/1‑149 (MQ=255)
                                                                aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCAGAGGGGGGTTATTTTTCCTTACTGAAC‑TGGGGATTaa                                                                                     >  1:1274609/1‑149 (MQ=255)
                                                                aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCCGAGAGGGGTTATTTTTCCTTACTGAAC‑TGGCGGtt                                                                                       >  1:259728/1‑147 (MQ=255)
                                                                aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAg                                                                                    <  1:1186150/150‑1 (MQ=255)
                                                                   ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTTAGAAGTGAGCGCAGAGAGGGGTTATTTTTCCTTAATGAAC‑TTGCGAGTAAGgtt                                                                                 >  1:419593/1‑150 (MQ=255)
                                                                    cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGGGGGCGCAGAGAGGGGTTATTTTTCCTGACTGAAC‑TGGCGGTTTAGttt                                                                                 >  1:1170647/1‑149 (MQ=255)
                                                                     ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGGGGGGGCCGAGGGGGGTTATCTTGCCTGGCTGAAC‑TGGGGGTTAAGatttc                                                                               >  1:1104968/1‑148 (MQ=255)
                                                                     ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCCTGTGAAGGGAGCGCCGAGAGGGGTTTTCTTGCCTGACTGAAC‑TTGCGGTTTAGttttc                                                                               >  1:853388/1‑148 (MQ=255)
                                                                      tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATTAGGAGGGAGCGCAGAGGGGGGTTATCTTTCCTGACTGCAC‑TGGGGGTTAAGTTTTcg                                                                              >  1:360527/1‑150 (MQ=255)
                                                                      tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTcc                                                                               <  1:212867/149‑1 (MQ=255)
                                                                         aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCATGGGGAGGGGGCGCCGGGAGGGGTTTTCTTTCCTGACTGGAA‑TGGCGATTAAGttt                                                                                 >  2:443302/1‑144 (MQ=255)
                                                                          aGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGGGGGGTTATCTTGCCTGACTGAAC‑CGGCGATTAAGTTTCcgcg                                                                            >  1:694375/1‑148 (MQ=255)
                                                                                aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAATGAGGGCAGAGGGGGGGTATTTTTGCTGACCGAAA‑TGGCGAATAAGATTACGCGGACCCCGc                                                                    >  2:614756/1‑150 (MQ=255)
                                                                                  aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGGAGTGAGCGCCGAGGGGGGTTATTTTTCCTGACTGAAC‑TGGGGGTTAAGTTTTCGCGTACCCgg                                                                     >  1:301165/1‑147 (MQ=255)
                                                                                  aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGCGAAGGGAGCGCAGAGAGGGGTTATTTTGGCTGACTGAAA‑TGGCGATTAAGTTTCCGCGGACCCCGCGtt                                                                 >  1:1223918/1‑151 (MQ=255)
                                                                                   ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGGAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAAC‑TGGGGGTTAAAGTTCCGCGTACCCGGCGgtt                                                                >  2:555932/1‑151 (MQ=255)
                                                                                   ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGGGAGGGGTTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTTCGCGTACCCGGCGgtt                                                                >  1:648079/1‑151 (MQ=255)
                                                                                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGGGCGCCGAGAGGGGTTATTTTTTCTGACTGAAA‑TGGCGGTTAAGATTTCGCGGAACCGGCGTTTTGt                                                             >  1:790984/1‑151 (MQ=255)
                                                                                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTa                                                               <  1:1109165/149‑1 (MQ=255)
                                                                                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCGGTGAAGAGAAGGCTGAGGGGGGTTATTTTGGCTGACTGAAC‑TGGGTGGTGAGTTTTCGCGTACCCGGCGGTTAGTTAtt                                                         >  1:370822/1‑151 (MQ=255)
                                                                                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTGCCCTGAGAAGGGAGGGCGGAGAGGGATTATTTTGCCTGACTGAAC‑TGGCGGTTAAGTTTCCGCGGACCCCGCGttttggtttt                                                         >  2:325472/1‑143 (MQ=255)
                                                                                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTACCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                                        <  1:1216792/151‑1 (MQ=255)
                                                                                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGATTTTTTTTCCTGACTGAAC‑TGGGGGTTAAAATTCCGCGTACCCGGCGTTTTGTtttt                                                         >  2:1066539/1‑149 (MQ=255)
                                                                                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                         <  1:719491/149‑1 (MQ=255)
                                                                                             tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTCTTTTTTTCTCCTGAGAAGGGGGCGCCACGGGGGAGTATCTTGCCTGACCGAACGTGGGCATTAGGTTTCC‑CGTACCCGCCGGTTTgttgtttgg                                                      >  2:53839/1‑148 (MQ=255)
                                                                                              cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGAACCCGGCGTTTAGTTATTGCg                                                      <  1:995214/150‑1 (MQ=255)
                                                                                                 gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatg                                                  <  2:1267778/151‑1 (MQ=255)
                                                                                                     cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATTTTTCCTGACTGAAC‑TGGCGGTTAAGTTTTCGCGTACCCCGCGGTTTGTTTTTGCGgatgat                                                >  1:919966/1‑149 (MQ=255)
                                                                                                     cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGGTTTGTTATTGCGggtg                                                  >  2:612697/1‑147 (MQ=255)
                                                                                                      aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgata                                               <  1:1211808/149‑2 (MQ=255)
                                                                                                          cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                                          <  2:87943/150‑1 (MQ=255)
                                                                                                             gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                      <  2:235486/151‑1 (MQ=255)
                                                                                                                       cgcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGa                             <  1:406141/150‑1 (MQ=255)
                                                                                                                         cgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGGTATTTTGTCTGAATGAAC‑TGGCGATTTAGTTTTCGCGTACCCGGCGTTTTGTTATTTCGGGTGATGATATTTAAGCTCGGCTGGtt                           >  1:600508/1‑150 (MQ=255)
                                                                                                                              ctcGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGAGGGGGGGGAGGGGGGCTATTTTTCCTTACTGAAC‑TGGGGGTTAAAGTTTCGCGGAGCCGGGGGTTTGTTATTGCGGTGGATGGTATTGAAACCCGGCTGATTTGtt                       >  1:277020/1‑149 (MQ=255)
                                                                                                                                cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTc                      <  2:1228705/148‑1 (MQ=255)
                                                                                                                                cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa                     <  2:908409/149‑1 (MQ=255)
                                                                                                                                cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa                     <  1:563065/149‑1 (MQ=255)
                                                                                                                                  tttGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGTGAGCGCCGAGAGGGGTTATTTTGCCTGACTGAAC‑TGGCGATTAAGTTTTCGCGTACCCGGCGTTTTGTTATTTCGGGTGATGATATTTAAGCTCGGCTGGTTTGTTCATTGt                 >  2:1043976/1‑151 (MQ=255)
                                                                                                                                      cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGcc             <  2:393482/151‑1 (MQ=255)
                                                                                                                                        gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt           >  1:800718/1‑151 (MQ=255)
                                                                                                                                            tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc       >  1:174366/1‑151 (MQ=255)
                                                                                                                                                tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt    <  2:708680/150‑1 (MQ=255)
                                                                                                                                                tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt    <  2:1087247/150‑1 (MQ=255)
                                                                                                                                                 ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg  >  1:541107/1‑151 (MQ=255)
                                                                                                                                                 ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg   <  2:895205/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAC‑TGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGG  >  NC_000913/4604195‑4604490

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