Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,960,930 G→C 5.2% S500C (TCC→TGC)  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,9300GC 5.2% 212.2 / 4.5 77S500C (TCC→TGC) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base C (2/2);  total (36/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 = 4.93e-01

AGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960782‑2961076
                                                                                                                                                    |                                                                                                                                                  
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                                                                       >  1:277515/1‑149 (MQ=255)
                                                                             aGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       <  2:438491/149‑1 (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)
                                                                                                                                                    |                                                                                                                                                  
AGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960782‑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: ‑