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.8% 367.6 / ‑1.4 82intergenic (‑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 . (34/46);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAG  >  NC_000913/4556611‑4556903
                                                                                                                                                      |                                                                                                                                              
tATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGt                                                                                                                                                >  2:391025/1‑151 (MQ=255)
   cAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttt                                                                                                                                             <  2:686272/151‑1 (MQ=255)
   cAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttt                                                                                                                                              <  2:521181/150‑1 (MQ=255)
        tataTTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttttttt                                                                                                                                         >  2:929945/1‑150 (MQ=255)
         atatTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttg                                                                                                                                        <  1:16853/150‑2 (MQ=255)
         atatTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttg                                                                                                                                        >  2:789366/1‑149 (MQ=255)
         atatTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttttttt                                                                                                                                         >  1:938351/1‑149 (MQ=255)
           attattTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttttgag                                                                                                                                      <  2:4828/150‑4 (MQ=255)
            ttattTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGtt                                                                                                                                   >  1:785849/1‑151 (MQ=255)
              attTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGttt                                                                                                                                  <  2:51561/150‑1 (MQ=255)
               tttACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTtg                                                                                                                                 >  1:881892/1‑150 (MQ=255)
               tttACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTtg                                                                                                                                 >  1:527738/1‑150 (MQ=255)
               tttACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACAGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTtg                                                                                                                                 >  2:573973/1‑150 (MQ=255)
                ttACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTtga                                                                                                                                <  2:611996/150‑1 (MQ=255)
                  aCAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTtgatg                                                                                                                              >  1:325668/1‑151 (MQ=255)
                     aaCTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGAttt                                                                                                                          >  2:204560/1‑151 (MQ=255)
                     aaCTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGAttt                                                                                                                          >  1:131469/1‑151 (MQ=255)
                      aCTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGAtttt                                                                                                                         <  1:666177/151‑1 (MQ=255)
                      aCTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGAtttt                                                                                                                         >  1:22280/1‑151 (MQ=255)
                         aaTTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTAt                                                                                                                      <  1:504139/151‑1 (MQ=255)
                               ttCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAAt                                                                                                                 <  1:132122/150‑1 (MQ=255)
                                tCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATc                                                                                                                >  2:923038/1‑150 (MQ=255)
                                  aaaTAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATg                                                                                                             <  2:261908/151‑1 (MQ=255)
                                      aataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACt                                                                                                          <  2:262833/150‑1 (MQ=255)
                                      aataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACt                                                                                                          <  1:880207/150‑1 (MQ=255)
                                       ataCATGGCTGATTATGCGGAAATAAATAACTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATGCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACttt                                                                                                        <  2:59284/151‑1 (MQ=255)
                                        taCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTa                                                                                                                       >  2:84900/1‑135 (MQ=255)
                                        taCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTa                                                                                                                       <  1:84911/135‑1 (MQ=255)
                                         aCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACttttt                                                                                                      >  2:117221/1‑151 (MQ=255)
                                            tGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTcc                                                                                                   >  1:398793/1‑151 (MQ=255)
                                              gCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTg                                                                                                 <  2:723827/151‑1 (MQ=255)
                                                 gATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaa                                                                                               <  1:635509/150‑1 (MQ=255)
                                                      tGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTc                                                                                          >  2:759880/1‑150 (MQ=255)
                                                       gCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAg                                                                                        <  2:838884/151‑1 (MQ=255)
                                                         ggAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTa                                                                                      <  1:636667/151‑1 (MQ=255)
                                                          gAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAg                                                                                     <  1:157959/151‑1 (MQ=255)
                                                           aaataaatAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGa                                                                                    <  2:957020/151‑1 (MQ=255)
                                                            aataaataATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGat                                                                                   <  2:786278/151‑1 (MQ=255)
                                                            aataaataATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGat                                                                                   <  1:497474/151‑1 (MQ=255)
                                                                 ataatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTc                                                                               >  1:480188/1‑150 (MQ=255)
                                                                   aatTTCCCACCGGAATTAAGTAGCAGCTGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCAt                                                                             >  1:384072/1‑150 (MQ=255)
                                                                     tttCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAgg                                                                          >  1:202980/1‑151 (MQ=255)
                                                                        cccACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCaa                                                                       <  1:640447/151‑1 (MQ=255)
                                                                        cccACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCa                                                                        >  2:973703/1‑150 (MQ=255)
                                                                          cACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTGTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGt                                                                     <  1:670322/151‑1 (MQ=255)
                                                                             cGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTa                                                                    >  2:336533/1‑149 (MQ=255)
                                                                                aaTTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGtt                                                               >  1:464914/1‑151 (MQ=255)
                                                                                    aaGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTtat                                                           <  1:392135/151‑1 (MQ=255)
                                                                                     aGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTtata                                                          <  2:59024/151‑1 (MQ=255)
                                                                                     aGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTtata                                                          <  1:245534/151‑1 (MQ=255)
                                                                                        agcagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACtt                                                       >  1:359769/1‑151 (MQ=255)
                                                                                        agcagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACtt                                                       <  1:342411/151‑1 (MQ=255)
                                                                                        agcagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACt                                                        <  2:842580/150‑1 (MQ=255)
                                                                                         gcagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACttt                                                      <  1:954554/151‑1 (MQ=255)
                                                                                          cagcGGTTACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTg                                                     <  2:887490/151‑1 (MQ=255)
                                                                                          cagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACttt                                                      >  2:341472/1‑150 (MQ=255)
                                                                                          cagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTg                                                     <  2:131458/151‑1 (MQ=255)
                                                                                           agcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGc                                                    <  1:292631/151‑1 (MQ=255)
                                                                                            gcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCt                                                   <  1:760276/151‑1 (MQ=255)
                                                                                               gTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGAc                                                <  2:31158/151‑1 (MQ=255)
                                                                                                tGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGAc                                                >  1:558463/1‑150 (MQ=255)
                                                                                                 gACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGAc                                                >  1:894246/1‑149 (MQ=255)
                                                                                                     aGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGAtt                                          >  2:530464/1‑151 (MQ=255)
                                                                                                     aGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGAtt                                          >  1:82947/1‑151 (MQ=255)
                                                                                                      gTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTc                                         <  1:930031/151‑1 (MQ=255)
                                                                                                          tttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAgg                                      <  1:922595/150‑1 (MQ=255)
                                                                                                             cATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTg                                  >  1:469498/1‑151 (MQ=255)
                                                                                                                cTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtct                               >  2:633965/1‑151 (MQ=255)
                                                                                                                cTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtc                                >  2:911112/1‑150 (MQ=255)
                                                                                                                   cGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtctc                              >  1:371565/1‑149 (MQ=255)
                                                                                                                       aCTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTg                        <  1:933889/151‑1 (MQ=255)
                                                                                                                            tCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACAt                    <  1:215767/150‑1 (MQ=255)
                                                                                                                            tCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACAt                    <  1:532695/150‑1 (MQ=255)
                                                                                                                             cGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATg                   <  1:788159/150‑1 (MQ=255)
                                                                                                                               gAATATTCAGAATATACTAGC‑GATTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGAt                >  1:92096/1‑151 (MQ=255)
                                                                                                                                 atatTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGt              <  2:938269/151‑1 (MQ=255)
                                                                                                                                   atTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTa             >  1:588243/1‑150 (MQ=255)
                                                                                                                                   atTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTa             <  1:320867/150‑1 (MQ=255)
                                                                                                                                    ttCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTaa            <  2:429508/150‑1 (MQ=255)
                                                                                                                                    ttCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAAt           >  1:282826/1‑151 (MQ=255)
                                                                                                                                     tCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAAt           <  1:67541/150‑1 (MQ=255)
                                                                                                                                      cAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGa         <  2:600723/151‑1 (MQ=255)
                                                                                                                                       aGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGaa        <  2:626446/151‑1 (MQ=255)
                                                                                                                                        gAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAAc       >  2:529225/1‑151 (MQ=255)
                                                                                                                                         aaTATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAAc       >  1:350842/1‑150 (MQ=255)
                                                                                                                                         aaTATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACa      <  1:759965/151‑1 (MQ=255)
                                                                                                                                           tataCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACaa     >  2:282391/1‑150 (MQ=255)
                                                                                                                                            ataCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACaaa    >  1:666782/1‑150 (MQ=255)
                                                                                                                                              aCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAg  <  1:300281/150‑1 (MQ=255)
                                                                                                                                              aCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAg  <  2:527683/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
TATCAATTTATATTATTTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAG  >  NC_000913/4556611‑4556903

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