Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,414,243 A→T 5.4% intergenic (+52/‑32) yjfC → / → aidB ATP‑Grasp family ATPase/DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,414,2430AT 5.4% 397.5 / 9.1 114intergenic (+52/‑32)yjfC/aidBATP‑Grasp family ATPase/DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (53/53);  new base T (2/4);  total (57/57)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 116 total aligned reads displayed.

GTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGA  >  NC_000913/4414096‑4414390
                                                                                                                                                   |                                                                                                                                                   
gTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTaa                                                                                                                                                    <  3:264819/149‑1 (MQ=255)
 tCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAAt                                                                                                                                                   <  1:12816/149‑1 (MQ=255)
  cGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATg                                                                                                                                                  <  1:194813/149‑1 (MQ=255)
   gatgatGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATgg                                                                                                                                                 >  4:428187/1‑149 (MQ=255)
   gatgatGAAGCGTGCGGAATGGGATCCGTGGAGATAAAACCCTTATCACCAAAGACACCTCACGTTTCGTTCCGCCTTACATTGCTTGATAAGAATGGTTTAGCAATCTTCTTCTGTCATGAATCCATTGCAGTAACCATACTAATTGt                                                                                                                                                >  1:107218/1‑149 (MQ=255)
    atgatgAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGt                                                                                                                                                >  3:120101/1‑149 (MQ=255)
    atgatgAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGt                                                                                                                                                >  2:629756/1‑149 (MQ=255)
           aGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGcc                                                                                                                                         <  2:271233/149‑1 (MQ=255)
             cGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCAt                                                                                                                                       <  2:230259/149‑1 (MQ=255)
             cGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCAt                                                                                                                                       <  4:535403/149‑1 (MQ=255)
              gTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCAtt                                                                                                                                      >  4:616922/1‑149 (MQ=255)
                  ggAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAAt                                                                                                                                                   <  3:376281/132‑1 (MQ=255)
                  ggAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCAATGGCAGTGACCATACTAATGGTGACTGCCATTGATg                                                                                                                                  >  2:364263/1‑149 (MQ=255)
                   gAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGTCCATACTAATGGTGACTGCCATTGATgg                                                                                                                                 >  4:656698/1‑149 (MQ=255)
                       ggggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggaggg                                                                                                                             >  4:410774/1‑149 (MQ=255)
                         ggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATgga                                                                                                                                <  4:575864/144‑1 (MQ=255)
                         ggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATgga                                                                                                                                >  3:575864/1‑144 (MQ=255)
                          gATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGgaga                                                                                                                          <  4:30052/149‑1 (MQ=255)
                           aTCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGac                                                                                                                         <  2:107218/149‑1 (MQ=255)
                              cGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGAcaca                                                                                                                      <  3:210617/149‑1 (MQ=255)
                               gTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAg                                                                                                                     >  1:394188/1‑149 (MQ=255)
                                tGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGt                                                                                                                    <  3:410774/149‑1 (MQ=255)
                                       aaCACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTgg                                                                                                             >  4:458823/1‑149 (MQ=255)
                                        acacacTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGc                                                                                                            >  4:307672/1‑149 (MQ=255)
                                          acacTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGCTAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCaa                                                                                                          >  2:204383/1‑149 (MQ=255)
                                           cacTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCaaa                                                                                                         >  3:204225/1‑149 (MQ=255)
                                           cacTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCaaa                                                                                                         <  1:464842/149‑1 (MQ=255)
                                           cacTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCaaa                                                                                                         <  1:641345/149‑1 (MQ=255)
                                              tgctCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGAACATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTc                                                                                                      >  2:382265/4‑149 (MQ=255)
                                              tGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTc                                                                                                      >  3:666810/1‑149 (MQ=255)
                                              tGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTc                                                                                                      >  1:606668/1‑149 (MQ=255)
                                                aTCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTcac                                                                                                    >  2:479032/1‑149 (MQ=255)
                                                aTCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTcac                                                                                                    >  1:457932/1‑149 (MQ=255)
                                                aTCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCATGGAttcatga                                                                                                                             <  1:82903/124‑7 (MQ=25)
                                                aTCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCATGGAttcatga                                                                                                                             >  2:82903/1‑118 (MQ=25)
                                                 tCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTcaca                                                                                                   >  4:286344/1‑149 (MQ=255)
                                                 tCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTcaca                                                                                                   <  1:553844/149‑1 (MQ=255)
                                                  cACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCacac                                                                                                  >  2:361434/1‑149 (MQ=255)
                                                    ccAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTCATGTTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCg                                                                                                >  3:385882/1‑149 (MQ=255)
                                                    ccAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCg                                                                                                >  4:307381/1‑149 (MQ=255)
                                                       aaGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGttt                                                                                             <  3:428187/149‑1 (MQ=255)
                                                          acacCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTa                                                                                          <  1:88394/149‑1 (MQ=255)
                                                            acCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCacac                                                                                                  >  3:681764/1‑139 (MQ=255)
                                                            acCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAAt                                                                                        <  4:204225/149‑1 (MQ=255)
                                                              cTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGGGGGAGACACAg                                                                                                                     <  3:574533/118‑1 (MQ=255)
                                                              cTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGGGGGAGACACAg                                                                                                                     >  4:574533/1‑118 (MQ=255)
                                                               tCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTGATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCaa                                                                                     >  3:501133/1‑149 (MQ=255)
                                                               tCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggag                                                                                                                               >  4:628407/1‑107 (MQ=255)
                                                               tCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggag                                                                                                                               <  3:628407/107‑1 (MQ=255)
                                                                cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAc                                                                                    <  3:542864/149‑1 (MQ=255)
                                                                cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAc                                                                                    <  4:584312/149‑1 (MQ=255)
                                                                  cGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCt                                                                                  <  2:457932/149‑1 (MQ=255)
                                                                         ttCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCAt                                                                           >  2:652772/1‑149 (MQ=255)
                                                                             gCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTaaa                                                                       <  2:3950/149‑1 (MQ=255)
                                                                                ttACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATaa                                                                    <  2:659299/149‑1 (MQ=255)
                                                                                ttACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATaa                                                                    <  4:385882/149‑1 (MQ=255)
                                                                                  aCATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACa                                                                  <  1:652772/149‑1 (MQ=255)
                                                                                  aCATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACa                                                                  <  1:479032/149‑1 (MQ=255)
                                                                                       gCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAAc                                                             <  1:594743/149‑1 (MQ=255)
                                                                                       gCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAAc                                                             <  3:307381/149‑1 (MQ=255)
                                                                                          ggATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTta                                                          <  1:629756/149‑1 (MQ=255)
                                                                                           gATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCtata                                                                               <  3:668477/127‑1 (MQ=255)
                                                                                           gATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCtata                                                                               >  4:668477/1‑127 (MQ=255)
                                                                                             tAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATAc                                                       <  4:666810/149‑1 (MQ=255)
                                                                                                gAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTg                                                    >  2:579119/1‑149 (MQ=255)
                                                                                                gAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTg                                                    >  4:565900/1‑149 (MQ=255)
                                                                                                 aaTGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGt                                                   >  4:317284/1‑149 (MQ=255)
                                                                                                 aaTGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGt                                                   >  1:486874/1‑149 (MQ=255)
                                                                                                 aaTGTTTTAGCAATATCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGAGGGGGCGAGACACAGTGCACTGGAAAATTCACAACGTTTTTAACCAACCTACACTATTAAATAACAACAACTTATACCTGt                                                   >  1:93249/1‑149 (MQ=255)
                                                                                                  aTGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCaa                                                              <  3:129650/137‑1 (MQ=255)
                                                                                                  aTGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCaa                                                              >  4:129650/1‑137 (MQ=255)
                                                                                                   tGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCt                                                 >  4:656564/1‑149 (MQ=255)
                                                                                                     ttttAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGa                                               >  3:140854/1‑149 (MQ=255)
                                                                                                     ttttAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGa                                               >  4:137837/1‑149 (MQ=255)
                                                                                                            aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTg                                                                                                                   <  2:185109/74‑1 (MQ=255)
                                                                                                            aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTg                                                                                                                   >  1:185109/1‑74 (MQ=255)
                                                                                                            aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                              <  1:170043/111‑1 (MQ=255)
                                                                                                            aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                              >  2:170043/1‑111 (MQ=255)
                                                                                                                tctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCACGctct                                    <  3:260684/149‑1 (MQ=255)
                                                                                                                  tttCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGGCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                              <  4:369306/105‑1 (MQ=255)
                                                                                                                  tttCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGGCACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                              >  3:369306/1‑105 (MQ=255)
                                                                                                                      tGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAg                              >  1:5988/1‑149 (MQ=255)
                                                                                                                        tCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                            <  4:426248/149‑1 (MQ=255)
                                                                                                                         cATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCgg                           <  1:501887/149‑1 (MQ=255)
                                                                                                                           tccaTCAATGGCAGTCACCATTAGTATGGTCACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgc                                           >  2:399491/4‑131 (MQ=37)
                                                                                                                           tccaTCAATGGCAGTCACCATTAGTATGGTCACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgc                                           <  1:399491/128‑1 (MQ=37)
                                                                                                                                 cATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGt                   >  4:442882/1‑149 (MQ=255)
                                                                                                                                 cATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGt                   <  1:204383/149‑1 (MQ=255)
                                                                                                                                    ggCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGaa                <  3:656564/149‑1 (MQ=255)
                                                                                                                                     gCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACtt                                                           <  2:358620/105‑1 (MQ=255)
                                                                                                                                     gCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACtt                                                           >  1:358620/1‑105 (MQ=255)
                                                                                                                                      cAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAgg              >  4:573129/1‑149 (MQ=255)
                                                                                                                                      cAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAgg              >  2:497662/1‑149 (MQ=255)
                                                                                                                                       aGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGt             >  2:687375/1‑149 (MQ=255)
                                                                                                                                       aGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGt             >  1:152241/1‑149 (MQ=255)
                                                                                                                                        gTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTg            <  2:500387/149‑1 (MQ=255)
                                                                                                                                          gACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgct          <  1:623110/149‑1 (MQ=255)
                                                                                                                                            ccATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctgg        >  4:350158/1‑149 (MQ=255)
                                                                                                                                               tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg     <  1:579119/149‑1 (MQ=255)
                                                                                                                                                  tAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGa  >  1:353239/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGA  >  NC_000913/4414096‑4414390

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