Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT40.4% 132.9 / 53.4 104S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (16/46);  new base T (42/0);  total (58/46)
Fisher's exact test for biased strand distribution p-value = 3.50e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604195‑4604492
                                                                                                                                                     |                                                                                                                                                     
aCTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttct                                                                                                                                                      <  4:206345/151‑1 (MQ=255)
   cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttt                                                                                                                                                    <  3:1043634/150‑1 (MQ=255)
    aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATTCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                                                                                  <  3:1233725/151‑1 (MQ=255)
      aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                <  2:103563/151‑1 (MQ=255)
      aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                                <  1:1273156/151‑1 (MQ=255)
       gAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                                                                               <  1:846294/151‑1 (MQ=255)
        aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                              <  4:422678/151‑1 (MQ=255)
        aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                              <  2:559811/151‑1 (MQ=255)
         aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCAtttttttttttt                                                                                                                                                  >  1:134273/1‑146 (MQ=255)
          aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttttttttct                                                                                                                                                >  1:871622/1‑147 (MQ=255)
              gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATga                                                                                                                                         <  1:262893/150‑1 (MQ=255)
                  cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAg                                                                                                                                     >  1:661402/1‑150 (MQ=255)
                   aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                                                                                   <  1:911275/151‑1 (MQ=255)
                     gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTTATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAg                                                                                                                                 <  4:569978/151‑1 (MQ=255)
                     gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGa                                                                                                                                  <  3:830661/150‑1 (MQ=255)
                       gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                                               <  4:294951/151‑1 (MQ=255)
                        tagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgag                                                                                                                               >  3:884647/1‑148 (MQ=255)
                          gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                                                                                            <  1:607148/151‑1 (MQ=255)
                             tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGTGAGCGCag                                                                                                                            >  4:479810/1‑148 (MQ=255)
                                 tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACCTTTTTTTTTTTTCTGCCAGGGGAAGTGAGCGGAGAGAggg                                                                                                                      >  4:32629/1‑150 (MQ=255)
                                  aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAggg                                                                                                                      >  2:565186/1‑149 (MQ=255)
                                  aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                                                                                     <  1:732714/150‑1 (MQ=255)
                                     gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTAt                                                                                                                 >  1:532887/1‑151 (MQ=255)
                                        agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGGGCGGAGGGAGGGAGTAtttt                                                                                                              >  1:1312609/1‑151 (MQ=255)
                                        agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGGAGAGggggggtatttt                                                                                                              >  1:934926/1‑143 (MQ=255)
                                         gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                              >  2:539180/1‑150 (MQ=255)
                                          tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTAGAAGGGAGCGCAGAGAGGGGTTATCtt                                                                                                              >  1:499808/1‑149 (MQ=255)
                                           ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                                                                            <  2:18542/150‑1 (MQ=255)
                                               tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCctg                                                                                                         >  3:885972/1‑149 (MQ=255)
                                                cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                      <  2:550752/151‑1 (MQ=255)
                                                 aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactg                                                                                                     <  3:528309/151‑1 (MQ=255)
                                                      gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTg                                                                                                <  4:389362/151‑1 (MQ=255)
                                                       tttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCACGAGAAGTGAGCGCGGAGGGGGGTTATCTTGCCTTACTGAAc                                                                                                  >  1:876849/1‑148 (MQ=255)
                                                        ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGGAGGGGGGGGAGGGAGGGGTTATCTTGGCTGGCTGAACTggg                                                                                              >  1:610064/1‑150 (MQ=255)
                                                        ttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAAc                                                                                                  >  4:260859/1‑147 (MQ=255)
                                                         tGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGTGAGCGCAGAGAGGGGGTTTCTTTTCTGTCTGaa                                                                                                   >  1:649272/1‑145 (MQ=255)
                                                         tGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGAGGGATTATCTTGCCTGACTGAACTgg                                                                                               >  3:709866/1‑149 (MQ=255)
                                                          gCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCg                                                                                             <  3:1117640/150‑1 (MQ=255)
                                                                 aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAAATGGCGATTAAGtt                                                                                    >  3:509317/1‑151 (MQ=255)
                                                                 aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGtt                                                                                    <  2:373618/151‑1 (MQ=255)
                                                                 aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCGGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACGGAACTGGCGATTAAGtt                                                                                    <  2:398324/151‑1 (MQ=255)
                                                                  cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCGATTAAGttt                                                                                   >  1:1029307/1‑151 (MQ=255)
                                                                  cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt                                                                                   >  2:330883/1‑151 (MQ=255)
                                                                    cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCACTTTTTTTTTTTCTGTCAAGAGCAGTGAAGGCCGGGAGGGGGTATTTTTCGCGACTGGACTGGCGATTAAGctttc                                                                                 >  3:1412544/1‑149 (MQ=255)
                                                                      tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                                 <  4:270562/149‑1 (MQ=255)
                                                                         aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGTATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGGAGGGAGGGGTTATTTTGCCTGTATGAACTGGGGGTTAAGTTTTcgc                                                                               >  1:1079189/1‑148 (MQ=255)
                                                                              ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                                       <  3:1304485/151‑1 (MQ=255)
                                                                                aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCg                                                                     <  3:72618/151‑1 (MQ=255)
                                                                                  aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTAGTAGGGAGCGCAGAGAGGGGGTATCTTGCCTGACTGAAATGGGGATTAAGTTTTCGCGTAccc                                                                         >  1:21876/1‑145 (MQ=255)
                                                                                  aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCCGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCggt                                                                   >  3:474897/1‑149 (MQ=255)
                                                                                        aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATTAGAAGTGAGCGCCGAGGGGGGTTATTTTGCCTGGCTGGACTGGGGATTAAAGTTTCGCGGAccccggggtttgtt                                                              >  1:186904/1‑139 (MQ=255)
                                                                                        aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGtt                                                              <  1:744915/150‑1 (MQ=255)
                                                                                         ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                            <  3:248447/151‑1 (MQ=255)
                                                                                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGGAGGGAGCGCAGAGAGGGGTTATTTTTCCTTGCTGGACCGGCGGTTAAGTTTCCCCGTACCCGGCGGTTAGTtttt                                                           >  1:1268109/1‑149 (MQ=255)
                                                                                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGGGAGCGCAGAGGGGGGTTATTTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCGTTTTGTTAtt                                                           >  2:37638/1‑149 (MQ=255)
                                                                                             tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATATTGCCTGGCTGAACTTGCGATTAAGTTTCCGCGTACCCGGCGTTTGGTtttttgg                                                        >  1:177349/1‑148 (MQ=255)
                                                                                               aGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                                      <  4:30306/151‑1 (MQ=255)
                                                                                                     cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgatg                                                 <  1:474821/150‑1 (MQ=255)
                                                                                                        ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCCGAGAGGGGTTATTTTGCCTGACTGCACTGGCGGTTAAGTTTTCGCGTACCCCGGGTTTTGTTATTGCGGGTGGTGgtat                                             >  1:331211/1‑151 (MQ=255)
                                                                                                        ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGata                                              <  1:41525/150‑1 (MQ=255)
                                                                                                           aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTCTTCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGGTGATGATAtt                                            >  2:312148/1‑149 (MQ=255)
                                                                                                           aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTGACTGAAATGGCGGTTAAGTTTTCGCGTACCCGGCGGTTTGTTATTGCGGGTGATGATAtt                                            >  3:162588/1‑149 (MQ=255)
                                                                                                            ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGaa                                         <  1:121518/151‑1 (MQ=255)
                                                                                                             gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGGAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTTAGTTTTCGCGTACCCCGCGGTTTGTTTTTTCGGGTGATGATAtt                                            >  1:873461/1‑147 (MQ=255)
                                                                                                             gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAg                                        <  1:859961/151‑1 (MQ=255)
                                                                                                               tGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGGTTAGTTATTGCGGGTGATGATAtt                                            >  3:239414/1‑145 (MQ=255)
                                                                                                               tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGGGAGCGCGGAGAGGGATTATCTTGCCTGACTGGACTGGCGGTTAAATTTTCGCGTACCCCGCGTTTTGTTATTTCGGGTGATTATATTGaa                                         >  3:667578/1‑148 (MQ=255)
                                                                                                                gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                                      <  3:84696/150‑1 (MQ=255)
                                                                                                                 ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGGAGGGAGCGCAGGGAGGGGTTATCTTTCCTGACTGAACTGGCGCTTAAGTTTTCGCGTACCCGGCGGTTTGTTATTGCGGATGATGATATTTAAACTc                                     >  4:481350/1‑150 (MQ=255)
                                                                                                                 ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGGGCGCAGAGGGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCCGCGTTTTGTTATTGCGGGTGGTGGTATTTAAGCTc                                     >  1:1029319/1‑150 (MQ=255)
                                                                                                                  cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGGTTAGTTATTGCGGGTGATGATATTGAAGCTCgg                                   >  2:53075/1‑151 (MQ=255)
                                                                                                                  cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCgg                                   <  3:1414950/151‑1 (MQ=255)
                                                                                                                   aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                  <  3:723412/151‑1 (MQ=255)
                                                                                                                     gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTg                                <  2:202677/151‑1 (MQ=255)
                                                                                                                      tCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGGGCGCAGGGAGGGATTATTTTTCCTGACTGGACTGGCGGTTAAATTTTCGCGGACCCGGCGTTTTGTTATTTCGGGTGATGATAATGAAACTCgg                                   >  3:200289/1‑147 (MQ=255)
                                                                                                                       cgcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGTGAGCGCAGAGAGGGGTTTTCTTTCCTGACTGAACTGGCGGTTAAAGTTTCCCGTACCCCGCGTTTTGTTATTGCGGGTGATGATAtt                                            >  1:716683/1‑137 (MQ=255)
                                                                                                                        gcgGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCATGAGGAGTGGGGGCCGAGAGGGGTTATTTTGCCTTACTGAACTTGCGGTTAAGGTTCCGCGGACCCCGGGCTTTGTTATTTCGGATGATGGTATTTTAGATCGGCTgg                               >  4:480957/1‑148 (MQ=255)
                                                                                                                         cgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATTAGAAGTGAGCGCCGGGAGGGGTTATCTTGCCTGTCTGAACTGGCGATTAAATTTTCGCGTAACCGGCGTTTTGTTATTTCGGATGATGATATTGAAGCTCGGCTGAt                              >  4:170579/1‑149 (MQ=255)
                                                                                                                          ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGGAGGGAGCGCAGAGAGGGGTTTTCTTGCCTGACTTAAATGGGGATTAAGTTTTCGCGGACCCGGCGGTTAGTTTTTTCGGGTGATGATATTGAAGCTCGGCTGGtttg                           >  3:1073431/1‑149 (MQ=255)
                                                                                                                           gATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCCCGTACCCGGCGGTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGGtttg                           >  2:135289/1‑148 (MQ=255)
                                                                                                                                cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGt                          <  3:1079318/146‑1 (MQ=255)
                                                                                                                                cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCa                       <  4:315784/149‑1 (MQ=255)
                                                                                                                                cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAt                      <  3:671467/150‑1 (MQ=255)
                                                                                                                                  tttGCGGCGGCCATTTTTTCTTTTTTTGCCATGAGAAGGGAGCGCCGGGAGGGGTTATCTTTCCTGACTGAAATGGCGAATAAATTTTCGCGTACCCCGCGTTTTGTTTTTTCGGGTGGTGATATTGAAGCTCGGCTGATTGGTTCAtt                     >  2:297650/1‑149 (MQ=255)
                                                                                                                                   ttGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGt                   <  1:5902/150‑1 (MQ=255)
                                                                                                                                     gcggcgGTCATTTTTTCTTTTTCTTCCATGAGAAGGGAGCGCAGGGAGGGGTTATCTTGCCTGACTGAAATGGCGATTAAGTTTTCGCGTACCCGGCGTTTAGTTATTTCGGATGATGATATTGAAGCTCGGCTGGTTGGTTCATTGTCg                 >  1:678930/1‑150 (MQ=255)
                                                                                                                                     gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGc                <  1:243374/151‑1 (MQ=255)
                                                                                                                                      cggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCTcc               >  1:929842/1‑151 (MQ=255)
                                                                                                                                        gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt             >  2:15652/1‑151 (MQ=255)
                                                                                                                                        gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAATGGTTTTCTTGCCTGACTGAACGGGCGATTAATTTTCTGCGTACCCCGCGTTTAGTTATTTCGTATTATTATATTGAAGTTCGGCTGATTGGTTCATTGTCGCCAt             <  1:141936/151‑1 (MQ=255)
                                                                                                                                           gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAc          >  3:1278/1‑151 (MQ=255)
                                                                                                                                            tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc         >  3:848663/1‑151 (MQ=255)
                                                                                                                                            tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc         >  1:1050790/1‑151 (MQ=255)
                                                                                                                                             cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc         <  4:227808/150‑1 (MQ=255)
                                                                                                                                             cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCg        >  2:451336/1‑151 (MQ=255)
                                                                                                                                              aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGc       <  3:1058778/151‑1 (MQ=255)
                                                                                                                                                tttttctttttcTGCCATGAGAAGTGAGCGCGGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGCTTGGTTCATTGTCGCCATCACCGCTg     >  3:306842/1‑151 (MQ=255)
                                                                                                                                                tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt      <  3:532962/150‑1 (MQ=255)
                                                                                                                                                tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt      <  3:510305/150‑1 (MQ=255)
                                                                                                                                                tttctctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg     >  1:752568/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604195‑4604492

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