Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,870,786 (C)5→6 coding (519/1050 nt) truD ← tRNA pseudouridine(13) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,870,7811.C98.2% 194.0 / ‑3.5 56coding (524/1050 nt)truDtRNA pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (27/28);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTCAACAACTTGATTAACGTCTGCTTTTTTGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGT  >  NC_000913/2870642‑2870911
                                                                                                                                            |                                                                                                                                    
gTCAACAACTTGATTAACGTCTGCTTTTTTGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTgcg                                                                                                                                       <  2:1309914/140‑1 (MQ=255)
gTCAACAACTTGATTAACGTCTCCTTTTTAGAGGCGCGCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTgcg                                                                                                                                       <  1:546621/140‑1 (MQ=255)
      aacTTGATTAACGTCTGCTTTTTTGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGcccccc                                                                                                                               <  1:269627/140‑2 (MQ=255)
      aacTTGATTAACGTCTGCTTTTTTGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGGATTCGTCTGCGCCGAGCGCTGCGcccccc                                                                                                                               >  1:719804/1‑139 (MQ=255)
           gATTAACGTCTGCTTTTTTGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCaa                                                                                                                          >  2:3158/1‑140 (MQ=255)
                aCGTCTGCTTTTTTGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGc                                                                                                                     >  1:760296/1‑140 (MQ=255)
                  gTCTGCTTTTTTGAGACGCTCAGCACCATCCTGATTACACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTa                                                                                                                   <  1:412599/140‑1 (MQ=255)
                   tCTGCTTTTTTGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTAc                                                                                                                  >  1:1286163/1‑140 (MQ=255)
                         tttttGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCaa                                                                                                            <  1:895495/140‑1 (MQ=255)
                            ttGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATcc                                                                                                         >  2:30199/1‑140 (MQ=255)
                              gAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCa                                                                                                       <  2:1473384/140‑1 (MQ=255)
                                 gcgcTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCaaaa                                                                                                    >  2:990298/1‑140 (MQ=255)
                                 gcgcTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGAATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCaaaa                                                                                                    >  2:224679/1‑140 (MQ=255)
                                    cTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCATGTGGCTACCGCCAATCCCAAAACGt                                                                                                 >  1:481901/1‑140 (MQ=255)
                                    cTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCaa                                                                                                      <  1:95709/135‑1 (MQ=255)
                                    cTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCaa                                                                                                      >  2:95709/1‑135 (MQ=255)
                                    cTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCa                                                                                                       >  2:1288164/1‑134 (MQ=255)
                                    cTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCa                                                                                                       <  1:1288164/134‑1 (MQ=255)
                                          aacaaTCTGATTAAACACCGCAGTGCGGGCTCCCGCCAAGCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCa                                                                                           <  2:481901/140‑1 (MQ=255)
                                            caaTCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCAcc                                                                                         >  2:826992/1‑140 (MQ=255)
                                              aTCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGa                                                                                       <  1:990298/140‑1 (MQ=255)
                                               tCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGaa                                                                                      >  1:1504656/1‑140 (MQ=255)
                                                 tGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAagt                                                                                    >  1:57055/1‑140 (MQ=255)
                                                 tGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAagt                                                                                    >  1:930286/1‑140 (MQ=255)
                                                   aTTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCACCCTGCAAGTTGGTCCGCCCAATCCAAAAATCTTGGGCAGCGAagtag                                                                                  >  2:751303/1‑140 (MQ=255)
                                                     aaaaCAACGCGCTGCGGACTGCAGACAACCAAAACTCCCCTTTTTTGCCATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGtt                                                                                <  2:547061/139‑1 (MQ=255)
                                                       aacaacGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCg                                                                              >  2:1411696/1‑140 (MQ=255)
                                                         caacGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTGCCGTCAACCCCAAAACATTGGCCACCGAAGTAGTTCGGt                                                                            >  2:417322/1‑140 (MQ=255)
                                                         caacGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGt                                                                            >  2:304421/1‑140 (MQ=255)
                                                           acGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGtt                                                                                >  2:818016/1‑134 (MQ=255)
                                                           acGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTGGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGtt                                                                                <  1:818016/134‑1 (MQ=255)
                                                                   cGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGGTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTa                                                                  >  1:888326/1‑140 (MQ=255)
                                                                         gCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCaaa                                                            >  2:1098674/1‑140 (MQ=255)
                                                                          ccGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAat                                                           <  1:1098674/140‑1 (MQ=255)
                                                                           cGACAACCAAAAACTCCGTTTTTTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAata                                                          <  2:760296/140‑1 (MQ=255)
                                                                           cGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAata                                                          <  1:304421/140‑1 (MQ=255)
                                                                           cGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAata                                                          <  1:368256/140‑1 (MQ=255)
                                                                           cGACAACCAAAAACTCCGTTTACTGCGATCGCGCACCCGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAata                                                          <  1:417322/140‑1 (MQ=255)
                                                                           cGACAACCAAAAACCCCGTTTATTGCGATCCCGCCCCGGAGTTTTGGTTTGCGCCCAGCGGTGGGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGGTTGGCCCCGAAGTAGTTCGGTACCCCTTTTACGCAAata                                                          <  2:1384824/140‑1 (MQ=255)
                                                                            gACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAatat                                                         <  1:805731/140‑1 (MQ=255)
                                                                               aaCCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGa                                                      >  2:171393/1‑140 (MQ=255)
                                                                                 ccAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATc                                                    <  2:57055/140‑1 (MQ=255)
                                                                                 ccAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATc                                                    <  1:1411696/140‑1 (MQ=255)
                                                                                 ccAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATc                                                    <  2:930286/140‑1 (MQ=255)
                                                                                    aaaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGa                                                 <  1:184272/140‑1 (MQ=255)
                                                                                    aaaaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGa                                                 <  1:146820/140‑1 (MQ=255)
                                                                                     aacccccGTTTATTGCGACCGCACACCGGAGTCTTGGTTTCCGCCCAGCGCTGCACCCCCCGCGAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGAc                                                <  2:152174/135‑1 (MQ=255)
                                                                                      aaCTCCGTTTATTGCGCGCGCGTACCGGAGTATGGCTTTACGCCCAGCCCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACg                                               <  2:719804/140‑1 (MQ=255)
                                                                                      aaCTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACg                                               <  1:766916/140‑1 (MQ=255)
                                                                                           cGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACgttgtt                                          >  2:314575/1‑140 (MQ=255)
                                                                                            gTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTc                                         >  1:277010/1‑140 (MQ=255)
                                                                                             tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATc                                                    >  2:1461005/1‑128 (MQ=255)
                                                                                             tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATc                                                    >  2:1341645/1‑128 (MQ=255)
                                                                                             tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATc                                                    <  1:1461005/128‑1 (MQ=255)
                                                                                             tttATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATc                                                    <  1:1341645/128‑1 (MQ=255)
                                                                                                       tCGCGCACCGGAGTATTGGTTTCCGCCCATCGCTGCGACCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATcgcg                             <  1:514517/140‑1 (MQ=255)
                                                                                                                       tGGTTTGCCCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgcg              <  1:224679/140‑1 (MQ=255)
                                                                                                                         gTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCa            <  2:816305/140‑1 (MQ=255)
                                                                                                                         gTTTGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCa            >  1:1379423/1‑140 (MQ=255)
                                                                                                                           ttGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGa          >  2:917702/1‑140 (MQ=255)
                                                                                                                            tGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgcg              >  2:128431/1‑135 (MQ=255)
                                                                                                                            tGCGCCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTcgcg              <  1:128431/135‑1 (MQ=255)
                                                                                                                             gcgcCCAGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAAc        <  2:277010/140‑1 (MQ=255)
                                                                                                                                   aGCGCTGCGCCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGt  <  1:30199/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                    
GTCAACAACTTGATTAACGTCTGCTTTTTTGAGGCGCTCAGCAACAATCTGATTAAACAACGCACTGCGGGCTGCCGACAACCAAAAACTCCGTTTATTGCGATCGCGCACCGGAGTATTGGTTTGCGCCCAGCGCTGCGCCCCCTGCAAGTTGCTACCGCCAATCCCAAAACGTTGGGCACCGAAGTAGTTCGGTACACCTTTTACGCAAATATCGATCAGACGTTGTTCAACGTCATCGCGATTGCTCACTTCGCGCAGAACCAGGGT  >  NC_000913/2870642‑2870911

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