Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,925,547 A→G 5.1% G29G (GGT→GGC mnmG ← 5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,925,5470AG 5.1% 237.0 / 3.0 79G29G (GGT→GGCmnmG5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base A (42/32);  new base G (2/2);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

TACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCGGGCCGTGGATTCTAC  >  NC_000913/3925399‑3925693
                                                                                                                                                    |                                                                                                                                                  
tACCCGCCGGATCGATCGCTTTCGCCATCAGACCTCCGAGGGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGa                                                                                                                                                    <  6:305/149‑1 (MQ=255)
   ccGCCTGCTCGCTCGCTTTCGCCATCAGACCGCCTAGGGCATCCACTGCTTTTACCAGATGGCCCTTCCCAATACCGCCGATCGCCGGGGTGCAGCTCATCTGCCCCAGAGTGGCGATATTGTGTGTCAAAAGCAGAGTCTGTTGAccc                                                                                                                                                 <  4:139077/149‑1 (MQ=255)
        tGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATAcg                                                                                                                                            >  4:170291/1‑149 (MQ=255)
             gatcgCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGcgccg                                                                                                                                       >  1:222251/1‑149 (MQ=255)
              atcgCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGcgccgc                                                                                                                                      <  5:90085/149‑1 (MQ=255)
               tcgCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCcgcg                                                                                                                                     >  3:199609/1‑149 (MQ=255)
                   tttCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTGCCCAATACCGCCGATCGCCGGGTTGCAGCGCATCTGCCCCAGAGTGTCGATATTGTGTGGCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCa                                                                                                                                 >  6:328036/1‑149 (MQ=255)
                         cATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCgg                                                                                                                           >  2:381644/1‑149 (MQ=255)
                              gcccgccgAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCg                                                                                                                      >  4:19694/3‑149 (MQ=255)
                               aCCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCgg                                                                                                                     <  3:316751/149‑1 (MQ=255)
                                ccgccgAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGt                                                                                                                    <  3:170291/149‑1 (MQ=255)
                                     gAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCCGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTg                                                                                                               <  4:331794/149‑1 (MQ=255)
                                     gAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTg                                                                                                               <  4:326301/149‑1 (MQ=255)
                                        tGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGGCCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCAt                                                                                                            >  8:29279/1‑149 (MQ=255)
                                        tGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCAt                                                                                                            <  2:222251/149‑1 (MQ=255)
                                          cATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGa                                                                                                          >  5:276590/1‑149 (MQ=255)
                                           aTCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGAc                                                                                                         <  3:19694/149‑1 (MQ=255)
                                            tCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATcgc                                                                                                                             <  7:415755/128‑1 (MQ=255)
                                            tCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATcgc                                                                                                                             >  8:415755/1‑128 (MQ=255)
                                            tCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCAt                                                                                                            >  2:271418/1‑145 (MQ=255)
                                            tCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCAt                                                                                                            <  1:271418/145‑1 (MQ=255)
                                              cACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGAcccc                                                                                                      >  5:346494/1‑149 (MQ=255)
                                               aCTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGAccccc                                                                                                     >  4:340311/1‑149 (MQ=255)
                                               aCTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGAccccc                                                                                                     >  2:89426/1‑149 (MQ=255)
                                                cttcttTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCg                                                                                                    >  1:157515/1‑149 (MQ=255)
                                                 ttcttTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGc                                                                                                   >  2:281203/1‑149 (MQ=255)
                                                   cttTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCcgcg                                                                                                                                     >  5:296233/1‑113 (MQ=255)
                                                   cttTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCcgcg                                                                                                                                     <  6:296233/113‑1 (MQ=255)
                                                    ttttACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCaa                                                                                                >  3:215698/1‑149 (MQ=255)
                                                      ttACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAat                                                                                               >  6:205390/1‑148 (MQ=255)
                                                      ttACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAat                                                                                               <  1:281203/148‑1 (MQ=255)
                                                        tcctgATGGCCCTTCCCAATAGGGCCGATCGCCGGGTTGGAGCGCTTCTGCACCTGAGTGTCGATATTGGGTGTCAAAAGCAGAGTCTGTTGACCGATACGCGCCGCGGACATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAatga                                                                                             <  6:346494/144‑1 (MQ=255)
                                                         ccAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGCCCCATACGCGCCGCGGCCCTCGCGGCCTCGGTGCCTGCATGCCCCCCGCCAatgat                                                                                            >  2:289761/1‑148 (MQ=255)
                                                         ccAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAatgat                                                                                            >  5:107365/1‑148 (MQ=255)
                                                          cAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAatgatga                                                                                          >  4:376305/1‑149 (MQ=255)
                                                           aGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAatgatgat                                                                                         <  4:199609/149‑1 (MQ=255)
                                                            gATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAatgatgat                                                                                         >  8:117184/1‑148 (MQ=255)
                                                             aTGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAAtgatgatga                                                                                       <  6:276590/149‑1 (MQ=255)
                                                              tGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGGCCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGAc                                                                                      <  7:29279/149‑1 (MQ=255)
                                                                tCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTGATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGGCTGCATGACCCCCGCCAATGATGCTGACGt                                                                                    >  8:35898/1‑149 (MQ=255)
                                                                 cccTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGt                                                                                    >  6:64021/1‑148 (MQ=255)
                                                                      cccAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGGCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCaaaa                                                                               <  6:157549/148‑1 (MQ=255)
                                                                        cAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGa                                                                            >  8:217218/1‑149 (MQ=255)
                                                                         aaTACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGAt                                                                           >  3:112063/1‑149 (MQ=255)
                                                                         aaTACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAACCAGAGTCTGTACACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGAt                                                                           >  5:133573/1‑149 (MQ=255)
                                                                          aTACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATc                                                                          <  1:89426/149‑1 (MQ=255)
                                                                             ccgccgATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAACAGCCGAGTCTGGTGACCCATCCGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGCCCCCCGCCAATGATGCTGACGTCAAAAGGATCCg                                                                        >  5:117205/1‑148 (MQ=255)
                                                                                  gATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATaa                                                                   >  1:317043/1‑148 (MQ=255)
                                                                                       ccGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATg                                                             <  3:340311/149‑1 (MQ=255)
                                                                                       ccGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATg                                                             <  6:107365/149‑1 (MQ=255)
                                                                                        cGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATgg                                                            <  3:376305/149‑1 (MQ=255)
                                                                                           gTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGa                                                         <  7:217218/149‑1 (MQ=255)
                                                                                            ttGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGAt                                                        <  4:112063/149‑1 (MQ=255)
                                                                                               cAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGc                                                     <  4:252542/149‑1 (MQ=255)
                                                                                                   tCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCg                                                 >  6:125023/1‑149 (MQ=255)
                                                                                                     aTCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCa                                               >  1:153163/1‑149 (MQ=255)
                                                                                                     aTCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCa                                               >  4:280689/1‑149 (MQ=255)
                                                                                                             cAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGAAGATGACGTCAAAAGGAGCCGGATAAAACATGGTGATTGCCTCGCATAACGCgg                                       >  4:114153/1‑149 (MQ=255)
                                                                                                                  tgtCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATg                                   >  3:369106/1‑148 (MQ=255)
                                                                                                                    tCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGaa                                 >  5:366100/1‑148 (MQ=255)
                                                                                                                          ttGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGa                           <  2:157515/148‑1 (MQ=255)
                                                                                                                            gtgtgtCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGAtt                         <  2:317043/148‑1 (MQ=255)
                                                                                                                              gtgtCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCACGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTg                        >  6:58500/1‑147 (MQ=255)
                                                                                                                               tgtCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCGCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGaa                      <  8:398947/148‑1 (MQ=255)
                                                                                                                                 tCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGcc                   >  8:184330/1‑149 (MQ=255)
                                                                                                                                   aaaaGCAGAGTCTGTTGGCCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCg                 >  5:379254/1‑149 (MQ=255)
                                                                                                                                   aaaaGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGCCCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCg                 >  7:72284/1‑149 (MQ=255)
                                                                                                                                   aaaaGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCg                 >  4:53576/1‑149 (MQ=255)
                                                                                                                                    aaaGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCaaa                                                                                <  5:177123/85‑1 (MQ=255)
                                                                                                                                    aaaGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCaaa                                                                                >  6:177123/1‑85 (MQ=255)
                                                                                                                                       gCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCGGGcc             >  2:364305/1‑149 (MQ=255)
                                                                                                                                            gTCTGTTGGCCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCGGGCCGTGGa        <  6:379254/149‑1 (MQ=255)
                                                                                                                                            gTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCGGGCCGTGGa        <  7:117184/149‑1 (MQ=255)
                                                                                                                                            gTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCACGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCGGGCCGTGGa        <  5:58500/149‑1 (MQ=255)
                                                                                                                                             tCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGccc                  >  1:104636/1‑138 (MQ=255)
                                                                                                                                             tCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGccc                  <  2:104636/138‑1 (MQ=255)
                                                                                                                                               tGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTTCTTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCACGGGCCGTGGATTc     >  4:235113/1‑149 (MQ=255)
                                                                                                                                                gTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCGGGCCGTGGATTc     >  4:292386/1‑148 (MQ=255)
                                                                                                                                                 ttGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCGGGCCGTGGATTc     <  5:125023/147‑1 (MQ=255)
                                                                                                                                                   gACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCATGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCGGGCCGTGGATTCTAc  <  1:381644/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TACCCGCCTGATCGATCGCTTTCGCCATCAGACCGCCGAGTGCATCCACTTCTTTTACCAGATGTCCCTTCCCAATACCGCCGATCGCCGGGTTGCAGCTCATCTGCCCCAGAGTGTCGATATTGTGTGTCAAAAGCAGAGTCTGTTGACCCATACGCGCCGCGGCCATCGCGGCCTCGGTGCCTGCATGACCCCCGCCAATGATGATGACGTCAAAAGGATCCGGATAAAACATGGTGATTGCCTCGCATAACGCGGTATGAAAATGGATTGAAGCCCGGGCCGTGGATTCTAC  >  NC_000913/3925399‑3925693

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