Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3500CT44.4% 64.5 / 34.5 65S64F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (8/27);  new base T (28/0);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 5.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAT  >  NC_000913/4604204‑4604499
                                                                                                                                                   |                                                                                                                                                      
aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGCTTGCGGCGGTCATTTTTTTCTTTTTTTGc                                                                                                                                                     >  1:125548/1‑151 (MQ=255)
    tGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGcc                                                                                                                                                    >  1:296362/1‑147 (MQ=255)
     gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag                                                                                                                                               <  1:1386286/151‑1 (MQ=255)
     gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATga                                                                                                                                                <  2:1109993/150‑1 (MQ=255)
       gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGaa                                                                                                                                             <  1:1016099/151‑1 (MQ=255)
         cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTAGaa                                                                                                                                             >  2:369140/1‑149 (MQ=255)
          aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                            <  1:760390/149‑1 (MQ=255)
          aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTATTCCATgatgagcg                                                                                                                                          >  1:621780/1‑145 (MQ=255)
           tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCATGAGaa                                                                                                                                             >  2:560215/1‑147 (MQ=255)
            gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATTAGGAggggg                                                                                                                                        >  1:929631/1‑147 (MQ=255)
            gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGa                                                                                                                                         <  2:1360372/150‑1 (MQ=255)
                  taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTTCCATGAGGAGTGAGGGCag                                                                                                                                   >  1:1059568/1‑150 (MQ=255)
                     cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTGGAAGTGGGGGCcgag                                                                                                                                 >  2:840293/1‑149 (MQ=255)
                      ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAgg                                                                                                                              >  1:476123/1‑151 (MQ=255)
                         aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCGGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCGTTTTTTTTTTTTTTGCCATGAGAAGGGGGGGCAGAGAggg                                                                                                                             >  2:1303577/1‑149 (MQ=255)
                         aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTACTTTTTTTTTTTTCTTCCATGAGAAGTGAGCCCAGAGAggg                                                                                                                             >  1:7905/1‑149 (MQ=255)
                            gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGggggggttt                                                                                                                        >  2:449405/1‑146 (MQ=255)
                            gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTAt                                                                                                                       <  1:1123726/151‑1 (MQ=255)
                              gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCGTTTTTTTTTTTTTTGGCGTTTGAAGGGAGGGCCGAGggggggt                                                                                                                          >  2:614902/1‑144 (MQ=255)
                               agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGAGAAGGGAGGGCAGAGGGGGATTATCtt                                                                                                                    >  1:282437/1‑151 (MQ=255)
                                        aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTTCTGTTTTTAATAATATACCATGGGAAGGGAGGGGGGAGAGGGAGTATCTTGCcttact                                                                                                            >  1:893352/1‑150 (MQ=255)
                                          aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGaa                                                                                                         <  2:587972/151‑1 (MQ=255)
                                           aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGAGTATCTTGCCTGACTGaa                                                                                                         >  2:1156315/1‑150 (MQ=255)
                                              tttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTg                                                                                                      <  1:1315806/150‑1 (MQ=255)
                                               ttGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGTATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTgg                                                                                                     <  1:1374642/150‑1 (MQ=255)
                                                   gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGGGGGCAGAGGGGGGTTATTTTTCcttac                                                                                                             >  2:512481/1‑138 (MQ=255)
                                                    ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTa                                                                                               <  2:784056/151‑1 (MQ=255)
                                                        aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTTCCTGAGAAGAGAGCGCAGCGAGGGGGTTTCTTGCctgactggattcggggtt                                                                                                >  2:28760/1‑135 (MQ=255)
                                                        aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTTGGAAAGGAGAGCCGAGAGGGATTATTTTTGCCCGCTGAAGTGGGGATTTAggt                                                                                           >  1:942388/1‑149 (MQ=255)
                                                         cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGttt                                                                                          <  1:143084/151‑1 (MQ=255)
                                                            ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGTGGGGGCAGAGAGGGATTTTTTTGCCTGACTGAACTGGGGATTTAGgtttc                                                                                        >  1:1277119/1‑148 (MQ=255)
                                                             tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                                        <  2:1064009/149‑1 (MQ=255)
                                                                   cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGGAGTGTGCGCAGAGAGGGGGTATTTTTTCTTATTCAACTTGGGATTAAGTTTCTTCGTAccc                                                                                >  2:31904/1‑151 (MQ=255)
                                                                   cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGGGATTTAGTTTCCGCGTAccc                                                                                >  1:930700/1‑151 (MQ=255)
                                                                   cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTTCAAGAGGAGTGAGGGGAGGGAGGGGTTTTTTTTCCTTACTCCAATGGCCGTTAAGTTTCCCCGTAcc                                                                                 >  2:948153/1‑150 (MQ=255)
                                                                     ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                                                <  2:1164161/149‑1 (MQ=255)
                                                                     ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                                              <  1:716800/151‑1 (MQ=255)
                                                                      cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGGACCCgg                                                                              >  1:394293/1‑150 (MQ=255)
                                                                        tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGAGAGGGCAGGGGGGGGGTTTTTTTCCTGACTGAACTGGCGGTTAAGTTTTCCCGTAACCGGCg                                                                            >  2:936600/1‑150 (MQ=255)
                                                                           cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGGCATTAGAAGTGGGGGCAGAGGGGGGTTATCTTGCCTGAATGAACTAGCGGTTAAGGTTCCGCGTACCCGGCGgtt                                                                         >  1:763663/1‑150 (MQ=255)
                                                                             ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg                                                                       <  2:466723/150‑1 (MQ=255)
                                                                                ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGTCTGACTGAACTGGCGGTTTAGTTTTCGCGGACCCCGCGTTTTGtttt                                                                   >  2:1354601/1‑149 (MQ=255)
                                                                                 cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                                  <  1:613185/151‑1 (MQ=255)
                                                                                    tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc                                                                <  2:273007/150‑1 (MQ=255)
                                                                                       ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCCGAGAGGGGTTATCTTTCCTGAATGAACTTGGGGTTAAGTTTTCCCGTAACCCGCGCTTTGGTTTTGCgg                                                              >  1:654615/1‑149 (MQ=255)
                                                                                       ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                                             <  1:186067/150‑1 (MQ=255)
                                                                                       ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCAATTTTTTTTTTTTTGGCATGGGAAGTGAGGGCCGGGGGGGGTTTTTTTTGCTGACTGAAATGGCGATTTAGTTTTCGCGTACCCGGCGGTTTGTTTTTGCgg                                                              >  2:5605/1‑149 (MQ=255)
                                                                                          gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                                          <  2:621765/150‑1 (MQ=255)
                                                                                                   ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGaa                                                <  2:597189/151‑1 (MQ=255)
                                                                                                          aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCGGCGTTTTGTTTTTGCGGGTGATGATATTTAAGCTCgg                                          >  1:180377/1‑150 (MQ=255)
                                                                                                          aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCGTTTTGTTTTTGCGGGTGGTGGTATTTAAACTCgg                                          >  1:566392/1‑150 (MQ=255)
                                                                                                          aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                         <  2:1311294/151‑1 (MQ=255)
                                                                                                            gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                                        <  2:1033370/150‑1 (MQ=255)
                                                                                                              cgcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGa                                      <  2:48032/150‑1 (MQ=255)
                                                                                                               gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                                    <  1:1247989/151‑1 (MQ=255)
                                                                                                                cgGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGTGAAGTGGGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGTTTTCCCGTACCCGGCGTTTTGTTTTTTCGGGTGATGATATTGAAGCTCGGCTGttt                                    >  1:950280/1‑150 (MQ=255)
                                                                                                                 ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATTGGAAGTGAGCGCAGAGAGGGGTTTTTTTTCCTTACTGAACTGGCGGTTAAGTTTTCGCGGAACCGGCGTTTTGTTTTTTCGGGTGATGATATTTAAGGCCGGCTgggt                                    >  2:142122/1‑146 (MQ=255)
                                                                                                                       cGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGtt                                <  1:583621/147‑1 (MQ=255)
                                                                                                                        gTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGGGAAGGGAGCGCCGGGAGGGGTTATTTTTTCTGGCTGCAATTGCGATTAAGGTTCCGCGGACCCGGCGGTTTGTTATTGCGGGTGGTGATATTGAAAGCCGGGGGGTTGGTTCAtt                            >  2:98539/1‑150 (MQ=255)
                                                                                                                               gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt                    >  2:334714/1‑151 (MQ=255)
                                                                                                                                cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc                   >  2:354876/1‑151 (MQ=255)
                                                                                                                                cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc                   >  1:119400/1‑151 (MQ=255)
                                                                                                                                   tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAc                 >  1:102345/1‑150 (MQ=255)
                                                                                                                                   tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTGCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAc                 <  2:338398/150‑1 (MQ=255)
                                                                                                                                     aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGc              <  2:1177367/151‑1 (MQ=255)
                                                                                                                                        tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt             <  2:134903/150‑1 (MQ=255)
                                                                                                                                         ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg           >  1:1166635/1‑151 (MQ=255)
                                                                                                                                          tttctttttctGCCATGAGAAGAGAGCGCAGCGAGGGGTTATTTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGCTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAACGCAgg           >  2:1068860/1‑150 (MQ=255)
                                                                                                                                            tctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc         <  1:215273/150‑1 (MQ=255)
                                                                                                                                               ttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtg     >  2:929646/1‑151 (MQ=255)
                                                                                                                                                  tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAt  >  2:422594/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
AAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAT  >  NC_000913/4604204‑4604499

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