Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,635,170 A→G 7.6% A313A (GCA→GCG xseA → exodeoxyribonuclease VII subunit XseA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,635,1700AG 7.6% 194.0 / 5.4 69A313A (GCA→GCGxseAexodeoxyribonuclease VII subunit XseA
Reads supporting (aligned to +/- strand):  ref base A (36/27);  new base G (3/2);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GCCGTAATCAGCAAGAGTTACTGCGCCAGGTGCAATCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGAAGATTCATCGCGCGCAAACGCG  >  NC_000913/2635022‑2635319
                                                                                                                                                    |                                                                                                                                                     
gCCGTAATCAGCAAGAGTTACTGCGCCAGGTGCAATCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCAc                                                                                                                                                      >  1:1947436/1‑150 (MQ=255)
 ccGTAATCAGCAAGAGTTACTGCGCCAGGTGCAATCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACg                                                                                                                                                     >  1:829784/1‑150 (MQ=255)
        cAGCAAGAGTTACTGCGCCAGGTGCAATCGACCCGTCCACGGCTGGAGATGGCGCTGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCAGCAGCTCCGGCTGGCCCGCCAGCac                                                                                                                                              >  2:1991255/1‑149 (MQ=255)
         aGCAAGAGTTACTGCGCCAGGTGCAATCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCGTGCCAGccgg                                                                                                                                             <  1:614196/150‑4 (MQ=255)
            aaGAGTTACTGCGCCAGGTGCAATCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCGTGCCAGccggagc                                                                                                                                          <  1:1575167/150‑7 (MQ=255)
            aaGAGTTACTGCGCCAGGTGCAATCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCa                                                                                                                                          <  2:189008/150‑1 (MQ=255)
                    cTGCGCCAGGTGCAATCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAg                                                                                                                                  <  1:721930/150‑1 (MQ=255)
                         ccAGGTGCAATCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCt                                                                                                                             >  2:525847/1‑150 (MQ=255)
                                  aTCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCg                                                                                                                    <  1:826304/150‑1 (MQ=255)
                                    cGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGaa                                                                                                                  >  2:487398/1‑150 (MQ=255)
                                      aCCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATg                                                                                                                <  2:120543/150‑1 (MQ=255)
                                      aCCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATg                                                                                                                <  2:1365279/150‑1 (MQ=255)
                                        ccGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAg                                                                                                              <  1:280512/150‑1 (MQ=255)
                                               cGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTgcg                                                                                                       >  1:461700/1‑150 (MQ=255)
                                                ggCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTgcgc                                                                                                      >  1:260980/1‑150 (MQ=255)
                                                 gCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCt                                                                                                     >  1:1745078/1‑150 (MQ=255)
                                                  cTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTg                                                                                                    >  1:1111544/1‑150 (MQ=255)
                                                  cTGGAGATGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTgg                                                                                                   >  2:1327644/1‑150 (MQ=255)
                                                     gagaTGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGaa                                                                                                 <  2:260980/150‑1 (MQ=255)
                                                     gagaTGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGaa                                                                                                 <  1:1387245/150‑1 (MQ=255)
                                                      agaTGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGaaa                                                                                                >  2:237091/1‑150 (MQ=255)
                                                           gCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCaa                                                                                           >  1:686113/1‑150 (MQ=255)
                                                            cGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAAc                                                                                          >  1:1924254/1‑150 (MQ=255)
                                                             gATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACt                                                                                         <  2:1321819/150‑1 (MQ=255)
                                                             gATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACt                                                                                         <  1:525847/150‑1 (MQ=255)
                                                                 gATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAg                                                                                     <  2:1745078/150‑1 (MQ=255)
                                                                       tatCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTgc                                                                                                        <  2:457253/125‑1 (MQ=255)
                                                                       tatCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTgc                                                                                                        >  1:457253/1‑125 (MQ=255)
                                                                       tatCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTAcc                                                                               <  2:1886419/150‑1 (MQ=255)
                                                                       tatCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTAcc                                                                               >  1:587491/1‑150 (MQ=255)
                                                                                aaCCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGCCAAccc                                                                      >  1:1436872/1‑148 (MQ=255)
                                                                                 aCCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAAcagc                                                                     >  1:369604/1‑150 (MQ=255)
                                                                                 aCCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAAcagc                                                                     <  1:487398/150‑1 (MQ=255)
                                                                                  ccGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAAcagcc                                                                    >  1:1593750/1‑149 (MQ=255)
                                                                                       cacGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACagcagcagc                                                               >  1:1090206/1‑150 (MQ=255)
                                                                                         cGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCgg                                                             >  1:967808/1‑150 (MQ=255)
                                                                                           tCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGtt                                                           >  1:98500/1‑150 (MQ=255)
                                                                                            cGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTa                                                          <  2:1947436/150‑1 (MQ=255)
                                                                                            cGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTa                                                          <  2:461700/150‑1 (MQ=255)
                                                                                              cTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAac                                                        <  2:888612/150‑1 (MQ=255)
                                                                                                ttACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAacac                                                      >  2:333367/1‑150 (MQ=255)
                                                                                                 tACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAAcaca                                                     >  2:1254941/1‑150 (MQ=255)
                                                                                                   cGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGc                                                   <  2:760497/150‑1 (MQ=255)
                                                                                                     cAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCa                                                                          <  1:1557917/125‑1 (MQ=255)
                                                                                                     cAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGAGTACCGGGCa                                                                          >  2:1557917/1‑125 (MQ=255)
                                                                                                     cAAAGTCATCACCGAAAACAGCAACAGCATCCGCAGCTACCGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCATTGCGCTGGAAAATCAACTTAAGCCTACCGGCCAACAGCAGCAGCGGTGACCACAGCgg                                                 >  2:859875/1‑150 (MQ=255)
                                                                                                          tcatcaCCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGaa                                            <  1:1699695/150‑1 (MQ=255)
                                                                                                             tcaCCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCa                                         >  2:1746217/1‑150 (MQ=255)
                                                                                                             tcaCCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCa                                         >  1:147741/1‑150 (MQ=255)
                                                                                                              caCCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAg                                        >  1:528864/1‑150 (MQ=255)
                                                                                                                    ttACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAAt                                  <  2:1431617/150‑1 (MQ=255)
                                                                                                                     tACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATc                                 >  2:1893938/1‑150 (MQ=255)
                                                                                                                      acagcaacagcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATcc                                <  2:587491/150‑1 (MQ=255)
                                                                                                                      acagcaacagcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATcc                                >  2:664842/1‑150 (MQ=255)
                                                                                                                      acagcaacagcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATcc                                >  2:671611/1‑150 (MQ=255)
                                                                                                                       cagcaacagcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCca                               >  1:1226449/1‑150 (MQ=255)
                                                                                                                           aacagcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGc                           <  2:1372977/150‑1 (MQ=255)
                                                                                                                           aacagcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGc                           >  2:1057036/1‑150 (MQ=255)
                                                                                                                            acagcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGcc                          <  2:31471/150‑1 (MQ=255)
                                                                                                                               gcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACTGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGaa                       <  1:1811865/150‑1 (MQ=255)
                                                                                                                               gcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTa                                                          >  2:1022258/1‑115 (MQ=255)
                                                                                                                               gcaTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTa                                                          <  1:1022258/115‑1 (MQ=255)
                                                                                                                                 aTCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGAAGa                     >  1:1878075/1‑150 (MQ=255)
                                                                                                                                   ccGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGAAGAtt                   <  1:510777/150‑1 (MQ=255)
                                                                                                                                       atggtttgCTGGCGTGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGAAGATTCATc               >  1:679179/8‑150 (MQ=255)
                                                                                                                                         cTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGAAGATTCATcgc             >  1:174433/1‑150 (MQ=255)
                                                                                                                                          gtttgCTGGCGTGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGAAGATTCATcgcg            >  2:1485706/5‑150 (MQ=255)
                                                                                                                                           tttgCTGGCGTGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCCGCAGCGGTTAACACAGCGGCTGACTCAGCAACATCCACAGCCGAAGATTCATcgcgc           >  2:1685227/4‑150 (MQ=255)
                                                                                                                                                   cACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGAAGATTCATCGCGCGCAAAcgc   >  1:975966/1‑150 (MQ=255)
                                                                                                                                                    aCGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGAAGATTCATCGCGCGCAAAcgcg  >  1:1630883/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGTAATCAGCAAGAGTTACTGCGCCAGGTGCAATCGACCCGTCAACGGCTGGAGATGGCGATGGATTATTATCTCGCCAACCGCACACGTCGCTTTACGCAAATTCATCACCGATTACAGCAACAGCATCCGCAGCTCCGGCTGGCACGCCAGCAAACCATGCTTGAGCGCCTGCAAAAGCGAATGAGCTTTGCGCTGGAAAATCAACTTAAGCGTACCGGGCAACAGCAGCAGCGGTTAACACAGCGGCTGAATCAGCAAAATCCACAGCCGAAGATTCATCGCGCGCAAACGCG  >  NC_000913/2635022‑2635319

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: