Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,483,659 T→A 5.1% D60V (GAT→GTT) ‡ valS ← valyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,483,6590TA 5.1% 321.8 / 5.6 97D60V (GAT→GTT) ‡valSvalyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (43/49);  new base A (3/2);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

GCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGTAAAGCGGCTGTTC  >  NC_000913/4483512‑4483807
                                                                                                                                                   |                                                                                                                                                    
gCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATc                                                                                                                                                     <  5:369648/149‑1 (MQ=255)
gCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATc                                                                                                                                                     <  4:110591/149‑1 (MQ=255)
  ttCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATcaa                                                                                                                                                   >  4:286256/1‑147 (MQ=255)
    cgcgGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCa                                                                                                                                                   <  4:354003/147‑1 (MQ=255)
     gcgGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatga                                                                                                                                                >  3:117055/1‑149 (MQ=255)
      cgGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgat                                                                                                                                               <  3:28079/149‑1 (MQ=255)
       ggCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgatg                                                                                                                                              <  7:220840/149‑1 (MQ=255)
       ggCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgatg                                                                                                                                              <  5:376579/149‑1 (MQ=255)
         ccGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGt                                                                                                                                            <  7:219587/149‑1 (MQ=255)
          cGTAGTCGTGACGGGTTGTTCCTTCTTCTGCGGCTATCTTGAGCTCAACGACCATGTGGTTAGCGATCCCGGCGTGGTCAGGACCGACCTGCCACAGGGTGTTTGTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGtt                                                                                                                                           <  7:369720/149‑1 (MQ=255)
            tAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTg                                                                                                                                         >  6:341555/1‑149 (MQ=255)
              gTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCt                                                                                                                                       >  8:249627/1‑149 (MQ=255)
                   gACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGAGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCa                                                                                                                                                   <  7:62508/132‑1 (MQ=255)
                   gACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGAGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCa                                                                                                                                                   >  8:62508/1‑132 (MQ=255)
                     cGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                                >  3:75454/1‑149 (MQ=255)
                     cGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGGTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                                <  5:224729/149‑1 (MQ=255)
                     cGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGGCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                                <  1:324420/149‑1 (MQ=255)
                       ggTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGt                                                                                                                              <  1:237281/149‑1 (MQ=255)
                         ttttACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGa                                                                                                                            <  6:90814/149‑1 (MQ=255)
                              ccTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCat                                                                                                                       <  1:116448/149‑1 (MQ=255)
                               cttcttCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCat                                                                                                                                                  <  5:13724/121‑1 (MQ=255)
                               cttcttCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCat                                                                                                                                                  >  6:13724/1‑121 (MQ=255)
                                  cttctGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCatat                                                                                                                     <  5:320088/147‑1 (MQ=255)
                                  cttctGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATg                                                                                                                    >  2:98527/1‑148 (MQ=255)
                                   ttctGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGc                                                                                                                   >  8:121708/1‑148 (MQ=255)
                                      tGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGa                                                                                                                                    >  4:205652/1‑128 (MQ=255)
                                         ggCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGc                                                                                                            >  1:185231/1‑149 (MQ=255)
                                          gCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                                >  1:181679/1‑128 (MQ=255)
                                          gCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGc                                                                                                                                <  2:181679/128‑1 (MQ=255)
                                          gCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGcc                                                                                                           >  4:149945/1‑149 (MQ=255)
                                                   cgcTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGt                                                                                                   <  3:149945/148‑1 (MQ=255)
                                                    gcTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTc                                                                                                 >  6:218155/1‑149 (MQ=255)
                                                     cTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGACCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcg                                                                                                >  8:81528/1‑149 (MQ=255)
                                                      tCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcgg                                                                                               >  2:233363/1‑149 (MQ=255)
                                                      tCAACGACCATCTGGGTAGCGATACCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcgg                                                                                               >  7:32313/1‑149 (MQ=255)
                                                            aCCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTcggcggcgg                                                                                         <  4:117055/149‑1 (MQ=255)
                                                               aTCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGAccc                                                                                                                         <  8:244700/114‑1 (MQ=255)
                                                               aTCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGAccc                                                                                                                         >  7:244700/1‑114 (MQ=255)
                                                                    ggTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatga                                                                                 <  4:276611/149‑1 (MQ=255)
                                                                    ggTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCACATGCAAACTGCCGGTGACGTTCGGCGGCGGGATc                                                                                     <  4:257295/145‑1 (MQ=255)
                                                                     gTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatgat                                                                                <  7:249627/149‑1 (MQ=255)
                                                                     gTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGGTCGGCGGCGGGATCatgat                                                                                <  1:233363/149‑1 (MQ=255)
                                                                     gTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTATGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatga                                                                                 <  4:275975/148‑1 (MQ=255)
                                                                     gTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTATGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatga                                                                                 >  3:275975/1‑148 (MQ=255)
                                                                      tAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCatgatg                                                                               >  6:295689/1‑149 (MQ=255)
                                                                        gCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGc                                                                              <  4:75454/148‑1 (MQ=255)
                                                                          gATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAg                                                                            >  4:36461/1‑148 (MQ=255)
                                                                                 gCGTGGTCAGTACCGACCTGCGACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGTATCATGATGCAGAAACTTTc                                                                    <  8:32313/149‑1 (MQ=255)
                                                                                 gCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTc                                                                    <  4:325663/149‑1 (MQ=255)
                                                                                 gCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTc                                                                    <  3:36461/149‑1 (MQ=255)
                                                                                   gTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCt                                                                  <  2:185231/149‑1 (MQ=255)
                                                                                      gTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgatg                                                                                                                                              >  6:262114/1‑70 (MQ=255)
                                                                                      gTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCatgatg                                                                                                                                              <  5:262114/70‑1 (MQ=255)
                                                                                      gTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGc                                                               >  2:95617/1‑149 (MQ=255)
                                                                                      gTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGc                                                               >  3:226260/1‑149 (MQ=255)
                                                                                       tCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGa                                                                                                      <  3:295710/109‑1 (MQ=255)
                                                                                       tCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGa                                                                                                      >  4:295710/1‑109 (MQ=255)
                                                                                         aGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGa                                                                                       <  5:298076/122‑1 (MQ=255)
                                                                                         aGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCttt                                                            <  5:295689/149‑1 (MQ=255)
                                                                                         aGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACACATATGCAAACTGCCGGTGACGTTCGGCGGCGGGa                                                                                       >  6:298076/1‑122 (MQ=255)
                                                                                          gTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTc                                                           <  5:341555/149‑1 (MQ=255)
                                                                                             ccGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATc                                                        <  1:147920/149‑1 (MQ=255)
                                                                                             ccGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATc                                                        <  1:95617/149‑1 (MQ=255)
                                                                                             ccGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATc                                                        >  2:373524/1‑149 (MQ=255)
                                                                                                  cTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCAt                                                   >  1:238740/1‑149 (MQ=255)
                                                                                                  cTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCAt                                                   <  3:286256/149‑1 (MQ=255)
                                                                                                  cTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCAATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGTAGAAACTTTCCTGGCTTTCATCGCCAt                                                   >  6:200814/1‑149 (MQ=255)
                                                                                                    gCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGGGACGCATATGCAAACTGCCGGGGCCGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGTCTTTCATCGCCAt                                                   >  3:38665/1‑147 (MQ=255)
                                                                                                     ccACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGTTCTTGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGTATCATGATCCAGAAACTTTCCTGGCTTTCATCGCCATTa                                                 <  2:7967/148‑1 (MQ=255)
                                                                                                     ccACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTa                                                 <  5:142927/148‑1 (MQ=255)
                                                                                                     ccACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTa                                                 <  7:121708/148‑1 (MQ=255)
                                                                                                      cacaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAg                                                <  1:274264/148‑1 (MQ=255)
                                                                                                      cacaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAg                                                >  2:274264/1‑148 (MQ=255)
                                                                                                       acaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGc                                              <  5:218155/149‑1 (MQ=255)
                                                                                                       acaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGc                                              >  1:373046/1‑149 (MQ=255)
                                                                                                       acaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGc                                              >  5:86587/1‑149 (MQ=255)
                                                                                                       acaGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGc                                              >  4:30767/1‑149 (MQ=255)
                                                                                                           ggTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTa                                           >  8:87081/1‑148 (MQ=255)
                                                                                                            gtgtTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTCTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCACCATTAGGCTTa                                           >  6:290348/1‑147 (MQ=255)
                                                                                                             tgtTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAg                                        <  7:81528/149‑1 (MQ=255)
                                                                                                                      ccTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGc                               >  6:65410/1‑149 (MQ=255)
                                                                                                                      ccTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGc                               >  4:85706/1‑149 (MQ=255)
                                                                                                                        tGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTc                                                                    >  2:435/1‑110 (MQ=255)
                                                                                                                        tGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTc                                                                    <  1:435/110‑1 (MQ=255)
                                                                                                                         gCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCttt                            >  7:213303/1‑149 (MQ=255)
                                                                                                                         gCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCttt                            >  5:29942/1‑149 (MQ=255)
                                                                                                                              cgcTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCa                       <  4:226260/149‑1 (MQ=255)
                                                                                                                              cgcTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCa                       <  7:87081/149‑1 (MQ=255)
                                                                                                                                 tGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTg                    >  5:59951/1‑149 (MQ=255)
                                                                                                                                 tGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTACGCTTAAAGTAGCCCTGCTTTTCCCAGTg                    >  1:209045/1‑149 (MQ=255)
                                                                                                                                      gCGGATCATGGTATCCATGATGGTCTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCACCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGt               <  5:290348/149‑1 (MQ=255)
                                                                                                                                      aaccaTCATGGATACCATGTTGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTTGc                                                               <  5:328475/97‑1 (MQ=255)
                                                                                                                                      aaccaTCATGGATACCATGATGGTTTGCTGGAAGGCGTGACACATATGCAAACAGCCGGTCACGAACGGCGGCGGGATCATGATACACAAACTAACCTGGc                                                               >  6:328475/5‑101 (MQ=255)
                                                                                                                                       accaTCATGGATACCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCt                                                                  <  2:189459/94‑1 (MQ=255)
                                                                                                                                       accaTCATGGATACCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCt                                                                  >  1:189459/4‑97 (MQ=255)
                                                                                                                                       accaTCATGGATACCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGt               >  4:113470/4‑148 (MQ=255)
                                                                                                                                           tCATGGTTTCCATGATGGTTTGCTGGATGGCGTGACCCATATGCTAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGTAAAGc          <  2:209045/149‑1 (MQ=255)
                                                                                                                                                   tCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGTAAAGCGGCTGtt   <  6:59951/148‑1 (MQ=255)
                                                                                                                                                   tCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGTAAAGCGGCTGTTc  >  8:202963/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTTCGCGGCCGTAGTCGTGACGGGTTTTACCTTCTTCTGCGGCAATCTTGCGCTCAACGACCATCTGGGTAGCGATCCCGGCGTGGTCAGTACCGACCTGCCACAGGGTGTTTTTGCCCTGCATGCGCTGATAGCGGATCATGGTATCCATGATGGTTTGCTGGAAGGCGTGACCCATATGCAAACTGCCGGTGACGTTCGGCGGCGGGATCATGATGCAGAAACTTTCCTGGCTTTCATCGCCATTAGGCTTAAAGTAGCCCTGCTTTTCCCAGTGCTCGTAAAGCGGCTGTTC  >  NC_000913/4483512‑4483807

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: