Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 13.4% intergenic (+34/‑183) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2640AC13.4% 329.7 / 43.6 119intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (53/50);  new base C (9/7);  total (62/57)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 121 total aligned reads displayed.

GTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGA  >  NC_000913/4060116‑4060412
                                                                                                                                                    |                                                                                                                                                     
    aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                   >  1:218848/1‑149 (MQ=255)
    aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                   >  3:631672/1‑149 (MQ=255)
     tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                                                                  >  2:153208/1‑149 (MQ=255)
      cGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAg                                                                                                                                                 <  4:761595/149‑1 (MQ=255)
      cGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                                 >  3:592293/1‑149 (MQ=255)
         tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                                                                             <  4:352927/149‑1 (MQ=255)
          gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttt                                                                                                                                            <  3:668829/149‑1 (MQ=255)
           atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                                                                           <  1:654646/149‑1 (MQ=255)
           atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCATTTTTTAGTACTAAAAATCCTTCCATGGCTGTCAGGtttt                                                                                                                                           <  2:847709/149‑1 (MQ=255)
            tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttt                                                                                                                                          <  1:42421/149‑1 (MQ=255)
            tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                                                                          >  3:166576/1‑149 (MQ=255)
            tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                                                                          >  1:314517/1‑149 (MQ=255)
             gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTAACGGAAAACGATCGTCGCCCCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                                                                         <  4:588868/149‑1 (MQ=255)
              aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                                                                        <  4:837622/149‑2 (MQ=255)
               cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                                                                       <  1:887673/149‑1 (MQ=255)
                gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                                                                      >  1:750216/1‑149 (MQ=255)
                gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                                                                      >  3:800564/1‑149 (MQ=255)
                  aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                                                                                    <  1:217009/149‑1 (MQ=255)
                         gAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGt                                                                                                                             <  1:350105/149‑1 (MQ=255)
                           aGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAATACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGa                                                                                                                           <  1:523137/149‑1 (MQ=255)
                             tGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAg                                                                                                                          >  3:128994/1‑149 (MQ=255)
                                 tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                                                                                    <  1:341375/134‑1 (MQ=255)
                                 tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                     >  3:68180/1‑149 (MQ=255)
                                 tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                     >  3:452869/1‑149 (MQ=255)
                                 tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcac                                                                                                                     >  3:236380/1‑149 (MQ=255)
                                       ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                                                                               <  3:233718/149‑1 (MQ=255)
                                         ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                                                                             >  1:401107/1‑149 (MQ=255)
                                         ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                                                                             <  2:314517/149‑1 (MQ=255)
                                          tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCGTGTCGTGATCGCCCACCAATGCg                                                                                                            >  4:280488/1‑149 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                           >  3:356726/1‑149 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                           >  3:641113/1‑149 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                                                                           >  3:89190/1‑149 (MQ=255)
                                            tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                          >  4:513996/1‑149 (MQ=255)
                                             ccGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                                                                         >  4:177587/1‑149 (MQ=255)
                                              cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                                                                        <  2:723030/149‑1 (MQ=255)
                                                 aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                                                                     <  2:331758/149‑1 (MQ=255)
                                                   tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                                                                   <  2:232347/149‑1 (MQ=255)
                                                    cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                                                                  <  3:348739/149‑1 (MQ=255)
                                                     gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                                                                 >  3:643282/1‑149 (MQ=255)
                                                         aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGa                                                                                                                           >  4:389792/1‑119 (MQ=255)
                                                         aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGa                                                                                                                           <  3:389792/119‑1 (MQ=255)
                                                                cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                                                                      <  4:479705/149‑1 (MQ=255)
                                                                cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                                                                      <  1:683809/149‑1 (MQ=255)
                                                                 tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                     <  4:620076/149‑1 (MQ=255)
                                                                     ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                                 >  1:330354/1‑149 (MQ=255)
                                                                      gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                                                                <  3:318998/149‑1 (MQ=255)
                                                                      gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                                                                <  4:839397/149‑1 (MQ=255)
                                                                      gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                                                                <  1:272619/149‑1 (MQ=255)
                                                                        aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGt                                                                              <  4:356726/149‑1 (MQ=255)
                                                                          aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTTGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                                                                 <  3:363120/81‑1 (MQ=255)
                                                                             aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAACACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                         <  4:631672/149‑1 (MQ=255)
                                                                             aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                                                                          <  4:128994/149‑1 (MQ=255)
                                                                             aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                         <  2:280352/149‑1 (MQ=255)
                                                                                 tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                                                                     <  1:471142/149‑1 (MQ=255)
                                                                                  cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                                                                    <  2:574982/149‑1 (MQ=255)
                                                                                    ccgcGCCAACCCCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                                                                  >  3:692684/1‑149 (MQ=255)
                                                                                      gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                                                                <  3:417631/149‑1 (MQ=255)
                                                                                      gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                                                                >  3:678029/1‑149 (MQ=255)
                                                                                      gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                                                                <  3:468963/149‑1 (MQ=255)
                                                                                      gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                                                                <  4:141349/149‑1 (MQ=255)
                                                                                           aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac                                                           <  4:838535/149‑1 (MQ=255)
                                                                                             ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                         <  3:197535/149‑1 (MQ=255)
                                                                                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                                   >  2:767976/1‑149 (MQ=255)
                                                                                                     cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGc                                                 >  3:113841/1‑149 (MQ=255)
                                                                                                         aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAg                                             <  2:198998/149‑1 (MQ=255)
                                                                                                          gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGt                                            >  1:168640/1‑149 (MQ=255)
                                                                                                           aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTc                                           <  3:376638/149‑1 (MQ=255)
                                                                                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGGACACCAATGCGAGTTATTCCCTTATAAGTCTGTGGTTACGTTAATAGAGATCTTTTTCACAATTATCTGAGAAGTCAt                                         >  2:614738/1‑149 (MQ=255)
                                                                                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCAt                                         >  1:41111/1‑149 (MQ=255)
                                                                                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACTAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTTAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATa                                        >  4:248960/1‑149 (MQ=255)
                                                                                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCGGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATa                                        >  1:658124/1‑149 (MQ=255)
                                                                                                               ttAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATAc                                       >  4:561094/1‑149 (MQ=255)
                                                                                                                  aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGcc                                    >  3:675568/1‑149 (MQ=255)
                                                                                                                   tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCg                                   >  4:627307/1‑149 (MQ=255)
                                                                                                                       gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTaa                               >  2:630256/1‑149 (MQ=255)
                                                                                                                        ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAt                              <  4:692684/149‑1 (MQ=255)
                                                                                                                         tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGACACAATTATCTGTAAAGTCATACGCCGTTAAtt                             >  1:286692/1‑149 (MQ=255)
                                                                                                                            ctAGTCATAAAACCCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCTTTAGCGCACACCAGTGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTAAATAGAGAGCTTTGTCACAATTCTCTGCAAAGTCATACGACGTTAAGTGct                          <  2:780331/148‑1 (MQ=255)
                                                                                                                              aGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGcttg                        <  3:512406/149‑2 (MQ=255)
                                                                                                                                tCATAAGAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttct                      >  2:633520/1‑149 (MQ=255)
                                                                                                                                tCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttct                      >  1:273933/1‑149 (MQ=255)
                                                                                                                                tCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTATGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGACGTTAATTGctttct                      >  3:782563/1‑149 (MQ=255)
                                                                                                                                 cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttctt                     <  1:281047/149‑1 (MQ=255)
                                                                                                                                  aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttcttt                    <  1:271497/149‑1 (MQ=255)
                                                                                                                                    aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCttttt                  <  1:302230/149‑1 (MQ=255)
                                                                                                                                     aaaacCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt                 <  1:116198/149‑1 (MQ=255)
                                                                                                                                     aaaaaCCTGCCATCGCTGGCAGGTTTTGTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt                 >  3:782942/1‑149 (MQ=255)
                                                                                                                                     aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt                 >  4:294543/1‑149 (MQ=255)
                                                                                                                                       aaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg               >  1:221835/1‑149 (MQ=255)
                                                                                                                                       aaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg               >  3:857910/1‑149 (MQ=255)
                                                                                                                                        aaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                                                                     <  4:40690/62‑1 (MQ=255)
                                                                                                                                        aaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc              >  2:574191/1‑149 (MQ=255)
                                                                                                                                         aCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcg             >  2:106582/1‑149 (MQ=255)
                                                                                                                                         aCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcg             >  1:407888/1‑149 (MQ=255)
                                                                                                                                          ccTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgt            >  2:697338/1‑149 (MQ=255)
                                                                                                                                           cTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgta           >  3:426125/1‑149 (MQ=255)
                                                                                                                                            tGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaa          >  2:823644/1‑149 (MQ=255)
                                                                                                                                            tGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaa          <  3:205056/149‑1 (MQ=255)
                                                                                                                                             gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaag         <  3:771039/149‑1 (MQ=255)
                                                                                                                                                    aTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGa  <  2:279547/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGA  >  NC_000913/4060116‑4060412

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: