Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,254 A→T 6.5% intergenic (+32/‑115) lysY → / → lysZ tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,2540AT 6.5% 1488.1 / 50.1 459intergenic (+32/‑115)lysY/lysZtRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base A (225/199);  new base T (15/19);  total (240/219)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 481 total aligned reads displayed.

AGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTT  >  NC_000913/781105‑781403
                                                                                                                                                     |                                                                                                                                                     
 ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                                      >  1:786324/1‑150 (MQ=255)
 ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                                      <  2:2775479/150‑1 (MQ=255)
 ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGaa                                                                                                                                                      <  2:3007400/150‑1 (MQ=255)
    gggggTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAAAACAGATGTTAAGCGTGAAGGa                                                                                                                                                   >  1:3167873/1‑150 (MQ=255)
      ggctCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCTTTAGCTCAGCTGGTAGAGCAGTTGACTTATAATCTATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATa                                                                                                                                                 <  2:6900187/147‑1 (MQ=14)
         tCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGGGAAGGATAAct                                                                                                                                              >  2:3220054/1‑149 (MQ=255)
             cggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgc                                                                                                                                          <  2:8670705/150‑1 (MQ=255)
             cggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCAGCACGACCCACCAGTTTTAAAATCAAACTCAGATGTTAATCGTTAAGGATAAACTTgc                                                                                                                                          >  1:6973989/1‑150 (MQ=14)
              ggTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTgcg                                                                                                                                         >  1:9732247/1‑150 (MQ=255)
               gTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTAAACTTGGTACTAGGATGTTAAGCGTGAGGGTTAACGtagagc                                                                                                                                        >  2:4908773/1‑145 (MQ=14)
                ttCGATCCCGTCATCACCCACCACTCGGGTCGTTATCTCAGTTGGTAGAGCAGTTGACTTTTAAAAAATTGGTCGCAGGTTAGAATCCAGCACGACCCACCAGGTTTAACATAAAAATCAGAAGTTAAGCGTGAAGGATAACGGTGAGcc                                                                                                                                       >  1:2192090/1‑150 (MQ=14)
                   gATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGATGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATAATGCAAGACCCACCAGTTTTAACATAAAAATAAGATGTTAAGAGTGAAGGATAAAGTTGAGACAGc                                                                                                                                    >  1:5442546/1‑150 (MQ=14)
                    aTCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTTTTAACATCAAACACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCa                                                                                                                                   >  1:667134/1‑150 (MQ=255)
                         gTCATCACCCACCACTCGGGTCGTTAGGTAAGTTGGTATAGCAGTTTACTTTTTATAAATTGGTGGGAGGTTCGAATACTGCACGACCAAACAGTTTTAACATCAAAATAATATTTTTTGGGTGTAGGATAAAGTTTGGGCAGAAACGGc                                                                                                                              >  1:9491512/1‑150 (MQ=255)
                          tcatcaCCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACgtta                                                                                                                             <  1:8632654/150‑4 (MQ=255)
                                     cacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTCAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGGAGGgcgag                                                                                                                  >  2:6443412/1‑150 (MQ=255)
                                     cacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgcgag                                                                                                                  <  2:7966056/150‑1 (MQ=255)
                                     cacTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgcgag                                                                                                                  >  2:1379309/1‑150 (MQ=255)
                                      acTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGGTAACGTTGCGCCAGCAACGGCCCGTAGGgcgagc                                                                                                                 <  2:9699499/150‑1 (MQ=255)
                                      acTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGgcgagc                                                                                                                 <  2:5176045/150‑1 (MQ=255)
                                          gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaag                                                                                                             <  2:8778707/150‑1 (MQ=255)
                                          gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaag                                                                                                             >  2:5771714/1‑150 (MQ=255)
                                             tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                                          <  2:7437362/150‑1 (MQ=255)
                                             tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                                          <  2:3554350/150‑1 (MQ=255)
                                             tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                                          >  1:5768284/1‑150 (MQ=255)
                                             tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCAGGTTCGAATCCAGCACGACCCACCATATTAAAAATAAAACAAAGATGATAACCGTGAAGGATAACGTTGCGAAAGCAACGGCCCGTAGGGCGagcgaagcga                                                                                                          >  1:3889222/1‑150 (MQ=11)
                                               gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGt                                                                                                        >  2:8971630/1‑150 (MQ=255)
                                                ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGa                                                                                                                                                       >  2:9603490/1‑102 (MQ=40)
                                                 tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCa                                                                                                      <  2:333755/150‑1 (MQ=255)
                                                    cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATcc                                                                                                   >  1:9859296/1‑150 (MQ=255)
                                                     tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCt                                                                                                  >  2:7268022/1‑150 (MQ=255)
                                                        gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCa                                                                                               >  1:1478558/1‑150 (MQ=255)
                                                         ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACAAAGATGTTAAGCGTGAAGCATAACGTTGCGACAGCAACGGCCCGAAGGGCGAGCGAAGCGAGTCATACTGCAc                                                                                              >  1:4625877/1‑150 (MQ=11)
                                                             tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                                                                                          >  2:7737537/1‑150 (MQ=255)
                                                                     atGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGTGAGTCATCCTGCACGACCCACCAATc                                                                                  >  2:2501593/2‑150 (MQ=255)
                                                                          ttttAATCAATTGGTCGCAGGTTCTAATCCTGCACGACCCACCATTTTTAACATCAAACTCAGATGTTAAGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                             <  1:1614721/150‑1 (MQ=255)
                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCTCGACCCACCAGTTTTAACATCAAACTATGATGTTAAGCGTGTCGGATAACGTTGCGCCAGCAACGTCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                             <  1:9052929/150‑1 (MQ=14)
                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                             >  1:10043726/1‑150 (MQ=255)
                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATACTGCACGACCCAACAATCTTaaa                                                                             >  2:1727943/1‑150 (MQ=255)
                                                                           tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAg                                                                            >  1:2823913/1‑150 (MQ=255)
                                                                           tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTATAACATCAAACTCAGATGTTAAACGTGAAGGATAAAGTTCAGCCAGCAACGGCACGAAGTGCGAGCGAAGCGAGACATCATGCACGACCCACCAATCAAAAAg                                                                            >  1:9578721/1‑150 (MQ=14)
                                                                               aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                                        <  2:7075361/150‑1 (MQ=255)
                                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTgg                                                                       <  1:6504361/150‑1 (MQ=255)
                                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGATAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCAGCACGACCCACCAAACTAAAATATTgg                                                                       >  2:4031796/1‑150 (MQ=14)
                                                                                tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACGCAGATGATAACCGGGAAGGATAACGTTCCGACAGCAACGGCCCGTAGGCCGAGCGAAGCGAGTCATCCTGCCCCACCCACCAATCTTAAAGGTTga                                                                       >  2:4749211/1‑149 (MQ=14)
                                                                                  aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcc                                                                     >  1:3119051/1‑150 (MQ=255)
                                                                                  aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACACAGAAGTTAACCGTGAAGGATAACGTTGCTCAAGCACCGGCCCGAAGGGCGGGCGAAGAGAGTCATACAGCACCACCCCCCAATCATAAAGGTTGGcc                                                                     >  1:3628278/1‑150 (MQ=14)
                                                                                      ggTCGCTGGTTCGAATCCTGCACGACCCACCATTTTTAACATCAAACTCAGATTTTATGCGTGAAGTATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                 <  1:9350184/150‑1 (MQ=18)
                                                                                      ggTCGCAGGCACGAATCCCGCACGCCCCACCAGCTTTAACTTAATTATCAGTTGTTAAGCGTGATGGATAACGTTGCGCCACCTACGGCCCGTAGGGAGCGCGTTGCGATAGATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                 <  2:396420/150‑1 (MQ=255)
                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                               >  1:1266636/1‑150 (MQ=255)
                                                                                        tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTAATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                               >  2:1218406/1‑150 (MQ=255)
                                                                                          gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGtca                                                             >  1:479443/1‑148 (MQ=255)
                                                                                           cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaa                                                            <  2:8053910/150‑1 (MQ=255)
                                                                                            aGGTTCGACTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaa                                                           <  2:4297606/150‑1 (MQ=255)
                                                                                            aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTATGGTTGTTGGTTATGGTTGCGCCAGCTTCGGCCTGTAGGGCGTGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTaaaa                                                           <  1:9148776/150‑1 (MQ=18)
                                                                                              gTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTGAAGATTGGCCCCGAGTAAAAAt                                                         >  1:2692742/1‑150 (MQ=255)
                                                                                                tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCCAGTAAAAATCt                                                       >  1:8783086/1‑150 (MQ=255)
                                                                                                   aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                    <  2:447056/150‑1 (MQ=255)
                                                                                                    accccgCACGACCCACCAGATTCAACATCCAACACATATGTTAAGCGTGAAGGATATCTTTGAGGCAGCTACGGCCGTTATGGCGTGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                   <  1:1970544/145‑1 (MQ=18)
                                                                                                    aTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCa                                                   <  1:8971630/150‑1 (MQ=255)
                                                                                                    aTCCAGCACGACCCACCACTTTTAACATCAAACACAGATCTTAAGCGTCAACGATAACGTAGCGCCAGCAACGCCCCGAAGCCCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAACATTGGCCCCGAGTAAAAATCTTTCa                                                   <  2:6108190/150‑1 (MQ=18)
                                                                                                       cTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGACAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                                                >  1:5211660/1‑150 (MQ=255)
                                                                                                        tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCTGGTa                                               >  1:5873375/1‑150 (MQ=255)
                                                                                                         gCACGACCCACCATTTTTAACATCAAAATCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCAGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                              <  2:1962981/150‑1 (MQ=255)
                                                                                                             gACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAcc                                          >  1:5590751/1‑150 (MQ=255)
                                                                                                             gACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAcc                                          >  1:3380431/1‑150 (MQ=255)
                                                                                                             gACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACAcc                                          >  1:7740563/1‑150 (MQ=255)
                                                                                                                  accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                     >  2:8210370/1‑150 (MQ=255)
                                                                                                                  accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTAt                                     <  1:1111045/150‑1 (MQ=255)
                                                                                                                   ccaGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                                    <  2:1020925/150‑1 (MQ=17)
                                                                                                                   ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                                    <  2:9346051/150‑1 (MQ=255)
                                                                                                                   ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATg                                    <  2:53617/150‑1 (MQ=255)
                                                                                                                    caGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATgg                                   >  1:7603537/1‑150 (MQ=255)
                                                                                                                     atttttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATggg                                  <  1:1715225/148‑1 (MQ=255)
                                                                                                                      gTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGt                                 <  1:6493765/150‑1 (MQ=255)
                                                                                                                         ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt                              >  1:4055013/1‑150 (MQ=255)
                                                                                                                         ttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGt                              <  2:2312956/150‑1 (MQ=255)
                                                                                                                           aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa                            <  1:6016317/150‑1 (MQ=255)
                                                                                                                           aaCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTa                            <  2:7484374/150‑1 (MQ=255)
                                                                                                                             cATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAACTTAAAGATTGGCCCCGAGTAAAAATCATTCAGGTCACACCCGAATGGGTCGTTAga                          >  2:3884294/1‑149 (MQ=255)
                                                                                                                              aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATGGGCCCCGAGTAAAAAACTGTCAGGTAACACCCGTATGGGTCGTTAGcg                         >  1:6280688/1‑149 (MQ=255)
                                                                                                                              aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGAAGGGCGAGAGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAgat                         >  1:7190334/1‑148 (MQ=255)
                                                                                                                                cAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa                       <  2:6093707/150‑1 (MQ=255)
                                                                                                                                cAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCa                       <  1:5045469/150‑1 (MQ=255)
                                                                                                                                  aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGt                     >  1:8201107/1‑150 (MQ=255)
                                                                                                                                     tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGAGTAGCGAGTCAACAAGCACGACACACCACACAATAAGATTTGCCCCGAGAAGAAGACTGTCGGGTAGCACACGTGTGTGTCGTTAGCtgtggatg                  >  2:9462312/1‑143 (MQ=9)
                                                                                                                                       aGATGTTAAGCGTGTAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                                      <  1:4923515/112‑1 (MQ=255)
                                                                                                                                       aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTa                <  2:9246667/150‑1 (MQ=255)
                                                                                                                                       aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGGTggga                >  2:5531308/1‑148 (MQ=255)
                                                                                                                                        gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTag               >  2:2872382/1‑150 (MQ=255)
                                                                                                                                          tGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTagag             <  1:24283/150‑1 (MQ=255)
                                                                                                                                          agatAAGCGAGAAGGATAACCTAGCGCCACCAACGCCCCGTAGGGCGAGCGAAGCGAGACATCCTGCACGACCCACCAAACAAAAACATTGGCCCCGAGTAAAAAACTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTagag             <  2:9898418/147‑1 (MQ=25)
                                                                                                                                            ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAATATTGGCCCCCAGTAAAAATCTTTCTGGTACCACCCGTATGGGTCGTGAGCACGGTGGGTagagga           >  2:5260449/1‑148 (MQ=21)
                                                                                                                                               aGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGtt        >  1:4304834/1‑150 (MQ=255)
                                                                                                                                                gCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTGGGTAGAGCAGTTg       >  1:4716275/1‑150 (MQ=255)
                                                                                                                                                    gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGTGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCCTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt   <  2:7825876/150‑1 (MQ=255)
                                                                                                                                                    gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt   >  1:5297531/1‑150 (MQ=255)
                                                                                                                                                    gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt   >  2:46505/1‑150 (MQ=255)
                                                                                                                                                    gAAGGATAACGGTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTCCAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt   <  1:8094973/150‑1 (MQ=255)
                                                                                                                                                     aaGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt  <  2:1600710/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTT  >  NC_000913/781105‑781403

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