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

AAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAT  >  NC_000913/4604204‑4604499
                                                                                                                                                   |                                                                                                                                                      
aaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGCTTGCGGCGGTCATTTTTTTCTTTTTTTGc                                                                                                                                                     >  1:125548/1‑151 (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‑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)
                                                                                                                                               ttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGgtg     >  2:929646/1‑151 (MQ=255)
                                                                                                                                                  tctGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAt  >  2:422594/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
AAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA‑TTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTA‑TCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAT  >  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: ‑