Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,811,334 A→G R28R (CGT→CGC)  rpmG ← 50S ribosomal subunit protein L33

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,811,3340AG98.8% 309.5 / NA 82R28R (CGT→CGC) rpmG50S ribosomal subunit protein L33
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (50/31);  minor base T (0/1);  total (50/32)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTG  >  NC_000913/3811187‑3811478
                                                                                                                                                   |                                                                                                                                                 
gggggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCg                                                                                                                                                 <  2:216595/150‑1 (MQ=255)
      aaGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGtt                                                                                                                                           <  1:341375/150‑1 (MQ=255)
             tAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGt                                                                                                                                     >  2:902097/1‑149 (MQ=255)
             tAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGt                                                                                                                                     >  2:246430/1‑149 (MQ=255)
                 aaaaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGt                                                                                                                               >  2:426350/1‑151 (MQ=255)
                   aaaaCCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGtat                                                                                                                             >  1:84700/1‑151 (MQ=255)
                         ccgccgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGt                                                                                                                       >  2:1021054/1‑151 (MQ=255)
                           gccgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGt                                                                                                                       >  2:514439/1‑149 (MQ=255)
                           gccgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTg                                                                                                                      >  1:550914/1‑150 (MQ=255)
                             cgGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGAcc                                                                                                                   >  1:19500/1‑151 (MQ=255)
                              ggAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGAcc                                                                                                                   >  1:568427/1‑150 (MQ=255)
                                   gAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAcc                                                                                                             >  2:958894/1‑151 (MQ=255)
                                        tttttGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACcagc                                                                                                          >  1:304360/1‑149 (MQ=255)
                                            tGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCagaag                                                                                                     <  2:553505/150‑1 (MQ=255)
                                             gTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCagaaga                                                                                                    >  1:128175/1‑150 (MQ=255)
                                             gTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCagaaga                                                                                                    <  2:817923/150‑1 (MQ=255)
                                             gTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAgaagaa                                                                                                   <  2:654615/151‑1 (MQ=255)
                                                    cAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCt                                                                                            <  1:407312/151‑1 (MQ=255)
                                                    cAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCt                                                                                            >  1:543985/1‑151 (MQ=255)
                                                     aaaGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCtt                                                                                           >  1:143083/1‑151 (MQ=255)
                                                     aaaGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCtt                                                                                           <  1:646731/151‑1 (MQ=255)
                                                      aaGCGAGGATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTg                                                                                          >  2:79223/1‑151 (MQ=255)
                                                         cGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGAtt                                                                                       <  1:988667/151‑1 (MQ=255)
                                                          gagaATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGAttt                                                                                      >  1:160440/1‑151 (MQ=255)
                                                          gagaATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGC‑AACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGAtttt                                                                                     >  1:211853/1‑151 (MQ=255)
                                                           agaATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGAtttt                                                                                     >  2:1063178/1‑151 (MQ=255)
                                                            gaATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGAtttt                                                                                     <  2:785076/150‑1 (MQ=255)
                                                            gaATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTtc                                                                                    >  1:789222/1‑151 (MQ=255)
                                                             aaTTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTtct                                                                                   >  2:262212/1‑151 (MQ=255)
                                                                tattTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCAc                                                                                >  2:557531/1‑151 (MQ=255)
                                                                 attTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTtctctc                                                                                <  2:883617/150‑3 (MQ=255)
                                                                 attTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCAc                                                                                <  1:816244/150‑1 (MQ=255)
                                                                   ttGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGa                                                                              >  2:56919/1‑150 (MQ=255)
                                                                    tGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGaa                                                                             <  2:480738/150‑1 (MQ=255)
                                                                    tGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAAt                                                                            >  2:977629/1‑151 (MQ=255)
                                                                     gATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTAGATTTTCTCACGAAt                                                                            >  2:601313/1‑150 (MQ=255)
                                                                      aTTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATAc                                                                          >  1:319808/1‑151 (MQ=255)
                                                                         ttCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCtt                                                                       >  2:815263/1‑151 (MQ=255)
                                                                           cGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCttt                                                                      >  1:215562/1‑150 (MQ=255)
                                                                               tcttTGTAGATCACGTGCTGGCGAACAACTGGATCGCATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGcc                                                                  >  2:243160/1‑150 (MQ=255)
                                                                                cttTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAt                                                                >  1:497067/1‑151 (MQ=255)
                                                                                 tttGTAGATCACGTGCTGGCGAACAACTGGATTGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑g                                                              >  1:799379/1‑151 (MQ=255)
                                                                                      aGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑Gattta                                                         <  1:246439/151‑1 (MQ=255)
                                                                                        aTCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑Gatttatt                                                       <  2:138531/151‑1 (MQ=255)
                                                                                         tCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑Gatttattt                                                      >  1:839223/1‑151 (MQ=255)
                                                                                          cACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑Gatttattt                                                      <  2:1000193/150‑1 (MQ=255)
                                                                                          cACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTc                                                     <  1:637503/151‑1 (MQ=255)
                                                                                           aCGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTcc                                                    >  2:33384/1‑151 (MQ=255)
                                                                                             gTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTTCGTTTGTTCTTCGTAGTGGTATAGAAGTTACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTACCCAT‑GATTTATTTCctc                                                  <  1:803510/151‑1 (MQ=255)
                                                                                             gTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCctc                                                  <  2:758234/151‑1 (MQ=255)
                                                                                             gTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCctc                                                  >  2:1052173/1‑151 (MQ=255)
                                                                                              tGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCctct                                                 >  2:202619/1‑151 (MQ=255)
                                                                                               gCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCctct                                                 <  2:649555/150‑1 (MQ=255)
                                                                                               gCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTa                                                <  2:557783/151‑1 (MQ=255)
                                                                                                 tGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTa                                                >  2:295172/1‑149 (MQ=255)
                                                                                                  ggCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAag                                              >  2:859091/1‑150 (MQ=255)
                                                                                                  ggCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTCGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAagt                                             >  2:871052/1‑151 (MQ=255)
                                                                                                     gAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAagtact                                          <  2:1012951/151‑1 (MQ=255)
                                                                                                        caacTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAagtacttag                                       <  2:260766/151‑1 (MQ=255)
                                                                                                         aacTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAagtacttagt                                      >  1:149150/1‑151 (MQ=255)
                                                                                                         aacTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAagtacttagt                                      <  2:535913/151‑1 (MQ=255)
                                                                                                             ggATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAagtacttagtac                                    <  1:661908/149‑1 (MQ=255)
                                                                                                               aTCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACtttt                                >  2:686565/1‑151 (MQ=255)
                                                                                                               aTCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACtttt                                <  1:964671/151‑1 (MQ=255)
                                                                                                                tCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTc                               >  2:640612/1‑151 (MQ=255)
                                                                                                                tCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTc                               >  2:728360/1‑151 (MQ=255)
                                                                                                                 cGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCg                              >  1:600793/1‑151 (MQ=255)
                                                                                                                    aTTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGcc                            >  1:543506/1‑150 (MQ=255)
                                                                                                                    aTTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCa                           <  1:539237/151‑1 (MQ=255)
                                                                                                                    aTTTTTTCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACC‑TTAGGCATGGTTTTATTTCTTTTAGGT‑CTTAGTACTTTTTGCCAc                          >  1:391103/1‑151 (MQ=255)
                                                                                                                       ttttCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACg                         <  1:179335/150‑1 (MQ=255)
                                                                                                                        tttCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACggg                       >  1:970396/1‑151 (MQ=255)
                                                                                                                         ttCAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACggg                       <  1:316532/150‑1 (MQ=255)
                                                                                                                           cAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCa                     <  2:262926/150‑1 (MQ=255)
                                                                                                                           cAGTTCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCAc                    >  1:66789/1‑151 (MQ=255)
                                                                                                                               tCCAGTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGc                  >  2:759566/1‑149 (MQ=255)
                                                                                                                                   gTTTTTCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCa            <  1:181931/151‑1 (MQ=255)
                                                                                                                                    tttttCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCa            <  2:485468/150‑1 (MQ=255)
                                                                                                                                     ttttCCGGCTTAGTGCGTTTGTTCTTCGTAGCGGTATAGTAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGc          >  1:424259/1‑151 (MQ=255)
                                                                                                                                      tttCCGGCTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGc          <  1:113207/150‑1 (MQ=255)
                                                                                                                                      tttCCGGCTTAGTGCG‑TTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGcc         >  2:936716/1‑150 (MQ=255)
                                                                                                                                             cTTAGTGCGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTg  >  1:159060/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCAT‑GATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTG  >  NC_000913/3811187‑3811478

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