Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,273,113 C→G 9.3% intergenic (+35/+256) PP_2003 → / ← PP_2004 hypothetical protein/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,273,1130CG 9.3% 234.4 / 13.2 75intergenic (+35/+256)PP_2003/PP_2004hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (33/35);  new base G (4/3);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGC  >  NC_002947/2272968‑2273259
                                                                                                                                                 |                                                                                                                                                  
gtggtgGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGa                                                                                                                                                >  1:588246/1‑150 (MQ=255)
   gtggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAg                                                                                                                                             <  1:1618956/150‑1 (MQ=255)
    tggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAgg                                                                                                                                            >  1:1220515/1‑150 (MQ=255)
    tggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAgg                                                                                                                                            >  1:637129/1‑150 (MQ=255)
      gATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGc                                                                                                                                          <  1:1420990/150‑1 (MQ=255)
        tGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCg                                                                                                                                        <  1:1553079/150‑1 (MQ=255)
          aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCGGt                                                                                                                                      <  1:743222/150‑1 (MQ=38)
          aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGt                                                                                                                                      >  2:112852/1‑150 (MQ=255)
             gCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCa                                                                                                                                   >  1:973117/1‑150 (MQ=255)
              ccAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAg                                                                                                                                  >  2:290407/1‑150 (MQ=255)
               cAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAgg                                                                                                                                 <  2:263254/150‑1 (MQ=255)
               cAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAgg                                                                                                                                 <  1:157725/150‑1 (MQ=255)
                       gCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAAc                                                                                                                         >  1:555994/1‑150 (MQ=255)
                       gCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAAc                                                                                                                         <  2:1220515/150‑1 (MQ=255)
                              gtgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCt                                                                                                                  >  1:1261905/1‑150 (MQ=255)
                                ggCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGt                                                                                                                <  1:279902/150‑1 (MQ=255)
                                      gccggtcgAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGa                                                                                                          <  2:205881/150‑1 (MQ=255)
                                              aaGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCcgc                                                                                                  <  1:290344/150‑1 (MQ=255)
                                                       ccATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACttt                                                                                         >  2:770110/1‑150 (MQ=255)
                                                         aTTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCt                                                                                       <  2:637129/150‑1 (MQ=255)
                                                              ggTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAaca                                                                                  >  1:191034/1‑150 (MQ=255)
                                                                  cagcagGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTg                                                                              >  1:943693/1‑150 (MQ=255)
                                                                    gcagGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTCGATGGGTCGCCGGAGCCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCt                                                                            >  1:31542/1‑150 (MQ=255)
                                                                    gcagGATGCCCAGCCGACCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCt                                                                            >  1:1512132/1‑150 (MQ=255)
                                                                           gCCCAGCCGACCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTa                                                                     <  2:1512132/150‑1 (MQ=255)
                                                                            cccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAc                                                                    <  1:290407/150‑1 (MQ=255)
                                                                             ccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGCTAACTTTCTGAACACCTGCTACGTTTAcc                                                                   >  1:375091/1‑150 (MQ=255)
                                                                                gccggccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACt                                                                >  1:1497022/1‑150 (MQ=255)
                                                                                 ccggccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTc                                                               <  1:1110953/150‑1 (MQ=255)
                                                                                  cggccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCa                                                              >  1:216536/1‑150 (MQ=255)
                                                                                     ccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATaa                                                           <  1:112852/150‑1 (MQ=255)
                                                                                       gAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGc                                                         >  2:736463/1‑150 (MQ=255)
                                                                                          gccgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTc                                                                                        >  1:948200/1‑116 (MQ=255)
                                                                                          gccgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTCTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTc                                                                                        <  2:948200/116‑1 (MQ=255)
                                                                                           ccgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGt                                                     <  2:943693/150‑1 (MQ=255)
                                                                                           ccgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGt                                                     <  1:800636/150‑1 (MQ=255)
                                                                                              ccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGc                                                         <  2:1315437/143‑1 (MQ=255)
                                                                                              ccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGc                                                         >  1:1315437/1‑143 (MQ=255)
                                                                                                  ccccTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACCACTGCTACGTTTACCACTCATAAGCTTGTCatcatc                                              >  2:428611/1‑150 (MQ=255)
                                                                                                      tccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCtt                                          <  2:588246/150‑1 (MQ=255)
                                                                                                              aTGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTaaa                                                                                                                          <  1:1261362/62‑1 (MQ=255)
                                                                                                              aTGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTaaa                                                                                                                          >  2:1261362/1‑62 (MQ=255)
                                                                                                              aTGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCAc                                  >  2:341466/1‑150 (MQ=255)
                                                                                                                  gTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCtt                                          <  1:436522/138‑1 (MQ=255)
                                                                                                                  gTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCtt                                          >  2:436522/1‑138 (MQ=255)
                                                                                                                  gTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGc                              >  1:1552572/1‑150 (MQ=255)
                                                                                                                     gccgggaCCGTCCCTTTCGCCGGCAGCCCGTGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGGGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTa                           <  2:1059641/144‑1 (MQ=255)
                                                                                                                     gCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTa                           >  1:3378/1‑150 (MQ=255)
                                                                                                                     gCCTGAACCGACCCTTTCGCCGCCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCt                                                                                       >  2:781068/1‑90 (MQ=255)
                                                                                                                     gCCTGAACCGACCCTTTCGCCGCCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCt                                                                                       <  1:781068/90‑1 (MQ=255)
                                                                                                                       cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACttt                                                                                         <  1:302064/86‑1 (MQ=255)
                                                                                                                       cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACttt                                                                                         >  2:302064/1‑86 (MQ=255)
                                                                                                                       cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAgg                         <  2:1518459/150‑1 (MQ=255)
                                                                                                                        tGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGc                        <  1:770110/150‑1 (MQ=255)
                                                                                                                         gAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCt                       >  1:1509950/1‑150 (MQ=255)
                                                                                                                             cggcccTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTa                           >  2:1618245/4‑142 (MQ=38)
                                                                                                                             cggcccTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTa                           <  1:1618245/139‑1 (MQ=38)
                                                                                                                             cggcccTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTAAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAg                   >  1:380476/4‑150 (MQ=38)
                                                                                                                             cGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAg                   <  1:1332508/150‑1 (MQ=255)
                                                                                                                              ggcccTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAgg                  >  2:454795/3‑150 (MQ=39)
                                                                                                                                 ccTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCt                                                                            >  2:1451914/1‑89 (MQ=255)
                                                                                                                                 ccTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCt                                                                            <  1:1451914/89‑1 (MQ=255)
                                                                                                                                 ccTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGa               >  1:981197/1‑150 (MQ=255)
                                                                                                                                 ccTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGa               <  1:213368/150‑1 (MQ=255)
                                                                                                                                   tttCGCGGGCTGCCGGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTc                                                                                        <  2:47489/75‑1 (MQ=25)
                                                                                                                                   tttCGCGGGCTGCCGGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTc                                                                                        >  1:47489/1‑75 (MQ=37)
                                                                                                                                   tttCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTg                                                      >  2:454530/1‑109 (MQ=255)
                                                                                                                                   tttCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTg                                                      <  1:454530/109‑1 (MQ=255)
                                                                                                                                   tttCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATa             >  1:591223/1‑150 (MQ=255)
                                                                                                                                       gCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTa         <  2:31542/150‑1 (MQ=255)
                                                                                                                                       gCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTa         <  2:216536/150‑1 (MQ=255)
                                                                                                                                          ggCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCcgc                                                                                                  >  1:924083/1‑58 (MQ=255)
                                                                                                                                          ggCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCcgc                                                                                                  <  2:924083/58‑1 (MQ=255)
                                                                                                                                            cAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAg    <  2:1261905/150‑1 (MQ=255)
                                                                                                                                              gCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGc  <  1:341466/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAGGC  >  NC_002947/2272968‑2273259

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