Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,261,359 T→A 5.2% intergenic (+53/‑310) yjbM → / → dusA uncharacterized protein YjbM/tRNA‑dihydrouridine synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,261,3590TA 5.2% 1539.6 / 42.4 459intergenic (+53/‑310)yjbM/dusAuncharacterized protein YjbM/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (230/205);  new base A (12/12);  total (242/217)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 472 total aligned reads displayed.

GAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTC  >  NC_000913/4261210‑4261508
                                                                                                                                                     |                                                                                                                                                     
  aGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtattta                                                                                                                                                     >  2:572193/1‑150 (MQ=255)
   gTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAAtatttat                                                                                                                                                    >  2:1401907/1‑150 (MQ=255)
      gTTTGACCCAACCTCAAAGGAAGATATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTgg                                                                                                                                                 <  2:902167/150‑1 (MQ=255)
          gACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGtcc                                                                                                                                             <  1:164563/150‑3 (MQ=255)
           aCCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCt                                                                                                                                            >  1:106963/1‑150 (MQ=255)
           aCCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCt                                                                                                                                            <  2:938821/150‑1 (MQ=255)
            cccAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtt                                                                                                                                           >  1:1400288/1‑150 (MQ=255)
            cccAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTAtt                                                                                                                                           >  2:677784/1‑150 (MQ=255)
                aCCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                                                                                       <  1:61847/150‑3 (MQ=255)
                aCCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGAACCCAATAAATATTGGGTCCttttat                                                                                                                                       <  2:1383635/150‑3 (MQ=255)
                 ccTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCtttttta                                                                                                                                      >  2:1208879/1‑149 (MQ=255)
                 ccTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                                                                      <  1:475613/150‑1 (MQ=255)
                 ccTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                                                                      >  1:348133/1‑150 (MQ=255)
                 ccTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtc                                                                                                                                      <  1:249248/150‑1 (MQ=255)
                  cTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                                                                     <  1:1257589/150‑1 (MQ=255)
                  cTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtct                                                                                                                                     <  2:806193/150‑1 (MQ=255)
                   tCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                                                                    >  1:231874/1‑150 (MQ=255)
                   tCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTtctc                                                                                                                                    >  1:585584/1‑150 (MQ=255)
                    cAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTctct                                                                                                                                   <  2:263336/150‑1 (MQ=255)
                      aaGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCTGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAt                                                                                                                                 >  2:1154801/1‑150 (MQ=255)
                          aaGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGCCCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                                                                             <  1:1087913/150‑1 (MQ=255)
                          aaGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAAt                                                                                                                             <  1:1209976/150‑1 (MQ=255)
                          aaGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTAGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGAATTATATAAAAGGAACCAAAATTTATTGGGTTCTTTTTTATCTATCAAt                                                                                                                             >  2:1103302/1‑150 (MQ=255)
                              taTTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTa                                                                                                                         <  2:642666/149‑1 (MQ=255)
                                 tAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTa                                                                                                                      <  1:1299430/150‑1 (MQ=255)
                                  aaaaTACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAg                                                                                                                     >  1:940538/1‑150 (MQ=255)
                                     aTACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAg                                                                                                                  >  1:38954/1‑150 (MQ=255)
                                      tACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAgt                                                                                                                 <  2:627470/150‑2 (MQ=255)
                                       aCTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                                                                                >  2:117534/1‑150 (MQ=255)
                                       aCTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAggg                                                                                                                <  2:1004593/150‑1 (MQ=255)
                                        cTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGga                                                                                                               <  1:436235/150‑1 (MQ=255)
                                         ttACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGgag                                                                                                              <  1:762940/150‑1 (MQ=255)
                                           aCAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGat                                                                                                            <  1:821559/150‑1 (MQ=255)
                                              gACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatata                                                                                                         >  1:697073/1‑150 (MQ=255)
                                                cacTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGatatata                                                                                                       >  2:21146/1‑149 (MQ=255)
                                                cacTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                                                                       <  1:1339857/150‑1 (MQ=255)
                                                cacTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATc                                                                                                       <  2:180471/150‑1 (MQ=255)
                                                       aaaTTAATTTATGATTTGTGCGGGGTTAAATCGCAGAAATAAAATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTAACAATACTATTAGCAGGGAGATATATCACCagag                                                                                                >  2:1250135/1‑150 (MQ=255)
                                                        aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGt                                                                                               >  2:720053/1‑150 (MQ=255)
                                                        aattaattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGt                                                                                               >  2:1284801/1‑150 (MQ=255)
                                                            aattTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTaa                                                                                           <  1:1241254/150‑1 (MQ=255)
                                                             attTATGATTTGTTCTGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAt                                                                                          >  2:760420/1‑150 (MQ=255)
                                                                tATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtg                                                                                       <  1:1427948/150‑1 (MQ=255)
                                                                 aTGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgt                                                                                      >  1:121875/1‑150 (MQ=255)
                                                                  tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGTTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                                                                     <  2:348133/150‑1 (MQ=255)
                                                                  tGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAAtgtgtg                                                                                     <  2:585584/150‑1 (MQ=255)
                                                                        gTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAttttt                                                                               <  1:928168/150‑1 (MQ=255)
                                                                         ttCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGAtttttt                                                                              <  2:1091080/150‑1 (MQ=255)
                                                                          tCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTttta                                                                             <  1:863915/150‑1 (MQ=255)
                                                                             gggTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTtttattt                                                                          >  2:428188/1‑150 (MQ=255)
                                                                                 tAAATCGGTTAAATAATATTTTAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                                                                      <  2:869959/150‑1 (MQ=255)
                                                                                  aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                                                                     <  2:889469/150‑1 (MQ=255)
                                                                                  aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                                                                     <  2:1186959/150‑1 (MQ=255)
                                                                                  aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                                                                     <  1:1137141/150‑1 (MQ=255)
                                                                                  aaaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgt                                                                     >  1:641587/1‑150 (MQ=255)
                                                                                   aaTCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcgtc                                                                    <  2:820669/150‑1 (MQ=255)
                                                                                     tCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGa                                                                  <  2:481420/150‑1 (MQ=255)
                                                                                      cGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGaa                                                                 >  1:600986/1‑150 (MQ=255)
                                                                                            aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa                                                           >  1:40301/1‑150 (MQ=255)
                                                                                            aataatATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGa                                                           >  1:984736/1‑150 (MQ=255)
                                                                                              taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGAtt                                                         <  1:809071/150‑1 (MQ=255)
                                                                                              taataTGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACACATTATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATAACCAGAGTTTAATGTGTGATTTTATATTTATCGTCGAACCTGGAtt                                                         >  1:1383147/1‑150 (MQ=255)
                                                                                                atatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGt                                                       >  2:133742/1‑150 (MQ=255)
                                                                                                 tatGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGtt                                                      >  2:954432/1‑150 (MQ=255)
                                                                                                    gtgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTATTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTAt                                                   >  2:227521/1‑150 (MQ=255)
                                                                                                     tgAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATc                                                  >  2:126954/1‑150 (MQ=255)
                                                                                                           tcttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGc                                            >  1:838884/1‑150 (MQ=255)
                                                                                                            cttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTTGcc                                           >  1:727305/1‑150 (MQ=255)
                                                                                                            cttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGcc                                           <  2:144704/150‑1 (MQ=255)
                                                                                                            cttcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGca                                           >  1:1042873/1‑149 (MQ=255)
                                                                                                              tcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCtt                                         >  2:544095/1‑150 (MQ=255)
                                                                                                              tcGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCtt                                         >  2:394887/1‑150 (MQ=255)
                                                                                                                gATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTaa                                       >  2:405614/1‑150 (MQ=255)
                                                                                                                  tAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACa                                     >  2:775221/1‑150 (MQ=255)
                                                                                                                    gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                                                                      >  1:792737/1‑115 (MQ=255)
                                                                                                                    gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTAtcg                                                                      <  2:792737/115‑1 (MQ=255)
                                                                                                                    gTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACaaa                                   >  1:904652/1‑150 (MQ=255)
                                                                                                                     tATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAg                                  >  2:1059405/1‑150 (MQ=255)
                                                                                                                      aTGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGt                                 <  2:165190/150‑1 (MQ=255)
                                                                                                                       tGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTATATCAATACTATTAGCAGGGGGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATAATTGGCCTTAACAAAGtt                                >  2:613102/1‑150 (MQ=255)
                                                                                                                        ggAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTa                               >  1:196072/1‑150 (MQ=255)
                                                                                                                              cATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCt                         >  1:1001876/1‑150 (MQ=255)
                                                                                                                                ttATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCtaa                       >  1:969500/1‑150 (MQ=255)
                                                                                                                                   tataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCACGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCtaataa                    >  1:539929/1‑150 (MQ=255)
                                                                                                                                    ataAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAg                   <  2:442961/150‑1 (MQ=255)
                                                                                                                                     taAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGa                  >  1:560245/1‑150 (MQ=255)
                                                                                                                                     taAAAGGACCAAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTTGCCTTAACAAAGTTAACGGCTAATAAGa                  >  2:1365192/1‑150 (MQ=255)
                                                                                                                                      aaaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGat                 <  1:1186245/150‑1 (MQ=255)
                                                                                                                                       aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGatt                >  1:515151/1‑150 (MQ=255)
                                                                                                                                       aaaGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGatt                <  1:33436/150‑1 (MQ=255)
                                                                                                                                         aGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGattat              >  1:1365057/1‑150 (MQ=255)
                                                                                                                                          ggACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGattatt             >  1:1106981/1‑150 (MQ=255)
                                                                                                                                           gACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTAttt            <  1:1172360/150‑1 (MQ=255)
                                                                                                                                           gACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTCTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTAttt            <  1:667/150‑1 (MQ=255)
                                                                                                                                            aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTc           <  2:836553/150‑1 (MQ=255)
                                                                                                                                             cccAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTcc          <  1:6314/150‑1 (MQ=255)
                                                                                                                                               cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCAt        >  2:377747/1‑150 (MQ=255)
                                                                                                                                               cAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCAt        <  2:1017104/150‑1 (MQ=255)
                                                                                                                                                  tatttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCaa     >  1:1269881/1‑149 (MQ=255)
                                                                                                                                                   atttattGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATAATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACt    >  2:1316032/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTACAGACACTTCTAAATTAATTTATGATTTGTTCGGGGTTAAATCGGAGAAATAATATGTGAAGTTCTTCGATAGTATGGAAGGCATTATATAAAAGGACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTC  >  NC_000913/4261210‑4261508

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