Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,745,237 T→C S35G (AGC→GGC)  trmD ← tRNA m(1)G37 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,745,2370TC97.1% 234.0 / ‑2.0 69S35G (AGC→GGC) trmDtRNA m(1)G37 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (38/29);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

CACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGTCTTTTACCGTTTATCCGGTGGTTTA  >  NC_000913/2745102‑2745372
                                                                                                                                       |                                                                                                                                       
cccccGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTgg                                                                                                                                     <  1:960481/138‑1 (MQ=255)
cACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTgg                                                                                                                                     >  2:1050217/1‑140 (MQ=255)
cACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTgg                                                                                                                                     >  1:746714/1‑140 (MQ=255)
         cgGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTcagc                                                                                                                            >  1:60813/1‑140 (MQ=255)
         cgGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTcagc                                                                                                                            >  2:144551/1‑140 (MQ=255)
            cTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAgg                                                                                                                         <  2:73403/140‑1 (MQ=255)
                           tGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGcc                                                                                                          <  2:1033428/140‑1 (MQ=255)
                            ggCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGccc                                                                                                         >  2:812933/1‑140 (MQ=255)
                               gTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGc                                                                                                      >  1:109664/1‑140 (MQ=255)
                                tCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGcc                                                                                                     >  1:243976/1‑140 (MQ=255)
                                tCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGcc                                                                                                     >  2:626951/1‑140 (MQ=255)
                                   cGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGt                                                                                                  <  2:852610/140‑1 (MQ=255)
                                   cGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGt                                                                                                  >  2:329625/1‑140 (MQ=255)
                                   cGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGt                                                                                                  <  2:746714/140‑1 (MQ=255)
                                       aGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTAcc                                                                                              <  1:335130/140‑1 (MQ=255)
                                          gTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCg                                                                                           <  2:243976/140‑1 (MQ=255)
                                          gTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCg                                                                                           <  1:224062/140‑1 (MQ=255)
                                          gTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCg                                                                                           <  1:144551/140‑1 (MQ=255)
                                             gCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTaa                                                                                        >  1:491639/1‑140 (MQ=255)
                                                    catTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTaa                                                                                 <  2:491639/140‑1 (MQ=255)
                                                           aTCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGcgcg                                                                          >  2:858417/1‑140 (MQ=255)
                                                           aTCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGcgcg                                                                          >  1:327830/1‑140 (MQ=255)
                                                             ccccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGa                                                                        >  1:619413/1‑140 (MQ=255)
                                                             ccccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGa                                                                        >  1:726143/1‑140 (MQ=255)
                                                             ccccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGa                                                                        >  2:515115/1‑140 (MQ=255)
                                                             ccccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGa                                                                        >  1:487027/1‑140 (MQ=255)
                                                             ccccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGa                                                                        >  2:709616/1‑140 (MQ=255)
                                                             ccccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGa                                                                        >  2:831534/1‑140 (MQ=255)
                                                             ccccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGa                                                                        >  1:942181/1‑140 (MQ=255)
                                                              cccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGaa                                                                       <  2:284457/140‑1 (MQ=255)
                                                              cccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGaa                                                                       >  1:1038694/1‑140 (MQ=255)
                                                              cccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGaa                                                                       <  1:1050217/140‑1 (MQ=255)
                                                              cccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGaa                                                                       >  1:464913/1‑140 (MQ=255)
                                                               ccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAAc                                                                      >  2:991273/1‑140 (MQ=255)
                                                               ccGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAAc                                                                      >  2:825119/1‑140 (MQ=255)
                                                                cGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACa                                                                     <  1:460091/140‑1 (MQ=255)
                                                                 ggTCCGCCGCCGTAAGGACAATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACAt                                                                    >  2:986362/1‑140 (MQ=255)
                                                                     cgccgccgTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCa                                                                <  2:60813/140‑1 (MQ=255)
                                                                         gccgTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGaa                                                            <  2:487027/140‑1 (MQ=255)
                                                                              aaGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAgg                                                       <  1:825119/140‑1 (MQ=255)
                                                                              aaGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAgg                                                       <  2:109664/140‑1 (MQ=255)
                                                                               aGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCt                                                     >  2:603591/1‑140 (MQ=255)
                                                                                     aTCGTCCACGGTACGGTGCCTGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTa                                                <  1:1023520/140‑1 (MQ=255)
                                                                                       cGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATg                                              <  1:812933/140‑1 (MQ=255)
                                                                                       cGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATg                                              <  2:327830/140‑1 (MQ=255)
                                                                                       cGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATg                                              <  1:592639/140‑1 (MQ=255)
                                                                                         tCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGcc                                            >  1:655545/1‑140 (MQ=255)
                                                                                                tacggtGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCcac                                     >  1:245078/1‑140 (MQ=255)
                                                                                                tacggtGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCcac                                     >  1:392080/1‑140 (MQ=255)
                                                                                                tacggtGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCcac                                     >  1:627827/1‑140 (MQ=255)
                                                                                                tacggtGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCcac                                     >  1:846866/1‑140 (MQ=255)
                                                                                                     tGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAgc                                >  1:821421/1‑140 (MQ=255)
                                                                                                      gCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAgcg                               <  1:991273/140‑1 (MQ=255)
                                                                                                       ccGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAgcgc                              >  1:959055/1‑140 (MQ=255)
                                                                                                        cGGTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGcc                             >  1:857355/1‑140 (MQ=255)
                                                                                                          gTCATGCGTGAAGTCGCGAGGCCTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGt                           >  1:936935/1‑140 (MQ=255)
                                                                                                          gTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCTGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGt                           >  2:598155/1‑140 (MQ=255)
                                                                                                          gTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGt                           >  1:158087/1‑140 (MQ=255)
                                                                                                          gTCATGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGt                           >  2:674968/1‑140 (MQ=255)
                                                                                                              tGCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGTCttt                       <  1:858417/140‑1 (MQ=255)
                                                                                                               gCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAgcg                               <  1:536133/131‑1 (MQ=255)
                                                                                                               gCGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAgcg                               >  2:536133/1‑131 (MQ=255)
                                                                                                                cGTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGTCTTTTa                     <  2:464913/140‑1 (MQ=255)
                                                                                                                 gTGAAGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGTCTTTTAc                    <  2:821421/140‑1 (MQ=255)
                                                                                                                    aaGTCGCGAGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGTCTTTTACCGt                 <  2:619413/140‑1 (MQ=255)
                                                                                                                            aGGACTCCAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGTCTTTTACCGTTTATCCgg         <  2:392080/140‑1 (MQ=255)
                                                                                                                            aGGACTCCAGCCCTGGATGCTCAGCAGGCCATATTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATACCAATCCACATAGCGCCGTCTTTTACCGTTTATCCgg         <  2:959055/140‑1 (MQ=255)
                                                                                                                                   cAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGTCTTTTACCGTTTATCCGGTGGTTTa  >  1:793401/1‑140 (MQ=255)
                                                                                                                                   cAGCCCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGTCTTTTACCGTTTATCCGGTGGTTTa  >  2:215550/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCGTAAGGACGATCGTCCACGGTACGGTGCCGGTCATGCGTGAAGTCGCGAGGACTCCAGCTCTGGATGCTCAGCAGGCCATTTTTAACTGCCCGGCCAGTTACCCCGTAATCGGTAATTGCGCGGAACATTTCAGGAAACAGGCTAATTATGCCAATCCACATAGCGCCGTCTTTTACCGTTTATCCGGTGGTTTA  >  NC_000913/2745102‑2745372

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