Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,034,756 C→T 5.2% C140Y (TGC→TAC)  lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,034,7560CT 5.2% 245.5 / 2.6 78C140Y (TGC→TAC) lysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base T (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

CGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGC  >  NC_000913/3034611‑3034894
                                                                                                                                                 |                                                                                                                                          
cGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGt                                                                                                                                         <  8:86333/149‑1 (MQ=255)
cGAGCGAACTTTTAAGGTGTTGCGGGATTCATCGGTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCTGCGGACGCAGTGCTTTGGTCAGCAGACGCAAGTCGGTGCAGt                                                                                                                                         <  5:260674/149‑1 (MQ=255)
    cgcACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGa                                                                                                                                      >  3:296760/1‑148 (MQ=255)
      cACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTgg                                                                                                                                       <  1:244060/145‑1 (MQ=255)
      cACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTgg                                                                                                                                       >  2:244060/1‑145 (MQ=255)
            aaaGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCTCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGAGCAGTGGATAGACAGtt                                                                                                                             >  6:245430/1‑149 (MQ=255)
             aaGGTGTTGCGGTATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTGAGCAGACGCAACTCGGTGCAGTGGATAGACAGt                                                                                                                              >  5:89004/1‑147 (MQ=255)
             aaGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGAt                                                                                                                                     <  7:269259/140‑1 (MQ=255)
             aaGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGAt                                                                                                                                     >  8:269259/1‑140 (MQ=255)
             aaGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATACCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGtt                                                                                                                             <  1:162023/148‑1 (MQ=255)
              aGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGTAGTGGATAGACAGTTCg                                                                                                                           >  8:49438/1‑149 (MQ=255)
                 tgtTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACTCAACTCGGTGCAGTGGATAGACAGTTCGCCg                                                                                                                        >  4:237479/1‑149 (MQ=255)
                     gCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttt                                                                                                                    >  2:279303/1‑149 (MQ=255)
                       gggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttg                                                                                                                  >  4:252436/1‑149 (MQ=255)
                       gggATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGACGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAACACACGCAACTCGCTGCAGTGGATAGACAGTTCGCCggttttg                                                                                                                  >  2:339432/1‑149 (MQ=255)
                          aTTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGGAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggtt                                                                                                               <  1:327025/149‑1 (MQ=255)
                          aTTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggtt                                                                                                               <  2:266294/149‑1 (MQ=255)
                           ttCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGGAGACTCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttggttt                                                                                                              <  3:237479/149‑1 (MQ=255)
                           ttCATAGTGGGATATGAGATCGAGATAAAGCTGACGATAGCGCGCTTCCTTATACTGCAAGCCGTGTAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACTCAACTCTGTGCAGTGGATAGACAGTTCTCCggttttggttt                                                                                                              <  6:89004/149‑1 (MQ=255)
                              aTCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAg                                                                                                                                          <  1:64461/118‑1 (MQ=255)
                              aTCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAg                                                                                                                                          >  2:64461/1‑118 (MQ=255)
                              aTCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttg                                                                                                                  <  7:119638/142‑1 (MQ=255)
                              aTCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttg                                                                                                                  >  8:119638/1‑142 (MQ=255)
                                     gagatgagatCGAGATAACGCTGACGATAGGGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGAAGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGGAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCt                                                                                                     <  5:245430/148‑1 (MQ=255)
                                           agatCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCtt                                                                                               <  1:279303/148‑1 (MQ=255)
                                             atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGAc                                                                                                                                 <  8:137578/112‑1 (MQ=255)
                                             atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGAc                                                                                                                                 >  7:137578/1‑112 (MQ=255)
                                             atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTcg                                                                                            >  7:223809/1‑149 (MQ=255)
                                             atCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCATTcg                                                                                            >  2:334674/1‑149 (MQ=255)
                                                gagaTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGTAGTGGATAGACAGTTCGCCggttttggtttt                                                                                                             >  2:17557/1‑129 (MQ=255)
                                                gagaTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGTAGTGGATAGACAGTTCGCCggttttggtttt                                                                                                             <  1:17557/129‑1 (MQ=255)
                                                   aTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGcc                                                                                        >  2:39345/1‑147 (MQ=255)
                                                   aTAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGAGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGcc                                                                                        >  6:222319/1‑147 (MQ=255)
                                                    tAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGTAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGa                                                                                      <  7:49438/148‑1 (MQ=255)
                                                     aaCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttgg                                                                                                                 <  8:338674/120‑1 (MQ=255)
                                                     aaCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttgg                                                                                                                 >  7:338674/1‑120 (MQ=255)
                                                      aCGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGCGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGa                                                                                   <  2:213444/149‑1 (MQ=255)
                                                        gCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATg                                                                                 >  8:214007/1‑149 (MQ=255)
                                                         cTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGt                                                                                <  4:296760/149‑1 (MQ=255)
                                                                   cgcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTc                                                                      <  5:222319/149‑1 (MQ=255)
                                                                    gcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg                                                                                     >  4:217282/1‑133 (MQ=255)
                                                                    gcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAg                                                                                     <  3:217282/133‑1 (MQ=255)
                                                                    gcgcTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGATCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCAGAGGATGTCGCCGAGGTcc                                                                     >  4:76326/1‑149 (MQ=255)
                                                                        ttCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGAAGTTCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGCTTTTGTGTTTGAACAGCTTCACTTTCGCGCCGAGGATGTCTCCTCGGTCACAtt                                                                 >  5:239572/1‑149 (MQ=255)
                                                                           cTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTGTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCAttttt                                                              <  6:7328/149‑1 (MQ=255)
                                                                            tGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttt                                                                                                                    <  7:309856/94‑1 (MQ=255)
                                                                            tGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttt                                                                                                                    >  8:309856/1‑94 (MQ=255)
                                                                            tGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCAtttttt                                                             >  1:99465/1‑149 (MQ=255)
                                                                                 cTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAAc                                                         <  7:214007/148‑1 (MQ=255)
                                                                                    ccaGCCGTTGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTGTGTGCAGTGGATAGACAGTGCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGc                                                      <  6:239572/146‑1 (MQ=255)
                                                                                    cAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGc                                                      <  4:160864/148‑1 (MQ=255)
                                                                                         cGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttg                                                                                                                  <  7:91814/83‑1 (MQ=255)
                                                                                         cGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCggttttg                                                                                                                  >  8:91814/1‑83 (MQ=255)
                                                                                         cGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGcc                                                                            <  5:147778/121‑1 (MQ=255)
                                                                                         cGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGcc                                                                            >  6:147778/1‑121 (MQ=255)
                                                                                             gAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGGTTTGACCAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATaaa                                            <  1:334674/149‑1 (MQ=255)
                                                                                               aTTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACg                                          >  1:126433/1‑149 (MQ=255)
                                                                                                       ggCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGa                                                                                                           >  4:78897/1‑76 (MQ=255)
                                                                                                       ggCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGa                                                                                                           <  3:78897/76‑1 (MQ=255)
                                                                                                           gCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGgaga                              <  1:339432/149‑1 (MQ=255)
                                                                                                             ggACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATc                            <  2:99465/149‑1 (MQ=255)
                                                                                                                  cAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCAc                       <  1:39345/149‑1 (MQ=255)
                                                                                                                   aGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACg                      <  8:223809/149‑1 (MQ=255)
                                                                                                                     tGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGGTCGCCGGTTTTGGGTTTTAACAGCTTACCTTTCGCGCCGAGGATGTCGCCTAGGTCACATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACggg                    <  3:76326/149‑1 (MQ=255)
                                                                                                                        tttGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCaa                 >  3:212774/1‑149 (MQ=255)
                                                                                                                            gTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTa             <  4:212774/149‑1 (MQ=255)
                                                                                                                            gTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTa             >  8:37237/1‑149 (MQ=255)
                                                                                                                               agcagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATc                            >  1:258910/1‑131 (MQ=255)
                                                                                                                               agcagACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATc                            <  2:258910/131‑1 (MQ=255)
                                                                                                                                   gACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATaa                                             <  4:251570/110‑1 (MQ=255)
                                                                                                                                   gACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATaa                                             >  3:251570/1‑110 (MQ=255)
                                                                                                                                     cGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAAt    <  2:126433/149‑1 (MQ=255)
                                                                                                                                      gCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCGTGTTTTGAAGTGGTCGTTATAACCGCCTTGCGGGAGATGGGCACGGGCAGCGTAGAGCTGAATg   >  2:317669/1‑149 (MQ=255)
                                                                                                                                      gCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGaa                                                          >  3:25361/1‑94 (MQ=255)
                                                                                                                                      gCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGaa                                                          <  4:25361/94‑1 (MQ=255)
                                                                                                                                        aaCTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCgg                                  >  4:312469/1‑116 (MQ=255)
                                                                                                                                        aaCTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCgg                                  <  3:312469/116‑1 (MQ=255)
                                                                                                                                        aaCTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGc  >  2:339050/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
CGAGCGCACTTTAAAGGTGTTGCGGGATTCATCGTTGGAGATGAGATCGAGATAACGCTGACGATAGCGCGCTTCCTGATCCTGCAAGCCGTGGAATTTATCCGGCAGCGGACGCAGTGCTTTGGTCAGCAGACGCAACTCGGTGCAGTGGATAGACAGTTCGCCGGTTTTGGTTTTGAACAGCTTACCTTTCGCGCCGAGGATGTCGCCGAGGTCCCATTTTTTGAACTGCTCGTTATAAACGCCTTCCGGGAGATCGTCACGGGCAACGTACAGCTGAATGC  >  NC_000913/3034611‑3034894

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