Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,960,928 C→G 6.0% G501R (GGC→CGC)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,9280CG 6.0% 215.0 / 5.0 78G501R (GGC→CGC) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base G (2/3);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

TAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960780‑2961076
                                                                                                                                                    |                                                                                                                                                    
tAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGcc                                                                                                                                                      <  1:320698/149‑1 (MQ=255)
tAAGGTAGAACCGATGGCCTCGCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTGGTCGTCATCTATGCCCCAACGAATGcc                                                                                                                                                      <  1:848257/149‑1 (MQ=255)
 acgggaGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGGTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCg                                                                                                                                                     <  1:486647/144‑1 (MQ=255)
 aaGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCCTCGCGTAGCCCCACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCg                                                                                                                                                     <  2:244788/149‑1 (MQ=255)
  aGGTAGAACCGATAGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAAAATACGCGACAGGCCAAAACGCCAGGAGTGTTCTCCGGTGGCGGGGAGTTCCATCTCGCGAACGAAGTCGTGAGCTATGCCCCAACCAATGCggg                                                                                                                                                    >  2:702735/1‑149 (MQ=255)
    gTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATAGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTAGCGATCGTTTTCGTCATCTATGCCCCAACGAATGCCggct                                                                                                                                                  >  1:581785/1‑147 (MQ=255)
     tAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAtt                                                                                                                                                 <  1:836295/149‑1 (MQ=255)
        aaCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGt                                                                                                                                              >  1:474046/1‑149 (MQ=255)
        aaCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGt                                                                                                                                              <  2:2255/149‑1 (MQ=255)
          ccGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGGGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTg                                                                                                                                            <  2:1246122/149‑1 (MQ=255)
           cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGGGGGGTGTCCGGGGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTgg                                                                                                                                           <  1:300628/149‑2 (MQ=255)
                 ccACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                                                                                                                     <  1:1320079/149‑1 (MQ=255)
                 ccACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGATCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                                                                                                                     <  2:688144/149‑1 (MQ=255)
                  cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCCGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTg                                                                                                                                    >  1:662286/1‑149 (MQ=255)
                  cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTGGTCGTCATCTATGCCCCACCGAATGCCGGAGTCGTTGACCCACTg                                                                                                                                    >  1:131893/1‑149 (MQ=255)
                   aCTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTgg                                                                                                                                   >  1:1103611/1‑149 (MQ=255)
                     ttgCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGGGGGTTGGCCGGGGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCg                                                                                                                                 <  1:986077/147‑1 (MQ=255)
                         ccTGCGCGCTCTCCATCGCGTAGGCCATCAAGATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTAATCTATGCCCCTACGAATGCCGGATTCGTTGACCCACTGGCGTaaa                                                                                                                             <  1:1208367/149‑1 (MQ=255)
                            gcgcgcTCTCCATCGCGTAGGCCAACAACATACGCGTCTGGCCAAATCGCGAGGTGTGTTGTCCGGTGGCGGGGAGTTGCAGCTCGCGAACGTTGTCGGCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATaa                                                                                                                          <  2:912595/149‑1 (MQ=255)
                            gcgcgcTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATaa                                                                                                                          <  1:572980/149‑1 (MQ=255)
                             cgcgcTCTCAATCTCGTAGCCCAACAACAAACGCGTCAGGGCAAATCGCCAGGTGTGTTTTCCGTTGGGGGTGAGTTCCTGCTCGCGAACGTTGTCGTCATATATGCCCCAACGTATGCCGGATTCGTTGACCCTCTGGCGTTAATatc                                                                                                                         <  2:546849/149‑3 (MQ=255)
                                gcTCTCCATCGCGTAGCCCAACCACCAACGCGGCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGGTGACCCACTGGCGTAAATAACGCa                                                                                                                      <  2:786669/149‑1 (MQ=255)
                                    tcCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGccc                                                                                                                  <  2:1264033/149‑1 (MQ=255)
                                     ccATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGCGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATACCGCAGcccg                                                                                                                 >  1:485331/1‑148 (MQ=255)
                                      cATCGCGTAGCCCAACAACATACGCGGCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCctt                                                                                                                <  1:634456/149‑1 (MQ=255)
                                          gcgTAGCCCAACAACATACGCGTCCGGCCAAAGCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTGCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCttcttc                                                                                                            <  1:1040160/149‑1 (MQ=255)
                                          gcgTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTTCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCttcgtc                                                                                                            >  2:705968/1‑149 (MQ=255)
                                          gcgTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAAAGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGTCCCACTGGCGTAAATAACGCAGCCCgtcttc                                                                                                            >  1:1463536/1‑149 (MQ=255)
                                           cgTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCg                                                                                                           <  1:1180505/149‑1 (MQ=255)
                                             gaGCCCAACAACATACGCGTGAGGCCAAATCTCCAGGTGTGTTGGCCGGGGGCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGt                                                                                                         <  1:1314994/148‑1 (MQ=255)
                                                cccAACAACATACGCGTCAGGCCAACTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGTATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGAt                                                                                                      <  1:816786/149‑1 (MQ=255)
                                                cccAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGAt                                                                                                      <  1:823506/149‑1 (MQ=255)
                                                accATCACCATACGCGTCTGGACTAATCGCAAGGTGTTTTGTGCGGTGGCGGGGAGTTCCAGCTCGAGATCGTTGTCTTCTTCTATGCCCCAACGAATTCCGGATTGGTTGACCAACTGGCGTAAATAACGCAGCCCTGCTTCGGTTAt                                                                                                      <  2:1384403/148‑1 (MQ=255)
                                                   aacaacATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCAGCTATGCCCCAACGAATGCCGGCTTCGGTGACCCACTGGCGTGAAGAACGCAGCCCTTCTGCGGTGCTGTc                                                                                                   >  2:489229/1‑149 (MQ=255)
                                                        caTACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGCGTTCCAGCTCGCTAACGTTGTCTTCATCTATGCCCCAACGAATGACGGATTCGTTGACCCACTGGCGTAAATACCGCAGCCCTTCTTCGGAGATGTCAAAcc                                                                                              >  1:1375127/1‑149 (MQ=255)
                                                            cgcgTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGGGCCAGCTCGCGAACGTTGTCGTCAGCTAGGCCCCAAAGAATGCCGGAGTCGTTGACCCACTGGCGTAAATACCGCAGCCCTTCTTCGGTGATGTCACCCCgcgc                                                                                          <  1:1268435/149‑1 (MQ=255)
                                                                  aGGCCAAATCGCCAGGTGAGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAg                                                                                    <  1:809405/149‑1 (MQ=255)
                                                                              cAGGTATGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATc                                                                        <  1:1533419/149‑1 (MQ=255)
                                                                               aGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       <  2:438491/149‑1 (MQ=255)
                                                                               aGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       >  1:277515/1‑149 (MQ=255)
                                                                                ggTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGCCCCACTGGCTTAAATACCGCCGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATccc                                                                      >  1:878811/1‑148 (MQ=255)
                                                                                ggTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCa                                                                      >  1:1391916/1‑149 (MQ=255)
                                                                                 gtgtgtTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGCCCCACTGGCGTAAATAACGCAGCCCGTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAg                                                                     >  1:648782/1‑149 (MQ=255)
                                                                                 gtgtgtTGTCCGGTGGCGGGGAGGTCCAGCTGGCGAACGTTGTCGTCATCTATGCCGCAACGAATGCCGGAGACGTTGACCCCCTGGCGTAAATAACGCAGCCCTTCTTCGGTGCTGGCAAACCGCGCCGCCAGCACCGGCACATCCAg                                                                     <  2:862397/149‑1 (MQ=255)
                                                                                    tgttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACACACTGGCGTAAATAACGCAGCCCTTCTTCGGTGAAGTCAAACCGCGCCGCCAGCACCGGCACAACCCGCac                                                                  >  1:1185783/1‑148 (MQ=255)
                                                                                     gttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAAc                                                                 >  1:448263/1‑149 (MQ=255)
                                                                                      ttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAACCCGCGCCGCCAGCACCGGCCCATCCAGCCACg                                                                >  1:489100/1‑149 (MQ=255)
                                                                                        gtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGcc                                                              >  1:765384/1‑149 (MQ=255)
                                                                                         tCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCa                                                             >  1:934572/1‑149 (MQ=255)
                                                                                           cGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATGCAGCAACGCCAGc                                                           >  2:784198/1‑149 (MQ=255)
                                                                                           cGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGc                                                           >  2:537382/1‑149 (MQ=255)
                                                                                            ggTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGTTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGca                                                          <  2:596599/149‑1 (MQ=255)
                                                                                                gCGGGGAGGGCCCGCTCGCGAACGTTGTCGTGAGCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGACCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATc                                                      <  2:1463536/149‑1 (MQ=255)
                                                                                                     gAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTg                                                 <  1:553411/149‑1 (MQ=255)
                                                                                                      aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTTATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGa                                                <  2:179622/149‑1 (MQ=255)
                                                                                                      aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGAAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGa                                                <  2:1103611/149‑1 (MQ=255)
                                                                                                       gTTCCAGCGCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                               <  2:194884/149‑1 (MQ=255)
                                                                                                        ttCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCCCATCCTCTGaca                                              >  1:1289632/1‑149 (MQ=255)
                                                                                                        ttCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGaca                                              <  2:415981/149‑1 (MQ=255)
                                                                                                         tCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGacac                                             <  2:809981/149‑1 (MQ=255)
                                                                                                         tCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGACGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGacac                                             <  2:662286/149‑1 (MQ=255)
                                                                                                          ccAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGAcaca                                            >  1:1553936/1‑149 (MQ=255)
                                                                                                                 gcgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCGCTGACACAAAACGAc                                     >  1:899644/1‑149 (MQ=255)
                                                                                                                   gAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGAGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACAACCAGCAGCGCCAGCCCATCCACAGACACAAACCGAcag                                   >  2:750666/1‑147 (MQ=255)
                                                                                                                   gAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAGTCGTTGACCCACTGGCGTAAATAGGGCAGCCCTTCTTCGGTGATGTCACACCGCGCCGCCAGCACCGGCAGAGCCAGCAACGCCAGCACGTCCTCTGACACAAACCGACTg                                   >  2:497416/1‑149 (MQ=255)
                                                                                                                    aaCGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCACCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGt                                  <  2:273493/149‑1 (MQ=255)
                                                                                                                        ttGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCATTCTTCGGTGATGTCAAACCGCGCCGCCAGCCCCGGCACATCCAGCACCGCCAGCACTTCCTCTGACACAAAACGACGGTCAgg                              >  1:907444/1‑149 (MQ=255)
                                                                                                                           tcgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCACCGCCAGCACATCCTCTGACACAAAACGCCTCTCAGGCAg                           >  2:906757/1‑149 (MQ=255)
                                                                                                                             gtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTg                         <  2:765384/149‑1 (MQ=255)
                                                                                                                               cATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGAt                       <  2:474046/149‑1 (MQ=255)
                                                                                                                                   tATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTACTTCGGTGATGTAAAACCGCGCCGACAGCACCGGCACAACCAGCAACGCCAGCACATCCTCTGAAAAAAAACGACTGTCATGCAGTGATaaaa                   >  2:447063/1‑147 (MQ=255)
                                                                                                                                    aTGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAg                  >  2:1556343/1‑149 (MQ=255)
                                                                                                                                      gCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGCCACAAAACGACTGTCAGGCAGTGATAACAGGc                >  2:115359/1‑149 (MQ=255)
                                                                                                                                       ggtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCt               >  2:1460394/4‑149 (MQ=255)
                                                                                                                                        cccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg              >  1:324198/1‑149 (MQ=255)
                                                                                                                                         tcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCCGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTATGACACCAAACGCCTGGCAGGCAGTGATAACCGGCAGa             >  1:1030352/2‑149 (MQ=255)
                                                                                                                                         ccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCACCGCCAGCACATCCTCTGACACAAAACGCCTGTCAGGCAGTGATAACAGGCTGa             >  1:1471763/1‑149 (MQ=255)
                                                                                                                                          cAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCGTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAATCGACTGTCAGGCAGTGATAACAGGCTGAt            >  2:787081/1‑149 (MQ=255)
                                                                                                                                             cgcaTGCGGGCTTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGGCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATaaa         <  1:1559751/146‑1 (MQ=255)
                                                                                                                                             cGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGCCTGTTATCACTGCCTGacagt         <  1:705968/149‑5 (MQ=255)
                                                                                                                                                  gCCGGATTCGTTGCCCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCCCATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCtt    >  1:1471947/1‑149 (MQ=255)
                                                                                                                                                   ccGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTc   <  2:277515/149‑1 (MQ=255)
                                                                                                                                                    cGGATTCGTTGAGCCAGTGGCGTAAATAACGCATGCCTTATTAGGTTATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACTAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTcc  <  2:1375127/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960780‑2961076

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