Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT46.9% 48.5 / 45.7 65S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (7/27);  new base T (30/0);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 7.46e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604198‑4604492
                                                                                                                                                   |                                                                                                                                                     
cATAGAAAAATGCGTCATGAGTAGTATCGGGATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTGTACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑Ttttttctttt                                                                                                                                                   <  1:748277/151‑1 (MQ=255)
     aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCC‑Ttttttcttttttt                                                                                                                                                >  2:95860/1‑147 (MQ=255)
      aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGCTTGCGGCGGTCATTTTTTTCTTTTTTTGc                                                                                                                                              >  1:125548/1‑151 (MQ=255)
       aaaTGCGTCATGAGTAGAATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCA‑tttttttttttttttcc                                                                                                                                             >  1:1092325/1‑145 (MQ=255)
          tGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGcc                                                                                                                                             >  1:296362/1‑147 (MQ=255)
           gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATgag                                                                                                                                        <  1:1386286/151‑1 (MQ=255)
           gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATga                                                                                                                                         <  2:1109993/150‑1 (MQ=255)
             gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGaa                                                                                                                                      <  1:1016099/151‑1 (MQ=255)
               cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑tttttttttttttttccctgcgga                                                                                                                                      >  1:725393/1‑137 (MQ=255)
               cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCCTTAGaa                                                                                                                                      >  2:369140/1‑149 (MQ=255)
                aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAg                                                                                                                                     <  1:760390/149‑1 (MQ=255)
                aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTATTCCATgatgagcg                                                                                                                                   >  1:621780/1‑145 (MQ=255)
                 tGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCT‑TTTTTTTTTTTTTTTCCATGAGaa                                                                                                                                      >  2:560215/1‑147 (MQ=255)
                  gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCATTAGGAggggg                                                                                                                                 >  1:929631/1‑147 (MQ=255)
                  gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGa                                                                                                                                  <  2:1360372/150‑1 (MQ=255)
                        taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCT‑TTTTTTTTTTTTTTTCCATGAGGAGTGAGGGCag                                                                                                                            >  1:1059568/1‑150 (MQ=255)
                           cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCCTTGGAAGTGGGGGCcgag                                                                                                                          >  2:840293/1‑149 (MQ=255)
                            ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAgg                                                                                                                       >  1:476123/1‑151 (MQ=255)
                               aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCGGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCG‑TTTTTTTTTTTTTTGCCATGAGAAGGGGGGGCAGAGAggg                                                                                                                      >  2:1303577/1‑149 (MQ=255)
                               aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAC‑TTTTTTTTTTTTCTTCCATGAGAAGTGAGCCCAGAGAggg                                                                                                                      >  1:7905/1‑149 (MQ=255)
                                  gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGggggggttt                                                                                                                 >  2:449405/1‑146 (MQ=255)
                                  gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑t                                                                                                                <  1:1123726/151‑1 (MQ=255)
                                    gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCG‑TTTTTTTTTTTTTTGGCGTTTGAAGGGAGGGCCGAGggggggt                                                                                                                   >  2:614902/1‑144 (MQ=255)
                                     agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTTTTCCCTGAGAAGGGAGGGCAGAGGGGGATTA‑TCtt                                                                                                             >  1:282437/1‑151 (MQ=255)
                                              aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGTTCT‑GTTTTTAATAATATACCATGGGAAGGGAGGGGGGAGAGGGAGTA‑TCTTGCcttact                                                                                                     >  1:893352/1‑150 (MQ=255)
                                                aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGaa                                                                                                  <  2:587972/151‑1 (MQ=255)
                                                 aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGAGTA‑TCTTGCCTGACTGaa                                                                                                  >  2:1156315/1‑150 (MQ=255)
                                                    tttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTg                                                                                               <  1:1315806/150‑1 (MQ=255)
                                                     ttGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGTATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTgg                                                                                              <  1:1374642/150‑1 (MQ=255)
                                                         gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTCTGCCCTGAGAAGTGGGGGCAGAGGGGGGTTA‑TTTTTCcttac                                                                                                      >  2:512481/1‑138 (MQ=255)
                                                          ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTa                                                                                        <  2:784056/151‑1 (MQ=255)
                                                              aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTTCCTGAGAAGAGAGCGCAGCGAGGGGGTT‑TCTTGCctgactggattcggggtt                                                                                         >  2:28760/1‑135 (MQ=255)
                                                              aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCCTTGGAAAGGAGAGCCGAGAGGGATTATTTTTGCCCG‑CTGAAGTGGGGATTTAggt                                                                                    >  1:942388/1‑149 (MQ=255)
                                                               cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGttt                                                                                   <  1:143084/151‑1 (MQ=255)
                                                                  ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCCTGAGGAGTGGGGGCAGAGAGGGATTT‑TTTTGCCTGACTGAACTGGGGATTTAGgtttc                                                                                 >  1:1277119/1‑148 (MQ=255)
                                                                   tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                                 <  2:1064009/149‑1 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCC‑TTTTTTTTTTTTTTTCCATGAGGAGTGTGCGCAGAGAGGGGGTA‑TTTTTTCTTATTCAACTTGGGATTAAGTTTCTTCGTAccc                                                                         >  2:31904/1‑151 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTA‑TCTTGCCTGACTGAACTGGGGATTTAGTTTCCGCGTAccc                                                                         >  1:930700/1‑151 (MQ=255)
                                                                         cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCC‑TTTTTTTTTTTTTTTTCAAGAGGAGTGAGGGGAGGGAGGGGTTT‑TTTTTCCTTACTCCAATGGCCGTTAAGTTTCCCCGTAcc                                                                          >  2:948153/1‑150 (MQ=255)
                                                                           ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                                         <  2:1164161/149‑1 (MQ=255)
                                                                           ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg                                                                       <  1:716800/151‑1 (MQ=255)
                                                                            cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGAGGGGTTA‑TCTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGGACCCgg                                                                       >  1:394293/1‑150 (MQ=255)
                                                                              tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCCTGAGAAGAGAGGGCAGGGGGGGGGTT‑TTTTTCCTGACTGAACTGGCGGTTAAGTTTTCCCGTAACCGGCg                                                                     >  2:936600/1‑150 (MQ=255)
                                                                                 cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCC‑TTTTTTTTTTTTTTGGCATTAGAAGTGGGGGCAGAGGGGGGTTA‑TCTTGCCTGAATGAACTAGCGGTTAAGGTTCCGCGTACCCGGCGgtt                                                                  >  1:763663/1‑150 (MQ=255)
                                                                                   ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg                                                                <  2:466723/150‑1 (MQ=255)
                                                                                      ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTA‑TCTTGTCTGACTGAACTGGCGGTTTAGTTTTCGCGGACCCCGCGTTTTGtttt                                                            >  2:1354601/1‑149 (MQ=255)
                                                                                       cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                           <  1:613185/151‑1 (MQ=255)
                                                                                          tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGc                                                         <  2:273007/150‑1 (MQ=255)
                                                                                             ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCT‑TTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCCGAGAGGGGTTA‑TCTTTCCTGAATGAACTTGGGGTTAAGTTTTCCCGTAACCCGCGCTTTGGTTTTGCgg                                                       >  1:654615/1‑149 (MQ=255)
                                                                                             ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                                      <  1:186067/150‑1 (MQ=255)
                                                                                             ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCA‑ATTTTTTTTTTTTTGGCATGGGAAGTGAGGGCCGGGGGGGGTTT‑TTTTTGCTGACTGAAATGGCGATTTAGTTTTCGCGTACCCGGCGGTTTGTTTTTGCgg                                                       >  2:5605/1‑149 (MQ=255)
                                                                                                gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                                   <  2:621765/150‑1 (MQ=255)
                                                                                                         ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGaa                                         <  2:597189/151‑1 (MQ=255)
                                                                                                                aTGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGGTTAAGTTTCCGCGTACCCGGCGTTTTGTTTTTGCGGGTGATGATATTTAAGCTCgg                                   >  1:180377/1‑150 (MQ=255)
                                                                                                                aTGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCGTTTTGTTTTTGCGGGTGGTGGTATTTAAACTCgg                                   >  1:566392/1‑150 (MQ=255)
                                                                                                                aTGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGc                                  <  2:1311294/151‑1 (MQ=255)
                                                                                                                  gTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCt                                 <  2:1033370/150‑1 (MQ=255)
                                                                                                                    cgcgGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGa                               <  2:48032/150‑1 (MQ=255)
                                                                                                                     gcgGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGAtt                             <  1:1247989/151‑1 (MQ=255)
                                                                                                                      cgGATCTCGTTTGCGGCGGTCC‑TTTTTTTTTTTTTTGCCCTGTGAAGTGGGCGCAGAGAGGGGTTA‑TCTTGCCTGACTGAACTGGCGGTTAAGTTTTCCCGTACCCGGCGTTTTGTTTTTTCGGGTGATGATATTGAAGCTCGGCTGttt                             >  1:950280/1‑150 (MQ=255)
                                                                                                                       ggATCTCGTTTGCGGCGGTCA‑TTTTTTTTTTTTTTGCCATTGGAAGTGAGCGCAGAGAGGGGTTT‑TTTTTCCTTACTGAACTGGCGGTTAAGTTTTCGCGGAACCGGCGTTTTGTTTTTTCGGGTGATGATATTTAAGGCCGGCTgggt                             >  2:142122/1‑146 (MQ=255)
                                                                                                                             cGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGtt                         <  1:583621/147‑1 (MQ=255)
                                                                                                                              gTTTGCGGCGGTCA‑TTTTTTTTTTTTCTGCCCTGGGAAGGGAGCGCCGGGAGGGGTTA‑TTTTTTCTGGCTGCAATTGCGATTAAGGTTCCGCGGACCCGGCGGTTTGTTATTGCGGGTGGTGATATTGAAAGCCGGGGGGTTGGTTCAtt                     >  2:98539/1‑150 (MQ=255)
                                                                                                                                     gcggTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt             >  2:334714/1‑151 (MQ=255)
                                                                                                                                      cggTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc            >  2:354876/1‑151 (MQ=255)
                                                                                                                                      cggTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc            >  1:119400/1‑151 (MQ=255)
                                                                                                                                         tCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAc          >  1:102345/1‑150 (MQ=255)
                                                                                                                                         tCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTGCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAc          <  2:338398/150‑1 (MQ=255)
                                                                                                                                           a‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGc       <  2:1177367/151‑1 (MQ=255)
                                                                                                                                              tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt      <  2:134903/150‑1 (MQ=255)
                                                                                                                                               ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg    >  1:1166635/1‑151 (MQ=255)
                                                                                                                                                tttctttttctGCCATGAGAAGAGAGCGCAGCGAGGGGTTA‑TTTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGCTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAACGCAgg    >  2:1068860/1‑150 (MQ=255)
                                                                                                                                                  tctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc  <  1:215273/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604198‑4604492

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: ‑