Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,020,252 A→T 5.4% E364D (GAA→GAT ssnA → putative chlorohydrolase/aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,2520AT 5.4% 335.1 / 3.9 104E364D (GAA→GATssnAputative chlorohydrolase/aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (51/47);  new base T (3/3);  total (54/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

TTCGGACATGTTTGAAGAGATGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGA  >  NC_000913/3020105‑3020398
                                                                                                                                                   |                                                                                                                                                  
ttCGGACATGTTTGAAGAGATGAAATTTGCCTGCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCACAGCCCTGACTAACGGTAACGACCTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTGTGGAAg                                                                                                                                                   >  6:114266/1‑149 (MQ=255)
   ggACATGTTTGAAGAGATGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCg                                                                                                                                                <  3:320776/149‑1 (MQ=255)
      cATGTTTGAAGAGATGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGc                                                                                                                                                  <  4:115050/144‑1 (MQ=255)
      cATGTTTGAAGAGATGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGc                                                                                                                                                  >  3:115050/1‑144 (MQ=255)
              aaGAGATGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCt                                                                                                                                     <  6:235968/149‑1 (MQ=255)
              aaGAGATGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCt                                                                                                                                     >  2:53275/1‑149 (MQ=255)
              aaGAGATGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCt                                                                                                                                     >  3:85229/1‑149 (MQ=255)
                gagaTGAAATTTGCCTTCTTTAAACGTCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGa                                                                                                                                   >  2:214327/1‑149 (MQ=255)
                   aTGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGAtt                                                                                                                                 >  2:4455/1‑148 (MQ=255)
                    tGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGAttt                                                                                                                                >  5:161027/1‑148 (MQ=255)
                    tGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGAttt                                                                                                                                >  5:114175/1‑148 (MQ=255)
                       aaTTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAAc                                                                                                                             >  3:73008/1‑148 (MQ=255)
                        aTTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCa                                                                                                                           >  6:126450/1‑149 (MQ=255)
                         tttGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCAt                                                                                                                          <  7:325663/149‑1 (MQ=255)
                          ttGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCAtt                                                                                                                         >  4:365502/1‑149 (MQ=255)
                             ccTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTg                                                                                                                       <  4:168941/148‑1 (MQ=255)
                               ttcttTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGATGCCGGTTACAAAGCTGAt                                                                                                                                  >  3:198794/1‑135 (MQ=255)
                               ttcttTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGATGCCGGTTACAAAGCTGAt                                                                                                                                  <  4:198794/135‑1 (MQ=255)
                               ttcttTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGATGCCGGTTACAAAGCTGAt                                                                                                                                  <  4:204377/135‑1 (MQ=255)
                               ttcttTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGATGCCGGTTACAAAGCTGAt                                                                                                                                  >  3:204377/1‑135 (MQ=255)
                                 cttTAAACATCGCGATGCTGGTGGTCCGCTGTGGCCTCACAGATTTGCCAAAGCCCTGACTAACGGAAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGAtat                                                                                                                                >  1:226671/1‑133 (MQ=255)
                                  tttAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTa                                                                                                                 >  6:29892/1‑149 (MQ=255)
                                    tAAACCTCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTAc                                                                                                                <  2:250803/148‑1 (MQ=255)
                                    tAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTa                                                                                                                 <  1:53275/147‑1 (MQ=255)
                                    tAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTAc                                                                                                                >  5:194655/1‑148 (MQ=255)
                                        cATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTc                                                                                                           >  6:183655/1‑149 (MQ=255)
                                         aTCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCg                                                                                                          >  4:89439/1‑149 (MQ=255)
                                            gcgATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGcc                                                                                                        >  5:107052/1‑148 (MQ=255)
                                                    ggtggtCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGtt                                                                                                                                             <  2:270157/103‑1 (MQ=255)
                                                    ggtggtCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGtt                                                                                                                                             >  1:270157/1‑103 (MQ=255)
                                                           cGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTgg                                                                                          >  3:190950/1‑147 (MQ=255)
                                                           cGCTGGGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGACCTGATGAGCCGCACCTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGc                                                                                         <  3:347461/148‑1 (MQ=255)
                                                             cTGTGGCCTGACAGGTTTGCCAAAGGCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTTGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGa                                                                                      <  6:237464/149‑1 (MQ=255)
                                                               gtgGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACa                                                                                    <  5:197089/149‑1 (MQ=255)
                                                                tgGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGAc                                                                                     <  5:126450/147‑1 (MQ=255)
                                                                tgGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACaa                                                                                   >  8:49407/1‑149 (MQ=255)
                                                                tgGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACa                                                                                    <  1:214327/148‑1 (MQ=255)
                                                                  gCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAata                                                                                 <  1:423568/149‑1 (MQ=255)
                                                                  gCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACCGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGCCAat                                                                                  >  2:271441/1‑148 (MQ=255)
                                                                    cTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCTTTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATc                                                                               >  2:313025/1‑149 (MQ=255)
                                                                    cTGACAGTTTTGCAAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATc                                                                               <  4:85229/149‑1 (MQ=255)
                                                                         aGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGCAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCgg                                                                          >  2:5214/1‑149 (MQ=255)
                                                                           ttttGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTTGCGCGAAATTTGGTCTTCTGGAATCCGTTTACAAAGCTGATTTaa                                                                                                                              >  8:170441/1‑95 (MQ=255)
                                                                           ttttGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCg                                                                                                          <  8:227790/115‑1 (MQ=255)
                                                                           ttttGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCg                                                                                                          >  7:227790/1‑115 (MQ=255)
                                                                           ttttGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCggg                                                                         >  1:99917/1‑148 (MQ=255)
                                                                           ttttGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGc                                                                        <  3:365502/149‑1 (MQ=255)
                                                                              tGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCat                                                                      <  5:352929/148‑1 (MQ=255)
                                                                                ccAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCatat                                                                    <  5:183655/148‑1 (MQ=255)
                                                                                 cAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCg                                                                  >  3:89272/1‑149 (MQ=255)
                                                                                  aaaGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGc                                                                 >  5:414924/1‑149 (MQ=255)
                                                                                    aGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCtt                                                               >  8:372266/1‑149 (MQ=255)
                                                                                       ccTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCg                                                            <  8:276212/149‑1 (MQ=255)
                                                                                           aCTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTCCAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATAGCGCTTTCGGTAt                                                        <  7:182855/149‑1 (MQ=255)
                                                                                           aCTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGt                                                          <  8:218405/147‑1 (MQ=255)
                                                                                           aCTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTa                                                         <  7:372266/148‑1 (MQ=255)
                                                                                           aCTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTa                                                         <  6:161027/148‑1 (MQ=255)
                                                                                           aCTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTa                                                         <  6:114175/148‑1 (MQ=255)
                                                                                           aCTAACGGTAACGAACTGAAGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTa                                                         <  1:4455/148‑1 (MQ=255)
                                                                                              aaCGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATggg                                                     >  2:143968/1‑149 (MQ=255)
                                                                                                 ggaaaCGAACTGATGAGCCGCAACTTTGGCGCGAAAATTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTAAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATAACGCCGGGAATATCGCTTTCg                                                            >  6:273734/4‑139 (MQ=255)
                                                                                                 ggTAAGGATCTGATGAGCCGCAAGTTTGGCGCGAAATTTGGGGTTCTGGATGCCGGTTACAAAGCTGATTTTTCCATTTGCGATTACAACTCGCCGACGCCGGTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCg                                                            <  5:273734/139‑1 (MQ=255)
                                                                                                 ggTAACGAACTGATGAGCCTCAACTTTGGCGCGAAATTTGGGCTTCTGGATGCCCGTTACAAAGCTGATTTAACCATTTGCGATTAAAACTCGCCGACGCCGCTGCTTGCAGACAATATAGCCGGGCATATCGCTTTCGGTATGGGCt                                                   >  3:137196/1‑148 (MQ=255)
                                                                                                   tAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGt                                                          >  1:281093/1‑139 (MQ=255)
                                                                                                   tAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGt                                                          <  2:281093/139‑1 (MQ=255)
                                                                                                   tAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAg                                                >  5:239923/1‑149 (MQ=255)
                                                                                                     aCGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAgg                                               <  4:190950/148‑1 (MQ=255)
                                                                                                       gAACTGATGAGACGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAg                                            <  6:194655/149‑1 (MQ=255)
                                                                                                        aaCTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGc                                           <  3:89439/149‑1 (MQ=255)
                                                                                                           tgatgaGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACATCTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGtt                                        >  8:295277/1‑149 (MQ=255)
                                                                                                           tgatgaGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGtt                                        >  2:60462/1‑149 (MQ=255)
                                                                                                           tgatgaGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGtt                                        >  6:360435/1‑149 (MQ=255)
                                                                                                            gatgaGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGACGCTGCTGGCAGACAATATCGCCGGGCATATCGCTGTCGGTATGGGCTCCGTCAGCGTTc                                       >  5:312172/1‑149 (MQ=255)
                                                                                                                aGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATTTCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAg                                   <  5:29892/149‑1 (MQ=255)
                                                                                                                   cGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGt                                <  6:239923/149‑1 (MQ=255)
                                                                                                                   cGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGt                                >  3:86685/1‑149 (MQ=255)
                                                                                                                   cGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGt                                >  2:62307/1‑149 (MQ=255)
                                                                                                                    gCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGt                                >  6:376567/1‑148 (MQ=255)
                                                                                                                    gCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGt                                >  2:103351/1‑148 (MQ=255)
                                                                                                                    gCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTg                               <  5:182178/149‑1 (MQ=255)
                                                                                                                       aCTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATg                            <  6:414924/149‑1 (MQ=255)
                                                                                                                           tGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGt                          <  1:143968/147‑1 (MQ=255)
                                                                                                                           tGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGt                          <  7:83090/147‑1 (MQ=255)
                                                                                                                             gcgcGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCaa                       <  2:99917/148‑1 (MQ=255)
                                                                                                                               gcgAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATgg                    >  6:147974/1‑149 (MQ=255)
                                                                                                                                cgAAATTTGGGCTTCTGGAATCCGGTTACAAAGCTGATTTAACAATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGgt                   <  5:331067/149‑1 (MQ=255)
                                                                                                                                cgAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGgt                   <  1:103351/149‑1 (MQ=255)
                                                                                                                                 gAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGgtg                  <  4:86685/149‑1 (MQ=255)
                                                                                                                                  aaaTTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCg                                                                  <  1:131517/100‑1 (MQ=255)
                                                                                                                                  aaaTTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCg                                                                  >  2:131517/1‑100 (MQ=255)
                                                                                                                                       tGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCTTTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATg             <  1:313025/148‑1 (MQ=255)
                                                                                                                                       ccggCTTCCAGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATg                            <  7:37747/131‑1 (MQ=255)
                                                                                                                                       ccggCTTCCAGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATg                            >  8:37747/3‑133 (MQ=255)
                                                                                                                                            ttCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATg             >  6:298997/1‑143 (MQ=255)
                                                                                                                                            ttCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATg             <  5:298997/143‑1 (MQ=255)
                                                                                                                                            ttCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTAt       >  2:394714/1‑149 (MQ=255)
                                                                                                                                              cTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATTGTGTGATGGTCTATg      <  5:376567/148‑1 (MQ=255)
                                                                                                                                                 gAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCCCAGCGTGATGGTCAATGGTGTGATGGTCTATGAAgc  >  1:292509/1‑148 (MQ=255)
                                                                                                                                                 gAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGaa    <  2:354091/147‑1 (MQ=255)
                                                                                                                                                 gAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGa  >  1:298308/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCGGACATGTTTGAAGAGATGAAATTTGCCTTCTTTAAACATCGCGATGCGGGTGGTCCGCTGTGGCCTGACAGTTTTGCCAAAGCCCTGACTAACGGTAACGAACTGATGAGCCGCAACTTTGGCGCGAAATTTGGGCTTCTGGAAGCCGGTTACAAAGCTGATTTAACCATTTGCGATTACAACTCGCCGACGCCGCTGCTGGCAGACAATATCGCCGGGCATATCGCTTTCGGTATGGGCTCAGGCAGCGTTCACAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGA  >  NC_000913/3020105‑3020398

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