Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 495,040 C→T 8.8% intergenic (+30/‑80) recR → / → htpG DNA repair protein RecR/chaperone protein HtpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913495,0400CT 8.8% 371.8 / 16.5 124intergenic (+30/‑80)recR/htpGDNA repair protein RecR/chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base C (55/58);  new base T (4/7);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTC  >  NC_000913/494894‑495189
                                                                                                                                                  |                                                                                                                                                     
tatGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTg                                                                                                                                                    <  1:775120/150‑1 (MQ=255)
 atGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCTCTCGTTTgc                                                                                                                                                   <  1:53459/150‑2 (MQ=255)
  tGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaa                                                                                                                                                  >  2:410945/1‑150 (MQ=255)
  tGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTGGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGTGATCCTGCTCTCGTTTgct                                                                                                                                                  <  1:299830/150‑3 (MQ=255)
    aCGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGTGATCCTGCTCTCGTTTgctta                                                                                                                                                <  1:1072193/150‑5 (MQ=255)
     cGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAatt                                                                                                                                               >  1:1626697/1‑150 (MQ=255)
     cGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAatt                                                                                                                                               <  1:1473273/150‑1 (MQ=255)
     cGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAatt                                                                                                                                               <  2:558086/150‑1 (MQ=255)
       gggAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAattat                                                                                                                                             <  1:951770/149‑1 (MQ=255)
        ggAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGCCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGGAattatt                                                                                                                                            >  1:1079869/1‑150 (MQ=255)
             ccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTccc                                                                                                                                       <  2:491096/150‑1 (MQ=255)
              cagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCTCTCGTTTgcttaaa                                                                                                                                              <  1:1361576/142‑7 (MQ=38)
              cagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCTCTCGTTTgcttaaa                                                                                                                                              >  2:1361576/1‑136 (MQ=38)
              cagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCTCTCGTTTgcttaaa                                                                                                                                              >  2:1367343/1‑136 (MQ=38)
              cagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCt                                                                                                                                      <  1:1458403/150‑1 (MQ=255)
              cagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCt                                                                                                                                      <  1:1486808/150‑1 (MQ=255)
              cagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGACTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCt                                                                                                                                      <  1:1097024/150‑1 (MQ=255)
                gccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTg                                                                                                                                    >  1:909125/1‑150 (MQ=255)
                 ccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGt                                                                                                                                   <  1:1350517/150‑1 (MQ=255)
                   gAATCGCTCATGGCGTTCCGGTTGGTGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTcc                                                                                                                                 <  2:631230/150‑1 (MQ=255)
                       cGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAt                                                                                                                             >  2:1089629/1‑150 (MQ=255)
                        gCTCATGGCGTTCCGGTTGGCGGCGAGCTGGCAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtc                                                                                                                            <  1:1582684/150‑1 (MQ=255)
                        gCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtc                                                                                                                            >  1:1282174/1‑150 (MQ=255)
                          tCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtctc                                                                                                                          >  1:1387136/1‑150 (MQ=255)
                           cATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtctct                                                                                                                         <  2:1143336/150‑1 (MQ=255)
                               gCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCca                                                                                                                     <  2:1425705/150‑1 (MQ=255)
                                    ccGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATcc                                                                                                                <  2:720941/150‑1 (MQ=255)
                                      ggTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTg                                                                                                              >  1:1436037/1‑150 (MQ=255)
                                      ggTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTg                                                                                                              >  1:390685/1‑150 (MQ=255)
                                        ttGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGtt                                                                                                            <  1:492191/150‑1 (MQ=255)
                                        ttGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGtt                                                                                                            >  2:1095492/1‑150 (MQ=255)
                                           gcggcgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGttttt                                                                                                         <  2:764844/150‑1 (MQ=255)
                                            cggcgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTa                                                                                                        >  2:1541407/1‑150 (MQ=255)
                                              gcgAGCTGGATATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAAc                                                                                                      <  1:265259/150‑1 (MQ=255)
                                               cgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAAcc                                                                                                     <  2:1216732/150‑1 (MQ=255)
                                                  gCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTa                                                                                                  >  1:1038412/1‑150 (MQ=255)
                                                      gAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAAt                                                                                              >  1:1645100/1‑150 (MQ=255)
                                                          tGGTCGACGGCACCACGTTGTCACACTCCCTTGCTGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCa                                                                                          >  2:899638/1‑150 (MQ=255)
                                                          tGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCa                                                                                          <  1:1006190/150‑1 (MQ=255)
                                                           ggTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCat                                                                                         <  1:1308916/150‑1 (MQ=255)
                                                           ggTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCat                                                                                         <  1:1051950/150‑1 (MQ=255)
                                                            gTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTTACCTTAAAATGGCatt                                                                                        >  1:1624845/1‑150 (MQ=255)
                                                               gACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGAATCGTTTTTAAGCAAACGAGAGCAGGATCACCTCCTCTCGCATGAAAATATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCattatt                                                                                     >  2:772324/1‑150 (MQ=255)
                                                                   gCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAgg                                                                                 <  1:410945/150‑1 (MQ=255)
                                                                   gCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAgg                                                                                 >  2:1355030/1‑150 (MQ=255)
                                                                   gCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAgg                                                                                 >  2:1572099/1‑150 (MQ=255)
                                                                    caccacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTTTTGAGGt                                                                                >  2:555968/1‑150 (MQ=255)
                                                                    caccacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGt                                                                                >  1:272089/1‑150 (MQ=255)
                                                                      ccacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAg                                                                              >  1:409548/1‑150 (MQ=255)
                                                                        acGTTGTCACACTACCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGAc                                                                            <  1:546641/150‑1 (MQ=255)
                                                                           ttGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTa                                                                         <  2:390685/150‑1 (MQ=255)
                                                                             gTCACACTCCCTTGCTGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACa                                                                       <  1:899638/150‑1 (MQ=255)
                                                                               cacaCTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATg                                                                     >  1:1120395/1‑150 (MQ=255)
                                                                               cacaCTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATg                                                                     <  2:409548/150‑1 (MQ=255)
                                                                                acacTCCCTTGCCGGGCGTCATAATATTATTTTTTAAGCAAACGAGAGCAGGATCACCTGCTATCGCTTGAAATTATTATCCCTTGTACACATCTCTCCAACATACTTTTTTTAACCTTAAAATTGAATAATTGAGGTAAAACTACAAGa                                                                    >  2:746135/1‑150 (MQ=255)
                                                                                   cTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAg                                                                 >  2:792172/1‑150 (MQ=255)
                                                                                   cTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAg                                                                 <  1:1541407/150‑1 (MQ=255)
                                                                                   cTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAg                                                                 >  2:1493265/1‑150 (MQ=255)
                                                                                    tCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTTTTGAGGTAGACCTACATGAAAgg                                                                <  1:555968/150‑1 (MQ=255)
                                                                                     cccTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGa                                                               <  2:826308/150‑1 (MQ=255)
                                                                                     cccTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGa                                                               >  2:763672/1‑150 (MQ=255)
                                                                                     cccTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGa                                                               >  2:1470420/1‑150 (MQ=255)
                                                                                     cccTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCCCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGa                                                               >  1:464811/1‑150 (MQ=255)
                                                                                      ccTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGAc                                                              <  2:1337339/150‑1 (MQ=255)
                                                                                       cTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACa                                                             <  2:833378/150‑1 (MQ=255)
                                                                                          gCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGa                                                          >  1:1163565/1‑150 (MQ=255)
                                                                                             gggCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCTATATTTTTTTTTTAACCTTAAAATGTCATTATTTATGTAGAACTACATTAAAGTACAAGAAAc                                                       <  2:398571/150‑1 (MQ=255)
                                                                                                cGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCg                                                    <  1:1187005/150‑1 (MQ=255)
                                                                                                  gcATAAGTTTCGTTTTTAAGCAAACGTGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATGCTATTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTg                                                  <  2:1410729/149‑1 (MQ=255)
                                                                                                   cATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTgg                                                 <  2:1282174/150‑1 (MQ=255)
                                                                                                   cATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTgg                                                 >  2:1210153/1‑150 (MQ=255)
                                                                                                   cATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTgg                                                 <  1:692105/150‑1 (MQ=255)
                                                                                                      aaGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCTCTTGAAATTATTCTCCCTTGTGCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGttt                                              <  1:1255980/150‑1 (MQ=255)
                                                                                                      aaGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGttt                                              >  2:1398718/1‑150 (MQ=255)
                                                                                                      aaGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGttt                                              <  1:1493265/150‑1 (MQ=255)
                                                                                                        gATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtc                                            <  2:1626697/150‑1 (MQ=255)
                                                                                                          ttCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                                          <  1:1095492/150‑1 (MQ=255)
                                                                                                          ttCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                                          <  2:884831/150‑1 (MQ=255)
                                                                                                          ttCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                                          <  1:429080/150‑1 (MQ=255)
                                                                                                            cGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAAAGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtc                                        >  2:627288/1‑150 (MQ=255)
                                                                                                             gTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtca                                       >  1:901606/1‑150 (MQ=255)
                                                                                                             gTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtca                                       >  1:1004120/1‑150 (MQ=255)
                                                                                                              tttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtcag                                      >  2:857445/1‑150 (MQ=255)
                                                                                                              tttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtcag                                      <  2:446800/150‑1 (MQ=255)
                                                                                                               ttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGCTTTTCAGTCAGa                                     >  2:1525952/1‑150 (MQ=255)
                                                                                                                tttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCAATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGaa                                    >  2:206553/1‑150 (MQ=255)
                                                                                                                 ttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAg                                   <  1:1355030/150‑1 (MQ=255)
                                                                                                                   tttcAAGCGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTg                                 >  2:1141218/4‑150 (MQ=255)
                                                                                                                   aaGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTg                                 >  1:1146253/1‑150 (MQ=255)
                                                                                                                      cAAGCGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCa                                                                                          <  1:1227480/90‑1 (MQ=37)
                                                                                                                      cAAGCGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCa                                                                                          >  2:1227480/1‑90 (MQ=37)
                                                                                                                       aaaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAAc                             <  2:1645100/150‑1 (MQ=255)
                                                                                                                        aaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTCCAGTCAGAAGTGAAACa                            <  2:1575439/150‑1 (MQ=255)
                                                                                                                               gCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTg                     >  1:122680/1‑150 (MQ=255)
                                                                                                                               gCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTg                     >  2:140324/1‑150 (MQ=255)
                                                                                                                                   gATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCAcc                 <  2:1004120/150‑1 (MQ=255)
                                                                                                                                            ctctCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATcc        <  2:1668455/150‑1 (MQ=255)
                                                                                                                                            ctctCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATcc        >  1:757358/1‑150 (MQ=255)
                                                                                                                                              ctcGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCAt      >  2:935197/1‑150 (MQ=255)
                                                                                                                                                cGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATtc    >  1:1547512/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTC  >  NC_000913/494894‑495189

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