Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,998,056 T→G 6.0% L25R (CTT→CGT)  uvrD → DNA‑dependent ATPase I and helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,998,0560TG 6.0% 1368.0 / 3.6 486L25R (CTT→CGT) uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base T (165/249);  new base G (40/2);  total (234/252)
Fisher's exact test for biased strand distribution p-value = 4.13e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 493 total aligned reads displayed.

ATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCG  >  NC_000913/3997907‑3998205
                                                                                                                                                     |                                                                                                                                                     
aTCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCt                                                                                                                                                       >  2:2052473/1‑150 (MQ=255)
 tCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCtt                                                                                                                                                      >  1:7139869/1‑150 (MQ=255)
  cTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATTACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTc                                                                                                                                                     >  1:933050/1‑150 (MQ=255)
    gTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTg                                                                                                                                                   <  2:2450673/150‑1 (MQ=255)
    gTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTg                                                                                                                                                   <  2:6586885/150‑1 (MQ=255)
    gTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTg                                                                                                                                                   <  2:5254119/150‑1 (MQ=255)
    gTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTg                                                                                                                                                   <  1:2549179/150‑1 (MQ=255)
    gTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTg                                                                                                                                                   <  1:851614/150‑1 (MQ=255)
     tatatataCCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTgg                                                                                                                                                  >  1:4052466/1‑150 (MQ=255)
       tatataCCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTg                                                                                                                                                >  1:361579/1‑150 (MQ=255)
       tatataCCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCAGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCACCAACCGTCTGGTg                                                                                                                                                >  1:4389223/1‑150 (MQ=255)
         tataCCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCt                                                                                                                                              <  1:6137598/150‑1 (MQ=255)
           taCCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTgg                                                                                                                                            <  1:1839101/150‑1 (MQ=255)
              ccAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTggcgg                                                                                                                                         >  1:5256647/1‑150 (MQ=255)
              ccAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTggcgg                                                                                                                                         >  2:5819427/1‑150 (MQ=255)
               cAGCTTTTTGGCGGAGGTCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGAGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAAAATCGCGAAGCGGTGGCCGCGCCAAGCAGCAACCATCTCGTGCTggcggg                                                                                                                                        >  2:6326445/1‑150 (MQ=255)
               cAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTGTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGGCGCGCCACGCAGCAACGTTCTGGTGCTggcggg                                                                                                                                        >  1:4230342/1‑150 (MQ=255)
               cAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTggcggg                                                                                                                                        >  1:1217059/1‑150 (MQ=255)
                   tttttGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGgcgggcgcg                                                                                                                                    >  1:3462244/1‑150 (MQ=255)
                   tttttGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCAGCTGGTGCAGgcgggcgcg                                                                                                                                    >  1:6057881/1‑150 (MQ=255)
                     tttGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcggg                                                                                                                                  <  1:2930168/150‑1 (MQ=255)
                     tttGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcggg                                                                                                                                  <  2:3342548/150‑1 (MQ=255)
                       tGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATAGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGgcgggcgcgggcc                                                                                                                                >  1:5754118/1‑149 (MQ=255)
                         gCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                                                                                                              <  1:2797859/150‑1 (MQ=255)
                         gCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                                                                                                              <  1:1014758/150‑1 (MQ=255)
                         gCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                                                                                                              <  2:4386455/150‑1 (MQ=255)
                          cGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                                                                                                             >  2:134994/1‑150 (MQ=255)
                            gAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGt                                                                                                                           <  2:3546324/150‑1 (MQ=255)
                             tggggGTTGCGCGTCTCCGCCCAACCTATTATTACGCGGCGGTGCGACTGGACGTTTCTTACCTGCTCGCCCGCCTGCCTGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGTGCGGGCGGGAGCGGGGATTGGTa                                                                                                                          <  2:1267653/149‑1 (MQ=255)
                             aGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTa                                                                                                                          <  1:6452752/150‑1 (MQ=255)
                             aGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTa                                                                                                                          <  1:667023/150‑1 (MQ=255)
                              gggCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTaa                                                                                                                         <  1:7005718/150‑1 (MQ=255)
                                gCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGACGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGa                                                                                                                       <  1:231732/150‑1 (MQ=255)
                                gCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGa                                                                                                                       <  2:3893661/150‑1 (MQ=255)
                                gCGTTGCGCATCTACGCCAAACCTAATTATACGCGGCGGTGCCAATGGACGTTTCTTAAATGCTCGACAGCCTAAAAGACAAACAGCTCGAAGCGGTGTCCCCGCCACGCAGCAACCGTCTGGTGCTGGCGGGCCCGGGCAggggacaga                                                                                                                       >  2:3822152/1‑145 (MQ=255)
                                 cGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAc                                                                                                                      <  2:4311637/150‑1 (MQ=255)
                                    tGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCCACCTTCTGGTGCTGGCGGGCGCGGGCCGTGGTAAGAcgcg                                                                                                                   >  1:5679511/1‑150 (MQ=255)
                                      cgcTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACTTTTCTTACCTGCTCGACAGCCTTAATGAAAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgcg                                                                                                                 <  2:2592453/150‑1 (MQ=255)
                                      cgcTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCGGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgcg                                                                                                                 >  1:688702/1‑150 (MQ=255)
                                       gcTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGt                                                                                                                >  2:1695578/1‑150 (MQ=255)
                                           ctcCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTg                                                                                                            <  1:4542703/150‑1 (MQ=255)
                                           ctcCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTg                                                                                                            <  2:6677318/150‑1 (MQ=255)
                                           ctcCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTg                                                                                                            <  2:1578057/150‑1 (MQ=255)
                                            tcCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTgg                                                                                                           >  2:6361577/1‑150 (MQ=255)
                                             ccGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGt                                                                                                          >  2:4556351/1‑150 (MQ=255)
                                             ccGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGt                                                                                                          <  1:3834431/150‑1 (MQ=255)
                                               gCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCATCTGGTGCTGGCGGGCGCGGGCCGTGGTAAGACGCGCGTACAGGTGc                                                                                                        >  1:5011111/1‑150 (MQ=255)
                                                  cAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATc                                                                                                     <  2:2287694/150‑1 (MQ=255)
                                                   aaCCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGGGGTAAGACGCGCGTACTGGTGCAGCg                                                                                                    >  2:6292787/1‑150 (MQ=255)
                                                      cTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTAt                                                                                                 >  2:4759534/1‑150 (MQ=255)
                                                      cTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCATATGGTGCTGGCGGGCGCGGGCCGTGGTAAGACGCGCGTAAGGGTGCAGCGTAt                                                                                                 >  1:6182792/1‑150 (MQ=255)
                                                       tATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCCCCGCAACCCACTGGTGGTGGGGGGCGCGGGCCGGGGTTAGACGCGCGGACAGGTGCaaaggata                                                                                                >  1:3132399/1‑143 (MQ=255)
                                                         tttttACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCATCAGGTGCTGGCGGGCGCGGGCAGGGGTAAGACGCGCGTACTGGTGCCTCGTATCGc                                                                                              >  1:5980745/1‑150 (MQ=255)
                                                             tACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCAGCGTATCGCCTgg                                                                                          >  2:5012490/1‑150 (MQ=255)
                                                              aCGCGGCGGGGCCATTGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGGGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCTGTGGTATGACGCGCGTACTGGTGCATCGTTTCGCCTGGt                                                                                         <  2:1498300/150‑1 (MQ=255)
                                                               cgcgGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACGGGTGCATCGTATCGCCTGGtt                                                                                        >  2:3121411/1‑150 (MQ=255)
                                                               cgcgGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGtt                                                                                        <  1:777755/150‑1 (MQ=255)
                                                                        gCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCTGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGt                                                                               <  1:1997984/150‑1 (MQ=255)
                                                                         ccAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTg                                                                              <  1:2681077/150‑1 (MQ=255)
                                                                          cAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTgg                                                                             >  2:5368547/1‑150 (MQ=255)
                                                                           aaTGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGa                                                                            >  2:4950224/1‑150 (MQ=255)
                                                                              ggACGTTTCTGACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGaaaa                                                                         <  2:5548682/150‑1 (MQ=255)
                                                                                  gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCGTCGGGTGCTGGCGGGCGCGGGCCGTGGTAAGACGCGCGTACTGGTGCCTCGTATCGCCGGGTTGATGGTCGTGGAAAACGGc                                                                     >  1:2375345/1‑150 (MQ=255)
                                                                                  gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCATCAGGTGCTGGCGGGCGCGGGCAGTGGTAAGACCCGCGTACTGGTGCATCGTATAGCCTGGTTGATGAGCGTGGAAAACttc                                                                     >  1:4104928/1‑148 (MQ=255)
                                                                                  gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAACGTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGc                                                                     >  2:3161594/1‑150 (MQ=255)
                                                                                     tCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCg                                                                  <  2:6857333/150‑1 (MQ=255)
                                                                                      cTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGc                                                                 <  2:6176460/150‑1 (MQ=255)
                                                                                       ttACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACGGCTCGcc                                                                >  2:5622237/1‑150 (MQ=255)
                                                                                          ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATa                                                             >  2:5137323/1‑150 (MQ=255)
                                                                                           cTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAGGACGCGCGTACTGGTGCATCGGATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAc                                                            <  1:1281267/150‑1 (MQ=255)
                                                                                             gCTCGACAGCCTTAATGACACACCGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCGGGGGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTc                                                          <  2:362107/150‑1 (MQ=255)
                                                                                             gCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCTTCGGGTGCTGGCGGGCGCGGGCAGGGGTAAGACGCGCGTACAGGTGCAACGTATCGCCGGGTTGATGAGCGTGGAAAACTGCTCGCCATAcac                                                          >  2:1987610/1‑148 (MQ=255)
                                                                                               tCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGa                                                        <  2:2766843/150‑1 (MQ=255)
                                                                                                cGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGAt                                                       <  2:4861545/150‑1 (MQ=255)
                                                                                                      ccTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCGGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGc                                                 >  2:1442178/1‑150 (MQ=255)
                                                                                                      ccTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGc                                                 <  1:5165870/150‑1 (MQ=255)
                                                                                                      ccTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACACGATTATGGc                                                 >  2:2020606/1‑150 (MQ=255)
                                                                                                        ttAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCgg                                               >  2:7205574/1‑150 (MQ=255)
                                                                                                        ttAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTATGGCga                                               <  2:602243/150‑2 (MQ=255)
                                                                                                          aaTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTg                                             >  2:6070739/1‑150 (MQ=255)
                                                                                                                  aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTAc                                     <  1:161260/150‑1 (MQ=255)
                                                                                                                    aGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACca                                   >  2:2851220/1‑150 (MQ=255)
                                                                                                                     gcgcGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACcaa                                  <  1:1486249/150‑1 (MQ=255)
                                                                                                                        cgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACcaacaa                               <  2:4079680/150‑1 (MQ=255)
                                                                                                                         gAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACaaa                              <  1:74663/150‑1 (MQ=255)
                                                                                                                          aaGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCCGTGGTAAGACGCGCGTACAGGTGCAACGTATCGCCTGGTTGATGAGCGTGGAAAAATGCTCGCCATAATCGATTATGGCGGTGACGTTTACCAACAAAg                             >  2:2225955/1‑150 (MQ=255)
                                                                                                                          aaGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCGGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGTCGGTTACCAACAAAg                             >  2:2909415/1‑150 (MQ=255)
                                                                                                                             cGGTGGCCGCCCCACGCAGCACCCCGCCGGCGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGGGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcgg                          <  1:5051834/150‑1 (MQ=255)
                                                                                                                               gTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggcg                        <  1:2171250/150‑1 (MQ=255)
                                                                                                                                   ccGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGcggcggcgg                    <  2:1181643/150‑1 (MQ=255)
                                                                                                                                     gcgcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGgag                  <  2:4442448/150‑1 (MQ=255)
                                                                                                                                     gcgcCACGCAGCAACCTGCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGgag                  >  2:800879/1‑150 (MQ=255)
                                                                                                                                      cgcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGgaga                 >  2:3037441/1‑150 (MQ=255)
                                                                                                                                        ccACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATg               <  2:7010131/150‑1 (MQ=255)
                                                                                                                                         cACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGc              >  1:4656099/1‑150 (MQ=255)
                                                                                                                                         cACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGc              <  1:630128/150‑1 (MQ=255)
                                                                                                                                             cagcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCa          >  1:3953684/1‑150 (MQ=255)
                                                                                                                                               gcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATc        <  1:5438175/150‑1 (MQ=255)
                                                                                                                                               gcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATc        <  1:7212252/150‑1 (MQ=255)
                                                                                                                                                  aCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTa     >  2:2188153/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCTGTATATATACCCAGCTTTTTGGCGGAGGGCGTTGCGCTTCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCG  >  NC_000913/3997907‑3998205

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑