Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,858,564 A→C 8.1% intergenic (‑25/+200) PP_mr50 ← / ← ligA Bacteria_large_SRP ncRNA/DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,858,5640AC 8.1% 195.4 / 6.2 73intergenic (‑25/+200)PP_mr50/ligABacteria_large_SRP ncRNA/DNA ligase
Reads supporting (aligned to +/- strand):  ref base A (37/30);  new base C (4/2);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

CGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGG  >  NC_002947/4858422‑4858713
                                                                                                                                              |                                                                                                                                                     
cGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTATCGCTGAGCAGCGa                                                                                                                                                >  2:1042546/1‑150 (MQ=255)
      cAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGa                                                                                                                                                       <  1:1078530/137‑1 (MQ=255)
      cAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGa                                                                                                                                                       >  2:1078530/1‑137 (MQ=255)
      cAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGgcg                                                                                                                                          >  1:1518211/1‑150 (MQ=255)
              gAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTa                                                                                                                                  >  1:1935503/1‑150 (MQ=255)
              gAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTa                                                                                                                                  >  1:1935424/1‑150 (MQ=255)
                      cccACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGCTCAGCGAgacat                                                                                                                                           >  1:1192970/1‑129 (MQ=255)
                      cccACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGCTCAGCGAgacat                                                                                                                                           <  2:1192970/133‑5 (MQ=255)
                      cccACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCtt                                                                                                                          >  1:750627/1‑150 (MQ=255)
                      cccACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCtt                                                                                                                          >  2:1029510/1‑150 (MQ=255)
                      cccACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCtt                                                                                                                          >  2:1922424/1‑150 (MQ=255)
                       ccaccaGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCttt                                                                                                                         >  1:2066892/1‑150 (MQ=255)
                       ccaccaGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCttt                                                                                                                         <  1:2164567/150‑1 (MQ=255)
                        caccaGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTg                                                                                                                        <  2:2194736/150‑1 (MQ=255)
                           caGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAg                                                                                                                     <  2:1541120/150‑1 (MQ=255)
                             tccACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCTGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCa                                                                                                                   <  1:262829/149‑1 (MQ=255)
                                acacCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCtt                                                                                                                          >  2:2282441/1‑140 (MQ=255)
                                acacCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCtt                                                                                                                          <  1:2282441/140‑1 (MQ=255)
                                  acCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGAACGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATAc                                                                                                              >  2:804958/1‑150 (MQ=255)
                                    cccGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACaa                                                                                                            <  1:1029510/150‑1 (MQ=255)
                                    cccGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACaa                                                                                                            <  2:2257034/150‑1 (MQ=255)
                                     ctggCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAAc                                                                                                           >  1:448862/3‑150 (MQ=255)
                                     ccGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAAc                                                                                                           <  2:1935503/150‑1 (MQ=255)
                                     ccGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAAc                                                                                                           <  2:1935424/150‑1 (MQ=255)
                                       ggCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCa                                                                                                         >  1:2954/1‑150 (MQ=255)
                                           cTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGc                                                                                                     >  2:1819262/1‑150 (MQ=255)
                                              gATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTg                                                                                                                                             <  2:907042/107‑1 (MQ=255)
                                              gATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTg                                                                                                                                             >  1:907042/1‑107 (MQ=255)
                                              gATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCtt                                                                                                                          <  2:1875439/126‑1 (MQ=255)
                                              gATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCtt                                                                                                                          >  1:1875439/1‑126 (MQ=255)
                                               cttttCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCa                                                                                                 <  2:976617/147‑1 (MQ=255)
                                                 gTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGa                                                                                               >  2:628071/1‑150 (MQ=255)
                                                  ttCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGaa                                                                                              >  2:92269/1‑150 (MQ=255)
                                                     ccGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGaaaaa                                                                                           >  2:14611/1‑150 (MQ=255)
                                                     ccGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGaaaaa                                                                                           >  1:1441594/1‑150 (MQ=255)
                                                     ccGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTACCAATCAAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGaaaaa                                                                                           >  1:25020/1‑150 (MQ=255)
                                                      cGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGCTCAGCGAgacat                                                                                                                                           <  2:233672/101‑5 (MQ=255)
                                                      cGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGCTCAGCGAgacat                                                                                                                                           >  1:233672/1‑97 (MQ=255)
                                                             gctgctCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTcatcat                                                                                   >  2:1173189/1‑150 (MQ=255)
                                                                 ctCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGgc                                                                                                                                           <  2:554773/90‑1 (MQ=255)
                                                                 ctCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGgc                                                                                                                                           >  1:554773/1‑90 (MQ=255)
                                                                    ccTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCaa                                                                            >  2:1332375/1‑150 (MQ=255)
                                                                     cTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAAt                                                                           <  2:1518211/150‑1 (MQ=255)
                                                                      ttCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATg                                                                          >  2:1280738/1‑150 (MQ=255)
                                                                         cGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACt                                                                       >  1:1502252/1‑150 (MQ=255)
                                                                            gCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGt                                                                    <  2:2066892/150‑1 (MQ=255)
                                                                              cTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCa                                                                  >  1:1043033/1‑150 (MQ=255)
                                                                                   cAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTc                                                             <  2:1502252/150‑1 (MQ=255)
                                                                                      gTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCgag                                                          >  2:2120885/1‑150 (MQ=255)
                                                                                      gTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCgag                                                          <  2:448862/150‑1 (MQ=255)
                                                                                         aaCGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATg                                                       <  2:1271574/150‑1 (MQ=255)
                                                                                            gggTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGcc                                                    >  1:465248/1‑150 (MQ=255)
                                                                                             ggTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCt                                                   >  1:1728234/1‑150 (MQ=255)
                                                                                              gTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTg                                                  <  2:25020/150‑1 (MQ=255)
                                                                                               tAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTgg                                                 <  1:628071/150‑1 (MQ=255)
                                                                                                  tCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACg                                              >  2:2079203/1‑150 (MQ=255)
                                                                                                   cAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACgg                                             <  1:1173189/150‑1 (MQ=255)
                                                                                                   cAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACgg                                             <  1:2079203/150‑1 (MQ=255)
                                                                                                    aaTGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACggg                                            >  1:2001093/1‑150 (MQ=255)
                                                                                                    aaTGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACggg                                            >  1:307762/1‑150 (MQ=255)
                                                                                                             ggACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAGCCACTGCGGCAGAAAAAGTCGTCATGCGTCAATGACGTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTcgc                                   >  1:541712/1‑150 (MQ=255)
                                                                                                               aCCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCgg                                 <  1:14611/150‑1 (MQ=255)
                                                                                                               aCCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCgg                                 <  1:694599/150‑1 (MQ=255)
                                                                                                                  gATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTa                              >  2:83422/1‑150 (MQ=255)
                                                                                                                     ggggCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAAc                           <  1:92269/150‑1 (MQ=255)
                                                                                                                     ggggCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCACCAttt                                                                                 <  2:2062789/96‑1 (MQ=255)
                                                                                                                     ggggCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCACCAttt                                                                                 >  1:2062789/1‑96 (MQ=255)
                                                                                                                            ccacTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAg                                                                                          <  2:1881283/80‑1 (MQ=255)
                                                                                                                            ccacTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAg                                                                                          >  1:1881283/1‑80 (MQ=255)
                                                                                                                                 aCATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACgg               <  2:1842012/150‑1 (MQ=255)
                                                                                                                                  cATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACggg              <  1:1280738/150‑1 (MQ=255)
                                                                                                                                   aTGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGc             >  1:141893/1‑150 (MQ=255)
                                                                                                                                    ccactcGCTGCTCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCg            >  2:80325/4‑150 (MQ=255)
                                                                                                                                      tctcGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGAc          <  1:1042546/150‑1 (MQ=255)
                                                                                                                                        tcGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACcg        >  2:1185031/1‑150 (MQ=255)
                                                                                                                                              agcagcGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAgg  <  2:465248/150‑1 (MQ=255)
                                                                                                                                              agcagcGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTACTACGGGCGACCGCGCAgg  <  1:83422/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
CGACCCCAGAAATGGAGGCGGCCCCACCAGCCACACCCCGGCACTCGATGTTCCCGCTATGGCTGCTCCCTTCCGGGCCTGACCAGGTTAACGGGTAATCAATGCGGGGGGACCGATGGGGCCACCACTACATGTCTCGCTGAGCAGCGAGTGGCGCCATTGTACCGGTCTTTGGAGCAGATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGG  >  NC_002947/4858422‑4858713

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: