Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,556,769 (T)9→8 intergenic (‑449/‑224) yjiC ← / → iraD uncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,556,7610T.98.7% 356.6 / ‑3.8 77intergenic (‑441/‑232)yjiC/iraDuncharacterized protein/RpoS stabilzer after DNA damage, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (30/46);  total (31/46)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGA  >  NC_000913/4556612‑4556906
                                                                                                                                                     |                                                                                                                                                 
aTCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtt                                                                                                                                                  <  1:228340/151‑1 (MQ=255)
   aaTTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCTGTGACAAGTATTTTCATCTACGTAACTATTCGTAATATTCAGAATATGCTAGCTGttt                                                                                                                                                 <  2:760226/149‑1 (MQ=255)
    aTTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttttt                                                                                                                                              <  2:193961/151‑1 (MQ=255)
    aTTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttttt                                                                                                                                              <  1:819702/151‑1 (MQ=255)
        atatTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAGATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttg                                                                                                                                           <  1:1212389/150‑2 (MQ=255)
        atatTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttga                                                                                                                                          <  2:868405/151‑3 (MQ=255)
         tattatTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttg                                                                                                                                           >  1:183262/1‑148 (MQ=255)
            tattTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGttt                                                                                                                                     >  1:1178272/1‑151 (MQ=255)
                tACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTtgat                                                                                                                                  <  1:1125569/150‑1 (MQ=255)
                tACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTtgat                                                                                                                                  <  2:369612/150‑1 (MQ=255)
                 aCAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTtgatg                                                                                                                                 >  1:1214204/1‑150 (MQ=255)
                  cAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTgatgat                                                                                                                               <  2:565749/151‑1 (MQ=255)
                     aCTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGAtttt                                                                                                                            <  1:707622/151‑1 (MQ=255)
                              ttCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAAt                                                                                                                    >  1:145552/1‑150 (MQ=255)
                              ttCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAAt                                                                                                                    <  1:1224550/150‑1 (MQ=255)
                                cAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCAt                                                                                                                 >  2:368830/1‑151 (MQ=255)
                                 aaaTAATACATGGCTGATTATGCGGAAATAAATCATTTCCCTCCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATg                                                                                                                <  2:603034/151‑1 (MQ=255)
                                   ataataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGAc                                                                                                              >  2:116822/1‑151 (MQ=255)
                                     aataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACtt                                                                                                            <  2:392124/151‑1 (MQ=255)
                                      ataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACttt                                                                                                           >  2:36946/1‑151 (MQ=255)
                                      ataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACttt                                                                                                           <  2:398502/151‑1 (MQ=255)
                                                gATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaa                                                                                                 <  1:1014389/151‑1 (MQ=255)
                                                gATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaa                                                                                                  >  2:713505/1‑150 (MQ=255)
                                                gATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaa                                                                                                  <  1:930014/150‑1 (MQ=255)
                                                gATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaa                                                                                                  <  2:496580/150‑1 (MQ=255)
                                                 aTTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaaa                                                                                                <  2:1165324/151‑1 (MQ=255)
                                                   tATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAg                                                                                               <  1:368841/150‑1 (MQ=255)
                                                        ggAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTa                                                                                         <  1:101098/151‑1 (MQ=255)
                                                             taaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGatat                                                                                    >  2:854133/1‑151 (MQ=255)
                                                              aaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGatat                                                                                    <  2:28148/150‑1 (MQ=255)
                                                               aataatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTc                                                                                  <  1:542766/151‑1 (MQ=255)
                                                                   atTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAg                                                                              <  2:1234696/151‑1 (MQ=255)
                                                                     ttCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAGTATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATa                                                                               >  2:349429/1‑148 (MQ=255)
                                                                       cccACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCaa                                                                          <  2:1211224/151‑1 (MQ=255)
                                                                       cccACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCa                                                                           <  2:827446/150‑1 (MQ=255)
                                                                       cccACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCa                                                                           <  1:829820/150‑1 (MQ=255)
                                                                         cACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGt                                                                        >  2:191914/1‑151 (MQ=255)
                                                                           ccGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTa                                                                       >  2:9061/1‑150 (MQ=255)
                                                                             ggAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAgg                                                                    <  1:368186/151‑1 (MQ=255)
                                                                              gAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGt                                                                   >  1:287775/1‑151 (MQ=255)
                                                                                 ttAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGtttt                                                                >  2:411260/1‑151 (MQ=255)
                                                                                  tAAGTAGCAGCCGCGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGtttt                                                                >  1:580923/1‑150 (MQ=255)
                                                                                     gTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGGTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTtat                                                              >  1:894472/1‑150 (MQ=255)
                                                                                     gTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATAc                                                            <  1:933421/151‑1 (MQ=255)
                                                                                     gTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATAc                                                            <  1:824065/151‑1 (MQ=255)
                                                                                        gcagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACttt                                                         <  2:618172/151‑1 (MQ=255)
                                                                                            cGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCt                                                      <  2:112774/150‑1 (MQ=255)
                                                                                                gACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACa                                                  <  2:77709/150‑1 (MQ=255)
                                                                                                gACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAg                                                 <  1:703542/151‑1 (MQ=255)
                                                                                                gACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAg                                                 <  1:336601/151‑1 (MQ=255)
                                                                                                  cAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAgg                                                >  1:1121137/1‑150 (MQ=255)
                                                                                                  cAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAgg                                                >  2:528103/1‑150 (MQ=255)
                                                                                                   aaGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGAt                                              <  2:429220/151‑1 (MQ=255)
                                                                                                    aGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGAtt                                             >  1:776114/1‑151 (MQ=255)
                                                                                                    aGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGAtt                                             >  1:452904/1‑151 (MQ=255)
                                                                                                     gTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTc                                            <  1:398919/151‑1 (MQ=255)
                                                                                                         tttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGATAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAgt                                         <  2:549591/150‑2 (MQ=255)
                                                                                                         tttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAgg                                         <  1:537949/150‑1 (MQ=255)
                                                                                                         tttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAgg                                         <  1:164018/150‑1 (MQ=255)
                                                                                                          ttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGcc                                       >  1:372831/1‑151 (MQ=255)
                                                                                                          ttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGc                                        <  2:114476/150‑1 (MQ=255)
                                                                                                            cATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTg                                     <  2:1117468/151‑1 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtctc                                 <  1:191901/150‑1 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtctc                                 <  1:218723/150‑1 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCa                                <  2:932890/151‑1 (MQ=255)
                                                                                                                 aCGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCa                                <  1:590083/151‑1 (MQ=255)
                                                                                                                   gTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtctc                                 >  1:171257/1‑148 (MQ=255)
                                                                                                                     aaCTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACt                            >  2:1064123/1‑151 (MQ=255)
                                                                                                                      aCTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTg                           <  2:686140/151‑1 (MQ=255)
                                                                                                                         aTTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGa                          <  2:1038964/149‑1 (MQ=255)
                                                                                                                          ttCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACAt                       >  2:581001/1‑151 (MQ=255)
                                                                                                                                atatTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGt                 <  1:528124/151‑1 (MQ=255)
                                                                                                                                   ttCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTaa               >  2:762489/1‑150 (MQ=255)
                                                                                                                                    tCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATg             >  2:971003/1‑151 (MQ=255)
                                                                                                                                     cAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATg             <  2:43559/150‑1 (MQ=255)
                                                                                                                                     cAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATg             >  2:1228242/1‑150 (MQ=255)
                                                                                                                                      aGAATATACTAGCGG‑TTTTTTTTGAGATTGATGATTTTTATCCAACCATGACTTTTTTCCTGAAAAGTCAGTAGAGATTCATAGCCAAGTAAGGTTTTATACTTTGCGGACAGGATTCAGGCCTGTCTCAGACTGACATGTATGTAATg             >  2:597826/1‑149 (MQ=255)
                                                                                                                                         atatACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACaa        <  1:261730/151‑1 (MQ=255)
                                                                                                                                         atatACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACaa        <  2:1228613/151‑1 (MQ=255)
                                                                                                                                          tataCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACaa        >  1:269978/1‑150 (MQ=255)
                                                                                                                                           ataCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACaaa       >  2:75839/1‑150 (MQ=255)
                                                                                                                                           ataCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACaaa       >  1:1183803/1‑150 (MQ=255)
                                                                                                                                             aCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAgg    <  1:1204719/151‑1 (MQ=255)
                                                                                                                                                aGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGa  >  2:284332/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGA  >  NC_000913/4556612‑4556906

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