Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,717,186 G→A 5.1% R13W (CGG→TGG)  tyrS ← tyrosyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,717,1860GA 5.1% 265.6 / 3.3 78R13W (CGG→TGG) tyrStyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base A (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTAAAGCA  >  NC_000913/1717040‑1717334
                                                                                                                                                  |                                                                                                                                                    
gAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGct                                                                                                                                                    >  1:332124/1‑149 (MQ=255)
 aaGCGTTTCAGGCATAACAACGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGctc                                                                                                                                                   <  4:385010/149‑1 (MQ=255)
  aGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGctct                                                                                                                                                  <  4:232456/149‑1 (MQ=255)
         cAGGCATAACAATGGAACAAGACGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAAtt                                                                                                                                           <  1:699733/149‑1 (MQ=255)
          aGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCActc                                                                                                                                                   <  1:322024/140‑1 (MQ=255)
          aGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCActc                                                                                                                                                   >  2:70676/1‑140 (MQ=255)
          aGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCActc                                                                                                                                                   >  2:322024/1‑140 (MQ=255)
          aGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCActc                                                                                                                                                   <  1:70676/140‑1 (MQ=255)
             cATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGttt                                                                                                                                       >  2:679331/1‑149 (MQ=255)
             cATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGttt                                                                                                                                       <  4:606049/149‑1 (MQ=255)
              aTAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTa                                                                                                                                      >  3:177084/1‑149 (MQ=255)
                   aaTGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCaa                                                                                                                                 <  1:679331/149‑1 (MQ=255)
                    aTGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAg                                                                                                                                >  1:18126/1‑149 (MQ=255)
                     tGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGt                                                                                                                               <  3:706359/149‑1 (MQ=255)
                     tGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGt                                                                                                                               <  2:371530/149‑1 (MQ=255)
                        aaCAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTAc                                                                                                                            >  4:480071/1‑149 (MQ=255)
                             gATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCtt                                                                                                                       <  3:236438/149‑1 (MQ=255)
                              aTGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTg                                                                                                                      <  3:242417/149‑1 (MQ=255)
                                  cccAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCAt                                                                                                                  >  4:384951/1‑149 (MQ=255)
                                   ccAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATc                                                                                                                 >  1:69474/1‑149 (MQ=255)
                                    cAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATcc                                                                                                                >  4:543073/1‑148 (MQ=255)
                                    cAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCa                                                                                                                >  1:90117/1‑149 (MQ=255)
                                    cAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCa                                                                                                                >  3:478448/1‑149 (MQ=255)
                                            aGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTcc                                                                                                        <  2:561273/149‑1 (MQ=255)
                                             gCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCa                                                                                                       <  3:384951/149‑1 (MQ=255)
                                              cTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCAt                                                                                                      <  4:386372/149‑1 (MQ=255)
                                               tGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATg                                                                                                     >  3:708087/1‑149 (MQ=255)
                                                   aGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCAt                                                                                                      <  3:464772/144‑1 (MQ=255)
                                                   aGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTGGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCAt                                                                                                      >  4:464772/1‑144 (MQ=255)
                                                       gTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAAc                                                                                             >  1:602137/1‑149 (MQ=255)
                                                         aGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGa                                                                                           <  2:18126/149‑1 (MQ=255)
                                                          ggATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGAc                                                                                          <  4:84848/149‑1 (MQ=255)
                                                           gATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACt                                                                                         >  1:428012/1‑149 (MQ=255)
                                                            aTCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTg                                                                                        >  4:783537/1‑149 (MQ=255)
                                                            aTCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTg                                                                                        >  3:730185/1‑149 (MQ=255)
                                                            aTCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTg                                                                                        >  4:542633/1‑149 (MQ=255)
                                                            aTCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTg                                                                                        >  1:38358/1‑149 (MQ=255)
                                                                      cAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGcc                                                                              <  3:2608/149‑1 (MQ=255)
                                                                        aTAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCgg                                                                            >  1:460741/1‑149 (MQ=255)
                                                                                  aTCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACtt                                                                  <  4:730185/149‑1 (MQ=255)
                                                                                   tCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGc                                                                                                                     >  2:395354/1‑97 (MQ=255)
                                                                                   tCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGc                                                                                                                     <  1:395354/97‑1 (MQ=255)
                                                                                    cGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGATCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGCATATAACGACTGCACCTTTGCCGGTACACGACtttt                                                                >  4:472425/1‑149 (MQ=255)
                                                                                      ggCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCg                                                              >  2:516110/1‑149 (MQ=255)
                                                                                      ggCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCg                                                              >  3:537950/1‑149 (MQ=255)
                                                                                       gCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGc                                                             <  3:73555/149‑1 (MQ=255)
                                                                                         cTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCa                                                           >  1:244848/1‑149 (MQ=255)
                                                                                              gcCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGc                                                      <  2:332124/149‑1 (MQ=255)
                                                                                                  gTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGaaa                                                  >  4:223648/1‑149 (MQ=255)
                                                                                                    cGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAaga                                                <  2:69474/149‑1 (MQ=255)
                                                                                                          gCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATa                                          <  2:602137/149‑1 (MQ=255)
                                                                                                            tAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGaaa                                                  <  1:584823/139‑1 (MQ=255)
                                                                                                            tAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGaaa                                                  >  2:584823/1‑139 (MQ=255)
                                                                                                              aCGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAAt                                      >  3:193082/1‑149 (MQ=255)
                                                                                                              aCGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTCAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAAt                                      >  4:251206/1‑149 (MQ=255)
                                                                                                                gCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCGATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATaa                                    >  2:517222/1‑149 (MQ=255)
                                                                                                                 cTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATaaa                                   <  2:428012/149‑1 (MQ=255)
                                                                                                                    ccTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTTCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACAAAGAATAAAGTg                                >  1:482524/1‑149 (MQ=255)
                                                                                                                     cTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGc                               <  2:482524/149‑1 (MQ=255)
                                                                                                                      tCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGcc                              >  3:2913/1‑149 (MQ=255)
                                                                                                                      tCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGcc                              <  2:244848/149‑1 (MQ=255)
                                                                                                                           cgtcACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAAt                         <  3:480071/149‑1 (MQ=255)
                                                                                                                            gtcACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATc                        >  4:273154/1‑149 (MQ=255)
                                                                                                                             tcACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCa                       <  4:708087/149‑1 (MQ=255)
                                                                                                                                ccTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATTTAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGa                    <  3:542633/149‑1 (MQ=255)
                                                                                                                                ccTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGa                    <  4:177084/149‑1 (MQ=255)
                                                                                                                                ccTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGa                    <  4:2913/149‑1 (MQ=255)
                                                                                                                                         ccAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCgg                                                                            <  1:723237/84‑1 (MQ=255)
                                                                                                                                         ccAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCgg                                                                            >  2:723237/1‑84 (MQ=255)
                                                                                                                                         ccAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCg           >  2:286961/1‑149 (MQ=255)
                                                                                                                                           aGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCCCCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGAt         <  3:543073/149‑1 (MQ=255)
                                                                                                                                           aGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGAt         <  4:537950/149‑1 (MQ=255)
                                                                                                                                               cccGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGAc                                                                    >  3:767553/1‑86 (MQ=255)
                                                                                                                                               cccGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGAc                                                                    <  4:767553/86‑1 (MQ=255)
                                                                                                                                               cccGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTaaa     >  1:605749/1‑149 (MQ=255)
                                                                                                                                               cccGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTaaa     >  3:637827/1‑149 (MQ=255)
                                                                                                                                                ccGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAAt                                      <  2:157510/115‑1 (MQ=255)
                                                                                                                                                ccGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAAt                                      >  1:157510/1‑115 (MQ=255)
                                                                                                                                                  gCTCTTGCGATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTAAAGCa  <  1:517222/149‑1 (MQ=255)
                                                                                                                                                  gCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTAAAGCa  <  3:273154/149‑1 (MQ=255)
                                                                                                                                                  gCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTAAAGCa  >  4:380721/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAAAGAGACATAGAATAAAGTGCCAGAATCAGGAGTACCAGCGATTAAAGCA  >  NC_000913/1717040‑1717334

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