Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,125,059 A→C 6.0% R584R (CGT→CGG priA ← Primosome factor n' (replication factor Y)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,125,0590AC 6.0% 198.1 / 2.1 63R584R (CGT→CGGpriAPrimosome factor n' (replication factor Y)
Reads supporting (aligned to +/- strand):  ref base A (25/34);  new base C (2/2);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCTGCATCATTCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTAC  >  NC_000913/4124921‑4125198
                                                                                                                                          |                                                                                                                                           
tAGCTGCATCATTCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGa                                                                                                                                   >  6:118975/1‑149 (MQ=255)
tAGCTGCATCATTCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGa                                                                                                                                   <  3:243809/149‑1 (MQ=255)
  gCTGCATCATTCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCCACAGCGGATGTTCCGGATGGTGCGTTGGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCCCGACCGGCGAcc                                                                                                                                 >  1:165775/1‑149 (MQ=255)
  gCTGCATCATTCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACACGCACGCCCGGCCCGACCGGCGAcc                                                                                                                                 >  1:171054/1‑149 (MQ=255)
      catcatTCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTgg                                                                                                                              <  7:190756/148‑1 (MQ=255)
      catcatTCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTgg                                                                                                                              <  7:296901/148‑1 (MQ=255)
           ttCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTaa                                                                                                                        <  3:318555/149‑1 (MQ=255)
            tCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTa                                                                                                                         >  7:55709/1‑147 (MQ=255)
              gccgcTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGc                                                                                                                     <  1:373177/149‑1 (MQ=255)
                 gcTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGa                                                                                                                  <  5:39611/149‑1 (MQ=255)
                   tCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGc                                                                                                                >  5:362794/1‑149 (MQ=255)
                          gcgcCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGCGCGAAACGc                                                                                                         >  5:242118/1‑149 (MQ=255)
                             ccTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGctct                                                                                                      >  5:16894/1‑149 (MQ=255)
                                gTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGAGTAAAGCTGAGCGAAACGCTCTGcc                                                                                                   >  3:241705/1‑149 (MQ=255)
                                   cGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAg                                                                                                >  5:327229/1‑149 (MQ=255)
                                     gCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGc                                                                                               >  4:142748/1‑148 (MQ=255)
                                        aaGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGaaa                                                                                           <  5:195199/149‑1 (MQ=255)
                                          ggCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGa                                                                                                                  >  2:174496/1‑124 (MQ=255)
                                          ggCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGa                                                                                                                  <  1:174496/124‑1 (MQ=255)
                                           gCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATc                                                                                        <  3:142748/149‑1 (MQ=255)
                                           gCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATc                                                                                        <  6:327229/149‑1 (MQ=255)
                                                gtAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCa                                                                                    <  5:253451/148‑1 (MQ=255)
                                                   gCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGGGCGAAAATCGGCAGaa                                                                                 <  8:222012/148‑1 (MQ=255)
                                                    ccTTTATAGAGCAACGGTTGCAACAGAGGATGTTCCGGATGGTGCGTGTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCCCGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGaaaa                                                                               <  2:165775/149‑1 (MQ=255)
                                                    ccTTAATAGAGCAACGATTGCAACAGAGGATGTTCCGGATGGTGCGTTAGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCCCGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGaaaa                                                                               <  2:171054/149‑1 (MQ=255)
                                                     cTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTgg                                                                                                                              >  5:131514/1‑101 (MQ=255)
                                                     cTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTgg                                                                                                                              <  6:131514/101‑1 (MQ=255)
                                                     cTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGaaaa                                                                               >  2:326942/1‑148 (MQ=255)
                                                      tttATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACa                                                                             <  3:156714/149‑1 (MQ=255)
                                                      tttATAGAGCAACGATTGCAACAGAGGATGTTCCGGATGGTGCGATTGCAGCACAACAACGCACTGATTACCCCCACGCCCGGCACGACCGGCGACCAGGGTGTAAAGCTCAGCGAAACGCTCAGCCGAGCGAAAATCGGCAGaaaa                                                                               >  5:91700/1‑147 (MQ=255)
                                                        tataGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACa                                                                             <  8:55709/147‑1 (MQ=255)
                                                         ataGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAgc                                                                           <  5:118975/148‑1 (MQ=255)
                                                           agagCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAgcgc                                                                         <  7:92198/148‑1 (MQ=255)
                                                                   gTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCcgtccac                                                                <  6:16894/149‑1 (MQ=255)
                                                                    tttGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAgc                                                                           >  6:344261/1‑137 (MQ=255)
                                                                    tttGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAgc                                                                           <  5:344261/137‑1 (MQ=255)
                                                                         aaCAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGc                                                                                                                >  8:83159/1‑95 (MQ=255)
                                                                         aaCAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGc                                                                                                                <  7:83159/95‑1 (MQ=255)
                                                                         aaCAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAg                                                          >  1:80338/1‑149 (MQ=255)
                                                                          aCAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGa                                                                                                                                   <  8:306370/75‑1 (MQ=255)
                                                                            agagGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCACCCTGTTTACCCGCCCGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCcgtccacgt                                                              >  6:163209/1‑142 (MQ=255)
                                                                            agagGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCACCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCcgtccacgt                                                              <  5:163209/142‑1 (MQ=255)
                                                                              agGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATg                                                     >  8:40386/1‑149 (MQ=255)
                                                                                  tGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAAc                                                 <  7:40386/149‑1 (MQ=255)
                                                                                  tGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAAc                                                 >  8:233764/1‑149 (MQ=255)
                                                                                   gTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAAcc                                                >  5:239136/1‑149 (MQ=255)
                                                                                    ttCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCa                                               <  1:187082/149‑1 (MQ=255)
                                                                                       cGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGc                                             <  2:80338/148‑1 (MQ=255)
                                                                                           tGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGTAACCAGCGTcac                                        <  1:326942/149‑1 (MQ=255)
                                                                                                      cagcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAAt                             <  6:239136/149‑1 (MQ=255)
                                                                                                       agcaCCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATg                            <  3:244430/149‑1 (MQ=255)
                                                                                                         caccacTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGt                          >  6:223880/1‑149 (MQ=255)
                                                                                                           ccacTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGa                        <  6:242118/149‑1 (MQ=255)
                                                                                                            cacTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACACTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGAc                       >  7:88215/1‑149 (MQ=255)
                                                                                                                  gCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGAGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCGTTc                  <  4:241705/148‑1 (MQ=255)
                                                                                                                   cccTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGc                >  4:43478/1‑149 (MQ=255)
                                                                                                                     cTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCa              <  5:223880/149‑1 (MQ=255)
                                                                                                                      tGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTg                         <  3:274901/137‑1 (MQ=255)
                                                                                                                      tGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTg                         >  4:274901/1‑137 (MQ=255)
                                                                                                                          tACCCGCCAGCCCGGCACGACCGGCGACCTGGGGGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCAt          <  6:362794/148‑1 (MQ=255)
                                                                                                                             ccGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTtgt      >  6:20776/1‑149 (MQ=255)
                                                                                                                               gCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCAAGTAATGCAACCAGCGTCACATCCGTGAAATGGTGACCTTTCGCCAGAATTtgtg     >  2:334420/1‑148 (MQ=255)
                                                                                                                                 aCGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTAc  <  3:43478/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
TAGCTGCATCATTCGCCGCTCAGCCAGCGCCTGTTCGGCAAAGGCGTCGTAGCCTTTATAGAGCAACGTTTGCAACAGAGGATGTTCCGGATGGTGCGTTTGCAGCACCACTTCGCCCTGTTTACCCGCACGCCCGGCACGACCGGCGACCTGGGTGTAAAGCTGAGCGAAACGCTCTGCCGAGCGAAAATCGGCAGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTAC  >  NC_000913/4124921‑4125198

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