Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,951,475 T→A 5.8% L329L (CTA→CTT 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,4750TA 5.8% 360.6 / 9.9 104L329L (CTA→CTTrecDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base T (45/53);  new base A (2/4);  total (47/57)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

Only 100 of 105 total aligned reads displayed.

ATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGTCTC  >  NC_000913/2951328‑2951622
                                                                                                                                                   |                                                                                                                                                   
aTCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTa                                                                                                                                                    >  3:899695/1‑149 (MQ=255)
  caccacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGc                                                                                                                                                  <  1:558469/149‑1 (MQ=255)
    ccagGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTGATGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTg                                                                                                                                                >  4:788382/1‑149 (MQ=255)
     cacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGc                                                                                                                                               >  2:583926/1‑149 (MQ=255)
     cacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGc                                                                                                                                               >  3:108932/1‑149 (MQ=255)
     cacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGc                                                                                                                                               <  3:237652/149‑1 (MQ=255)
     cacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGc                                                                                                                                               >  4:553562/1‑149 (MQ=255)
      acGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGCCAGGCAGCTAAgccgcc                                                                                                                                              <  3:504345/149‑1 (MQ=255)
        ggTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTg                                                                                                                                            <  1:415146/149‑1 (MQ=255)
         gTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTgg                                                                                                                                           <  1:631176/149‑1 (MQ=255)
             aTCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGc                                                                                                                                               <  4:377847/141‑1 (MQ=255)
             aTCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGc                                                                                                                                               >  3:377847/1‑141 (MQ=255)
             aTCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACg                                                                                                                                       >  2:332023/1‑149 (MQ=255)
             aTCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACg                                                                                                                                       >  1:574703/1‑149 (MQ=255)
              tCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGCCAGGCAGCTAAGCCGCCTgacgggg                                                                                                                                      <  2:293101/149‑7 (MQ=255)
               cgccgcAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCAGCTAAGCCGCCTgacgggga                                                                                                                                     <  4:474829/149‑8 (MQ=255)
                   gcAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggc                                                                                                                                 <  2:294879/149‑1 (MQ=255)
                   gcAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggc                                                                                                                                 <  2:210361/149‑1 (MQ=255)
                    cAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggcg                                                                                                                                <  4:164939/149‑1 (MQ=255)
                     aGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggcgg                                                                                                                               >  4:150961/1‑149 (MQ=255)
                      gCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGt                                                                                                                              >  2:229525/1‑149 (MQ=255)
                        tAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTaa                                                                                                                            >  1:628590/1‑149 (MQ=255)
                         aaCTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTaaa                                                                                                                           <  1:332023/149‑1 (MQ=255)
                         aaCTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTaaa                                                                                                                           <  3:553562/149‑1 (MQ=255)
                         aaCTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTaaa                                                                                                                           <  1:866054/149‑1 (MQ=255)
                          aCTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggcgg                                                                                                                               >  2:153523/1‑144 (MQ=255)
                          aCTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggcgg                                                                                                                               <  1:153523/144‑1 (MQ=255)
                           cTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGc                                                                                                                         >  1:210565/1‑149 (MQ=255)
                           cTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGc                                                                                                                         <  4:67332/149‑1 (MQ=255)
                                cAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCcgcg                                                                                                                    >  4:609765/1‑149 (MQ=255)
                                cAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCcgcg                                                                                                                    >  2:39602/1‑149 (MQ=255)
                                 aaTGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGt                                                                                                                   >  1:494330/1‑149 (MQ=255)
                                  aTGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGtt                                                                                                                  <  3:733810/149‑1 (MQ=255)
                                  aTGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCCGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGtt                                                                                                                  <  2:930883/149‑1 (MQ=255)
                                   tGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTg                                                                                                                 >  3:686059/1‑149 (MQ=255)
                                   tGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTg                                                                                                                 <  2:415229/149‑1 (MQ=255)
                                     ccAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGc                                                                                                               <  3:150961/149‑1 (MQ=255)
                                         aaTCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATaa                                                                                                           <  1:297356/149‑1 (MQ=255)
                                         aaTCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATaa                                                                                                           <  1:173885/149‑1 (MQ=255)
                                                 ccGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGat                                                                                                   >  1:299604/1‑149 (MQ=255)
                                                   gAAACGATAGCTTTTTTGCTGCAGGCAGAGACTGTCGCGCAAAGATGCAGCTTCTGTGCCGATTCCTGCCGGAACGTGAGTCCCCGTCCGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGGATAAGCGCAGatat                                                                                                 >  4:881619/1‑149 (MQ=255)
                                                      aCGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTTCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGc                                                                                              >  4:679085/1‑149 (MQ=255)
                                                       cGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGCCATAAGCGCAGATATCGcc                                                                                             <  1:288586/149‑1 (MQ=255)
                                                              ttttttGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGcaca                                                                                      >  1:688599/1‑149 (MQ=255)
                                                              ttttttGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGcaca                                                                                      >  2:559969/1‑149 (MQ=255)
                                                              ttttttGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGcaca                                                                                      >  2:503336/1‑149 (MQ=255)
                                                                   tGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGcccc                                                                                 <  1:456048/149‑1 (MQ=255)
                                                                      agcagGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGc                                                                              <  4:681082/149‑1 (MQ=255)
                                                                      agcagGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGc                                                                              <  1:112139/149‑1 (MQ=255)
                                                                      agcagGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGc                                                                              <  3:161712/149‑1 (MQ=255)
                                                                         aggcagAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTc                                                                           >  3:791469/1‑149 (MQ=255)
                                                                             agagACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAAcc                                                                       >  3:376694/1‑149 (MQ=255)
                                                                               agaCTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTAAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGa                                                                     <  1:229525/149‑1 (MQ=255)
                                                                                gaCTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAg                                                                    <  2:75215/149‑1 (MQ=255)
                                                                                gaCTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAg                                                                    <  2:915706/149‑1 (MQ=255)
                                                                                       gcgcAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGtt                                                             <  4:767184/149‑1 (MQ=255)
                                                                                         gcAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGa                                                           >  2:548587/1‑149 (MQ=255)
                                                                                         gcAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGa                                                           <  2:688599/149‑1 (MQ=255)
                                                                                          cAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGAt                                                          <  2:726404/149‑1 (MQ=255)
                                                                                          cAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGAt                                                          <  4:108932/149‑1 (MQ=255)
                                                                                          cAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGAt                                                          <  4:686059/149‑1 (MQ=255)
                                                                                           aaaGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATc                                                         <  1:559969/149‑1 (MQ=255)
                                                                                            aaGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCa                                                        <  2:496046/149‑1 (MQ=255)
                                                                                                tGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGAt                                                    <  2:494330/149‑1 (MQ=255)
                                                                                                   cGCTTCTGTGCCGGTTCCTGCCGGCACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGc                                                 >  4:850310/1‑149 (MQ=255)
                                                                                                   cGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGc                                                 >  4:691646/1‑149 (MQ=255)
                                                                                                      ttCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCga                                              >  4:263477/1‑149 (MQ=255)
                                                                                                          gtgCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGTGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGaaa                                          <  4:577891/149‑1 (MQ=255)
                                                                                                           tgCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAg                                         <  3:51787/149‑1 (MQ=255)
                                                                                                           ctcacgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAg                                         >  2:56367/6‑149 (MQ=255)
                                                                                                             ccGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGAt                                       >  2:626830/1‑149 (MQ=255)
                                                                                                              acgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATc                                      >  1:144602/3‑149 (MQ=255)
                                                                                                              acgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATc                                      >  3:83747/3‑149 (MQ=255)
                                                                                                               ggTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCa                                     >  4:730855/1‑149 (MQ=255)
                                                                                                               ggTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCa                                     >  4:139936/1‑149 (MQ=255)
                                                                                                               ggTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCa                                     <  4:82317/149‑1 (MQ=255)
                                                                                                                gTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCAc                                    >  2:131429/1‑149 (MQ=255)
                                                                                                                 ttCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGcc                                                                             >  3:189228/1‑107 (MQ=255)
                                                                                                                 ttCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGcc                                                                             <  4:189228/107‑1 (MQ=255)
                                                                                                                     tGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTcgcg                               >  2:280234/1‑149 (MQ=255)
                                                                                                                               tGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCggg                     >  2:594969/1‑149 (MQ=255)
                                                                                                                               tGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCggg                     <  1:768392/149‑1 (MQ=255)
                                                                                                                                  gTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCaa                  <  3:251288/149‑1 (MQ=255)
                                                                                                                                   tCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAg                 <  1:62236/149‑1 (MQ=255)
                                                                                                                                   tCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAg                 >  1:882316/1‑149 (MQ=255)
                                                                                                                                    ccccGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAgg                >  3:758882/1‑149 (MQ=255)
                                                                                                                                    ccccGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAgg                <  3:691646/149‑1 (MQ=255)
                                                                                                                                     cccGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGc               >  4:144132/1‑149 (MQ=255)
                                                                                                                                        gccAGGCAGCTAAGCCGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTc            >  3:598886/3‑149 (MQ=255)
                                                                                                                                        gTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTc            <  2:882316/149‑1 (MQ=255)
                                                                                                                                        gTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTc            <  4:134223/149‑1 (MQ=255)
                                                                                                                                        gTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTc            <  1:280234/149‑1 (MQ=255)
                                                                                                                                          cAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGa          <  2:578333/149‑1 (MQ=255)
                                                                                                                                          cAGGCAGCTAAGCCGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGa          >  1:873987/1‑149 (MQ=255)
                                                                                                                                           aGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGAt         >  4:776679/1‑149 (MQ=255)
                                                                                                                                           aGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGAt         <  1:599683/149‑1 (MQ=255)
                                                                                                                                             gcggcTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCa       <  3:609765/149‑1 (MQ=255)
                                                                                                                                                 cTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGtct   <  3:730855/149‑1 (MQ=255)
                                                                                                                                                 cTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGtct   <  3:139936/149‑1 (MQ=255)
                                                                                                                                                  ttAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGtctc  >  1:118763/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGCATGATCGGGCAAGGCGTCGATCAGTCTC  >  NC_000913/2951328‑2951622

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