Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,433,830 C→T 9.7% V66M (GTG→ATG)  qorB ← NAD(P)H:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,433,8300CT 9.7% 145.1 / 4.3 48V66M (GTG→ATG) qorBNAD(P)H:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (17/26);  new base T (3/2);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGAAGTCTCCATCGTTTTGTCTTGTGTTTAAATACGCTAACACCTAAACTTACTTTTAGTAAGTACGTACAAAAAGGTAAGTATGAAA  >  NC_000913/4433552‑4434114
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                            
gATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGt                                                                                                                                                                                                                                                                                    <  2:206022/289‑1 (MQ=255)
               ttCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTCTAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATACCATTACGATGCCGCGGGGCACGTTGACCCACTTCGCTCGACGAGAGCAGCAGCAGTTCTTCCATTCCCTGAAGTCCCGATGTCCGTGCg                                                                                                                                                                                                                                                                     >  1:15580/1‑289 (MQ=255)
                      gctgcCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCa                                                                                                                                                                                                                                                               <  2:214236/288‑1 (MQ=255)
                      gctgcCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCa                                                                                                                                                                                                                                                               <  2:128392/288‑1 (MQ=255)
                         gcCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGaa                                                                                                                                                                                                                                                                                      <  2:349405/262‑1 (MQ=255)
                         gcCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGaa                                                                                                                                                                                                                                                                                      >  1:349405/1‑262 (MQ=255)
                         gcCGGGGCGCTGGCGAGGTAGTTTTCACGGTACCAGCCGTTGCGGAGCAGGGTGTATACGCTGCCAGAATCAGCCAGCATTTTCTCCGTCTCGCTGTCATCATCGGCCAGGCCGAGAGGGGCGTTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGTCAGCGCCCCTTGCAGCATTAATCACCTTAAGATGATGCGGGGCCCGTTGACCCACTTCGCTGGAAGAGATCAGCAGTCGTTTTTCCATTCCATGaattccatctataggataggcttcatcg                                                                                                                                                                                                                                                            >  1:517046/1‑262 (MQ=255)
                                cgcTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGt                                                                                                                                                                                                                                                     <  2:6606/288‑1 (MQ=255)
                                            ggTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCa                                                                                                                                                                                                                                         <  2:161344/287‑1 (MQ=255)
                                            aGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCa                                                                                                                                                                                                                                         <  2:413598/288‑1 (MQ=255)
                                            aGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCa                                                                                                                                                                                                                                         <  1:76170/288‑1 (MQ=255)
                                                ttCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTa                                                                                                                                                                                                                                    <  1:306701/289‑1 (MQ=255)
                                                          cAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTtgtg                                                                                                                                                                                                                          <  2:227275/289‑1 (MQ=255)
                                                            gCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTtgt                                                                                                                                                                                                                           >  2:352297/1‑286 (MQ=255)
                                                                 tGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAg                                                                                                                                                                                                                   >  2:391435/1‑289 (MQ=255)
                                                                 tGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGAATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAg                                                                                                                                                                                                                   >  2:49232/1‑289 (MQ=255)
                                                                             tgtAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCttt                                                                                                                                                                                                       >  2:370603/1‑289 (MQ=255)
                                                                                aaaCGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGCTGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGGAATGCCTTTTGCTGCCATGGCCTGGGATTTTg                                                                                                                                                                                                     >  1:488129/1‑288 (MQ=255)
                                                                                          tgAAGCAGCCAAGATTGTCTCCGGCTAGATGTGCTCATCTGCGAGGCCCAGCGGGGAGGTATCTGCATGTAGCAGGCGTGTATAAGCGTTAAATTTCACGCCAGCCGCCTTTGCGGCAGTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTGCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTacga                                                                                                                                                                                          <  1:384272/288‑1 (MQ=255)
                                                                                                                                     aGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAgg                                                                                                                                               >  2:330682/1‑289 (MQ=255)
                                                                                                                                       gCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGa                                                                                                                                              <  1:440470/288‑1 (MQ=255)
                                                                                                                                        ccGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGAt                                                                                                                                             >  1:422561/1‑288 (MQ=255)
                                                                                                                                         cGAGCGGGGAGGTAGCTGCATGGAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGAtt                                                                                                                                            <  1:469201/288‑1 (MQ=255)
                                                                                                                                              ggggAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATa                                                                                                                                       <  2:372422/288‑1 (MQ=255)
                                                                                                                                                      aTCTGAATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGa                                                                                                                              <  1:49232/289‑1 (MQ=255)
                                                                                                                                                      aTCGGCATGTAGCAGGCGGGTATAATGGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGa                                                                                                                              <  1:412512/289‑1 (MQ=255)
                                                                                                                                                               tAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTgg                                                                                                                     >  1:61512/1‑289 (MQ=255)
                                                                                                                                                                       tGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGc                                                                                                             <  1:352297/289‑1 (MQ=255)
                                                                                                                                                                        ggTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACGCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGc                                                                                                             <  1:213630/288‑1 (MQ=255)
                                                                                                                                                                          tataAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCAc                                                                                                           >  1:333571/1‑288 (MQ=255)
                                                                                                                                                                                   tAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAg                                                                                                  <  2:422561/288‑1 (MQ=255)
                                                                                                                                                                                   tAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGc                                                                                                 >  1:214946/1‑289 (MQ=255)
                                                                                                                                                                                            cGCCAGCCGGCTTTGCGGCATTAAGAACATTACGATGATGCGGTGCACGTTGACCCACTTCGCTGGAAGAGATGCACAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCACGGCTTGGGCTTTTGCCAGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCAACAGTAATAGCGATCATggg                                                                                        <  2:517046/289‑1 (MQ=255)
                                                                                                                                                                                               cAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCATTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGaa                                                                                      <  2:15580/288‑1 (MQ=255)
                                                                                                                                                                                                  ccgccTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAGTAGCGATCATGGGAAGtc                                                                                   >  2:49806/1‑288 (MQ=255)
                                                                                                                                                                                                          gCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGAAGTCTCCATCGtt                                                                          <  1:330682/289‑1 (MQ=255)
                                                                                                                                                                                                               aTTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGAAGTCTCCATCGTTTTGTc                                                                     <  2:214946/289‑1 (MQ=255)
                                                                                                                                                                                                                ttAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGAAGTCTCCATCGTTTTGTCt                                                                    >  2:531550/1‑289 (MQ=255)
                                                                                                                                                                                                                               tgctgcGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCATTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTtggctggc                                                                                                                                                                   <  2:68024/179‑1 (MQ=255)
                                                                                                                                                                                                                               tgctgcGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCATTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTtggctggc                                                                                                                                                                   >  1:68024/1‑179 (MQ=255)
                                                                                                                                                                                                                                  tgcGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCa                                                                                                                                                                                                                                         <  1:51689/106‑1 (MQ=255)
                                                                                                                                                                                                                                  tgcGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCa                                                                                                                                                                                                                                         >  2:51689/1‑106 (MQ=255)
                                                                                                                                                                                                                                              tGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTacgaacga                                                                                                                                                                                      >  2:309676/1‑145 (MQ=255)
                                                                                                                                                                                                                                              tGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTacgaacga                                                                                                                                                                                      <  1:309676/145‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGAAGTCTCCATCGTTTTGTCTTGTGTTTAAATACGCTAACACCTAAact                                        <  2:255445/279‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGAAGTCTCCATCGTTTTGTCTTGTGTTTAAATACGCTAACACCTAAact                                        >  1:255445/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                   aTCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGAAGTCTCCATCGTTTTGTCTTGTGTTTAAATACGCTAACACCTAAACTTCCTTTTAGTAAGTACGTACaaa                 >  1:20852/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                  tCCACGCCCCGAGGTGCGCATGTCAGTGCGGCTTCACCGCCGTAGTCAGCCTGACGAACGGTAATGCCTTGGGCTGCCAGGGCTTGGACTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGAAGTCTCCATCGTTTTGTCTTGTGTTTAAATACGCTAACACCTAAACTTACTTTTAGTAAGTACGTACAAAAAGGTAAGTATGaaa  <  1:370603/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                            
GATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAGCATTTTCTCCGTCTCGATGTGCTCATCGGCGAGGCCGAGCGGGGAGGTATCTGCATGTAGCAGGCTGGTATAAGCGATAAATTTCACGCCAGCCGCCTTTGCGGCATTAATAACATTACGATGCTGCGGGGCACGTTGACCCACTTCGCTGGAAGAGATCAGCAGTAGTTTTTCCACTCCCTGAAGTGCAGATGTCAGTGCGGCTTCATCGCCGTAGTCAGCCTGACGCACGGTAATGCCTTGTGCTGCCAGGGCTTGGGCTTTTGCCGGATTACGAACGATAGCCACTATTTGGCTGGCAGGAACCGTTTTCATCAAGGATTCAATAACATAGTGACCAAGTTGGCCAGTGGCACCAGTAATAGCGATCATGGGAAGTCTCCATCGTTTTGTCTTGTGTTTAAATACGCTAACACCTAAACTTACTTTTAGTAAGTACGTACAAAAAGGTAAGTATGAAA  >  NC_000913/4433552‑4434114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑