Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,951,479 T→C 9.0% Q328R (CAG→CGG)  recD ← exonuclease V (RecBCD complex), alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,951,4790TC 9.0% 332.7 / 18.6 111Q328R (CAG→CGG) recDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base T (46/55);  new base C (3/7);  total (49/62)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 115 total aligned reads displayed.

ACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGTCTCGACAT  >  NC_000913/2951331‑2951627
                                                                                                                                                    |                                                                                                                                                    
accacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCt                                                                                                                                                      >  4:237790/1‑149 (MQ=255)
accacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCt                                                                                                                                                      >  2:457560/1‑149 (MQ=255)
 ccacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTg                                                                                                                                                     >  1:322741/1‑149 (MQ=255)
  cacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGc                                                                                                                                                    >  2:845798/1‑149 (MQ=255)
   acGGTTTATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTATTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACAGAAGTCCCCGTCAGGCGGCTTAGCTGcc                                                                                                                                                   >  2:846076/1‑149 (MQ=255)
   acGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGcc                                                                                                                                                   <  3:629834/149‑1 (MQ=255)
    cGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCt                                                                                                                                                  <  3:857268/149‑1 (MQ=255)
      gTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTgg                                                                                                                                                >  4:786362/1‑149 (MQ=255)
      gTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCAGCTAAGCCGCCTga                                                                                                                                                <  3:18714/149‑2 (MQ=255)
        tGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGCCAGGCAGCTAAGCCGCCTgacg                                                                                                                                              <  1:520853/149‑4 (MQ=255)
         gATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCGGAAACGATAGCTTTGTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCAGCTAAGCCGCCTgacgg                                                                                                                                             <  1:709416/149‑5 (MQ=255)
           tCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGc                                                                                                                                           >  2:848502/1‑149 (MQ=255)
              ccgcAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGCCAGGCAGCTAAGCCGCCTgacggg                                                                                                                                            <  1:346272/145‑6 (MQ=255)
              ccgcAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGCCAGGCAGCTAAGCCGCCTgacggg                                                                                                                                            >  2:346272/1‑140 (MQ=255)
                 cAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggcg                                                                                                                                     <  1:454353/149‑1 (MQ=255)
                  aGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggcgg                                                                                                                                    <  3:7348/149‑1 (MQ=255)
                  aGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggcgg                                                                                                                                    <  1:482219/149‑1 (MQ=255)
                     tAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTaa                                                                                                                                 >  1:580627/1‑149 (MQ=255)
                      aaCTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTaaa                                                                                                                                >  1:284416/1‑149 (MQ=255)
                      aaCTGACCAATGCCAGAATCGCTGCCGAAAAGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAACATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGGaaa                                                                                                                                >  1:481643/1‑149 (MQ=255)
                        cTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGc                                                                                                                              <  1:303312/149‑1 (MQ=255)
                         tGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGcc                                                                                                                             >  4:174451/1‑149 (MQ=255)
                             cAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGCCAGGCAGCTAAGCCGCCTgacgg                                                                                                                                             <  2:669491/129‑5 (MQ=255)
                             cAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGCCAGGCAGCTAAGCCGCCTgacgg                                                                                                                                             <  2:48063/129‑5 (MQ=255)
                             cAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGCCAGGCAGCTAAGCCGCCTgacgg                                                                                                                                             >  1:48063/1‑125 (MQ=255)
                                 gCCAGAATCGCGGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTgg                                                                                                                     <  2:59446/149‑1 (MQ=255)
                                   cAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGGTTCTGTGCCGGTTTCTTCGGGTACGGGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGTGGTAAAGCCCGCGTTGGCa                                                                                                                   >  2:638942/1‑149 (MQ=255)
                                      aaTCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATaa                                                                                                                <  1:470775/149‑1 (MQ=255)
                                      aaTCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATaa                                                                                                                <  4:302526/149‑1 (MQ=255)
                                      aaTCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATaa                                                                                                                <  1:7428/149‑1 (MQ=255)
                                        tCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAgc                                                                                                              <  1:845798/149‑1 (MQ=255)
                                         cGCTACCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAgcg                                                                                                             >  2:150309/1‑149 (MQ=255)
                                                gAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggcg                                                                                                                                     <  1:596059/118‑1 (MQ=255)
                                                gAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGatat                                                                                                      >  1:608376/1‑149 (MQ=255)
                                                  aaCGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGTATAAGCGCAGATATCg                                                                                                    <  1:189209/149‑1 (MQ=255)
                                                     gATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGccc                                                                                                 <  2:41128/149‑1 (MQ=255)
                                                     gATAGCTGTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGccc                                                                                                 >  2:477041/1‑149 (MQ=255)
                                                       tAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAg                                                                                               >  2:849454/1‑149 (MQ=255)
                                                       tAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAg                                                                                               >  1:281936/1‑149 (MQ=255)
                                                             ttttGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCATCACAGc                                                                                         <  2:478807/149‑1 (MQ=255)
                                                              tttGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGcc                                                                                        <  3:204536/149‑1 (MQ=255)
                                                                  cagcagGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAg                                                                                    >  3:684999/1‑149 (MQ=255)
                                                                   agcagGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGc                                                                                   <  4:829846/149‑1 (MQ=255)
                                                                    gcaggcagAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGcc                                                                                  >  1:521518/1‑149 (MQ=255)
                                                                     caggcagAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCt                                                                                 <  2:749890/149‑1 (MQ=255)
                                                                          agagACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAAcc                                                                            >  2:686868/1‑149 (MQ=255)
                                                                             gaCTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAg                                                                         <  2:322741/149‑1 (MQ=255)
                                                                                tGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGcc                                                                      >  3:844572/1‑149 (MQ=255)
                                                                                      gcAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGa                                                                <  3:402100/149‑1 (MQ=255)
                                                                                        aaaGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATc                                                              >  1:208789/1‑149 (MQ=255)
                                                                                         aaGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCa                                                             >  1:197408/1‑149 (MQ=255)
                                                                                            aTGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGa                                                          <  1:711745/149‑1 (MQ=255)
                                                                                                 gCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGcc                                                     >  2:880832/1‑149 (MQ=255)
                                                                                                 gCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGcc                                                     >  3:242051/1‑149 (MQ=255)
                                                                                                   ttCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCga                                                   >  2:56672/1‑149 (MQ=255)
                                                                                                   ttCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCga                                                   >  3:156185/1‑149 (MQ=255)
                                                                                                   ttCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCga                                                   >  3:798179/1‑149 (MQ=255)
                                                                                                     cTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCgaga                                                 <  3:639645/149‑1 (MQ=255)
                                                                                                       gtgCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGaaa                                               <  1:630202/149‑1 (MQ=255)
                                                                                                       gtgCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGaaa                                               <  4:810204/149‑1 (MQ=255)
                                                                                                       gtgCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGaaa                                               <  1:364682/149‑1 (MQ=255)
                                                                                                       gtgCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGaaa                                               <  3:288023/149‑1 (MQ=255)
                                                                                                       gtgCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGaaa                                               <  2:197408/149‑1 (MQ=255)
                                                                                                        tgCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAg                                              <  2:521518/149‑1 (MQ=255)
                                                                                                         tcacgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCa                                                             <  4:741999/129‑1 (MQ=255)
                                                                                                         tcacgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCa                                                             >  3:741999/5‑133 (MQ=255)
                                                                                                         tcacgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGa                                             >  4:591589/5‑149 (MQ=255)
                                                                                                         tcacgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGa                                             >  2:855485/5‑149 (MQ=255)
                                                                                                         tcacgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGa                                             >  2:227867/5‑149 (MQ=255)
                                                                                                          ccGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGAt                                            <  2:580627/149‑1 (MQ=255)
                                                                                                           cGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATc                                           <  3:688953/149‑1 (MQ=255)
                                                                                                            ggTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCa                                          <  4:502268/149‑1 (MQ=255)
                                                                                                            ggTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCa                                          >  1:164066/1‑149 (MQ=255)
                                                                                                                 cTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTcgc                                     >  1:179921/1‑149 (MQ=255)
                                                                                                                 cTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTcgc                                     >  2:287897/1‑149 (MQ=255)
                                                                                                                  tGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTcgcg                                    >  1:122728/1‑149 (MQ=255)
                                                                                                                    ccGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCa                                  >  3:276021/1‑149 (MQ=255)
                                                                                                                    ccGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCa                                  >  2:508776/1‑149 (MQ=255)
                                                                                                                      ggAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATg                                >  2:136455/1‑149 (MQ=255)
                                                                                                                      ggAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGACCGCGTTGGCTTAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGGTCACGATCGTCGAGAAAGATCACTCGCGCATg                                >  2:50485/1‑149 (MQ=255)
                                                                                                                        aaCGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGAt                              <  4:798179/149‑1 (MQ=255)
                                                                                                                         aCGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATc                             <  3:129637/149‑1 (MQ=255)
                                                                                                                           gTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCgg                           >  3:805913/1‑149 (MQ=255)
                                                                                                                           gTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGACGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCgg                           <  4:842896/149‑1 (MQ=255)
                                                                                                                             gAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGc                         >  3:585412/1‑149 (MQ=255)
                                                                                                                                tCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAg                      <  1:428177/149‑1 (MQ=255)
                                                                                                                                  cccGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGc                    <  1:880832/149‑1 (MQ=255)
                                                                                                                                       cAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGa               <  3:82355/149‑1 (MQ=255)
                                                                                                                                       cAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGa               <  2:164338/149‑1 (MQ=255)
                                                                                                                                       cAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGa               <  1:227867/149‑1 (MQ=255)
                                                                                                                                       cAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGa               <  3:237790/149‑1 (MQ=255)
                                                                                                                                         ggcggcTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATc             <  2:608376/149‑1 (MQ=255)
                                                                                                                                          gcggcTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCa            >  1:728620/1‑149 (MQ=255)
                                                                                                                                          gcggcTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCa            <  1:715699/149‑1 (MQ=255)
                                                                                                                                          gcggcTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGAACACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCa            >  2:715699/1‑149 (MQ=255)
                                                                                                                                            ggcTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGt          <  4:156864/149‑1 (MQ=255)
                                                                                                                                              cTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGtct        <  1:533935/149‑1 (MQ=255)
                                                                                                                                                 aGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGTCTCGa     <  4:6667/149‑1 (MQ=255)
                                                                                                                                                    tGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGTCTCGAcat  <  2:164066/149‑1 (MQ=255)
                                                                                                                                                    tGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGTCTCGAcat  <  1:477041/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGTCTCGACAT  >  NC_000913/2951331‑2951627

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