Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,925 A→T 8.3% intergenic (‑112/+87) ppiC ← / ← yifN peptidyl‑prolyl cis‑trans isomerase C (rotamase C)/PemK toxin family pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,9250AT 8.3% 379.6 / 17.7 140intergenic (‑112/+87)ppiC/yifNpeptidyl‑prolyl cis‑trans isomerase C (rotamase C)/PemK toxin family pseudogene
Reads supporting (aligned to +/- strand):  ref base A (63/65);  new base T (8/4);  total (71/69)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 144 total aligned reads displayed.

TTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCA  >  NC_000913/3959776‑3960073
                                                                                                                                                     |                                                                                                                                                    
tttACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGaaa                                                                                                                                                      >  1:1587257/1‑150 (MQ=255)
  taaaaGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAGTAATTGTCAGGCGGCGATTTTAGAATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAAt                                                                                                                                                    >  1:527981/4‑150 (MQ=255)
    cAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATcc                                                                                                                                                  >  2:865231/1‑150 (MQ=255)
    cAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATcc                                                                                                                                                  <  2:714804/150‑1 (MQ=255)
     aaGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCg                                                                                                                                                 <  2:899975/150‑1 (MQ=255)
      aGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCgg                                                                                                                                                >  2:1427579/1‑150 (MQ=255)
        gATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGc                                                                                                                                              <  2:370045/150‑1 (MQ=255)
            tGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCg                                                                                                                                          >  2:1471700/1‑150 (MQ=255)
             gCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGc                                                                                                                                         >  1:1101029/1‑150 (MQ=255)
             gCAGTGCGTCTGCTGTTTTTTCCATGATCGTGCCTTGAGTGAATATTTGTCAGGCTGATATTTTAGCATTATCCTCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGATTTTCCGGGCttttg                                                                                                                                         <  2:6537/150‑4 (MQ=255)
              cAGTGCTGCTGCTGTTTTTCCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGAGTTTAGCATGATCCGCAGATGCTAACGTATTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCa                                                                                                                                        <  1:439132/150‑1 (MQ=255)
              cAGGGCTGCTGCTCTTTTTTCTATGATCGTTCCTTGAGTTAATAATTGGCAGGCGTCGATTTTAGAATGATCCGCAGATGCTAACGTAATTGGCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCa                                                                                                                                        <  2:158499/150‑1 (MQ=255)
                gTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTTTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAgg                                                                                                                                      <  1:313080/150‑1 (MQ=255)
                 tgctgctgcTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGt                                                                                                                                     >  1:1294789/1‑150 (MQ=255)
                    tgctgctgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATTATCCGCAGATGCTTACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTacc                                                                                                                                  >  2:908304/1‑150 (MQ=255)
                    tgctgctgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTccc                                                                                                                                  >  2:378391/1‑150 (MQ=255)
                         ctgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAg                                                                                                                             >  1:586968/1‑150 (MQ=255)
                          tgTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAgg                                                                                                                            <  2:1499701/150‑1 (MQ=255)
                           gTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAggg                                                                                                                           <  2:1212363/150‑1 (MQ=255)
                             ttttGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGga                                                                                                                         <  2:385337/150‑1 (MQ=255)
                                tGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaa                                                                                                                      >  1:460241/1‑150 (MQ=255)
                                 gCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                                                                                                                     >  1:1107859/1‑150 (MQ=255)
                                  ccATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaaa                                                                                                                    >  1:1162310/1‑150 (MQ=255)
                                  ccATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaaa                                                                                                                    >  1:417604/1‑150 (MQ=255)
                                   cATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCATATGCTAACGTAATTGTCAGTCAGGGAGTGTCTGGTACATATATGGCAAAAGCCCGGAAAAACCGGGCTTCGCAGGTCCCTCAAGTGGAGaaaaa                                                                                                                   >  1:650456/1‑149 (MQ=255)
                                       aTCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAACCTTAg                                                                                                               >  1:1520412/1‑150 (MQ=255)
                                         cGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGGCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAggg                                                                                                             <  2:1015881/150‑1 (MQ=255)
                                          gTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                                                                                                            >  2:397913/1‑150 (MQ=255)
                                           tGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGATTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                                                           >  1:88045/1‑150 (MQ=255)
                                           tGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                                                           <  2:134772/150‑1 (MQ=255)
                                           tGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                                                           <  2:498761/150‑1 (MQ=255)
                                           tGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAACCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGcc                                                                                                           <  2:674986/150‑1 (MQ=255)
                                             ccTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                                                                                                         >  1:849955/1‑150 (MQ=255)
                                                     gAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCt                                                                                                 <  1:1265224/150‑1 (MQ=255)
                                                       ataatTGTCAGGCGGCGATTTTAGCTTGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTAc                                                                                               <  1:947587/150‑1 (MQ=255)
                                                        taatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa                                                                                              <  1:1404620/150‑1 (MQ=255)
                                                         aatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                                                                                             <  1:612850/150‑1 (MQ=255)
                                                         aatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCAGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACaa                                                                                             >  2:334892/1‑150 (MQ=255)
                                                            tGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGATTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGt                                                                                          <  2:88045/150‑1 (MQ=255)
                                                                      cgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                                                                <  2:454766/150‑1 (MQ=255)
                                                                      cgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCggg                                                                                <  1:993938/150‑1 (MQ=255)
                                                                       gATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCgggg                                                                               >  1:648664/1‑150 (MQ=255)
                                                                            tAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGt                                                                          >  1:1350585/1‑150 (MQ=255)
                                                                             aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtt                                                                         >  2:358769/1‑150 (MQ=255)
                                                                             aGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtt                                                                         >  1:1456189/1‑150 (MQ=255)
                                                                              gCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGtta                                                                        >  1:1377060/1‑150 (MQ=255)
                                                                                    tCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGttattttat                                                                  <  1:1380547/150‑1 (MQ=255)
                                                                                      cGCAGATGCTAACGTAATTGTCCGTCAGGGAGTGGCAGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                                                                <  2:1047790/150‑1 (MQ=255)
                                                                                      cGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTg                                                                <  2:586968/150‑1 (MQ=255)
                                                                                              cTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTa                                                        <  2:21231/150‑1 (MQ=255)
                                                                                              atAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCTGGGATGGTTATTTTATTGTCAATTTa                                                        <  2:1325080/149‑1 (MQ=255)
                                                                                                aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATg                                                      <  2:648664/150‑1 (MQ=255)
                                                                                                aaCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATg                                                      >  1:1021262/1‑150 (MQ=255)
                                                                                                   gTAATTTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCg                                                  <  1:1605128/150‑1 (MQ=255)
                                                                                                    tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGTGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCg                                                  <  2:1221149/150‑1 (MQ=255)
                                                                                                    tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCg                                                  >  2:809933/1‑150 (MQ=255)
                                                                                                       ttGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGata                                               <  2:1588981/150‑1 (MQ=255)
                                                                                                        tGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGatat                                              <  1:1385924/150‑1 (MQ=255)
                                                                                                           cagtcagGGAGTGGCTGGTACATATATGTCAAAAGCCCGGAAAATCCGGGCTTCGCAGTACCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGc                                           <  1:908304/150‑1 (MQ=255)
                                                                                                           cagtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGc                                           >  1:766044/1‑150 (MQ=255)
                                                                                                           cagtcagGGAGGGGCTTGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGc                                           <  2:1256778/150‑1 (MQ=255)
                                                                                                            agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCa                                          >  1:299918/1‑150 (MQ=255)
                                                                                                                  ggAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCattgatt                                    <  2:1591087/150‑1 (MQ=255)
                                                                                                                   gAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCattgattg                                   >  2:150286/1‑150 (MQ=255)
                                                                                                                     gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGAc                                 <  1:814880/150‑1 (MQ=255)
                                                                                                                     gTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGAc                                 <  2:113701/150‑1 (MQ=255)
                                                                                                                      tGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACt                                >  1:1712056/1‑150 (MQ=255)
                                                                                                                         cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTtga                             >  2:1253558/1‑150 (MQ=255)
                                                                                                                         cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTtga                             >  1:500032/1‑150 (MQ=255)
                                                                                                                         cTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGGCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTtga                             <  1:617263/150‑1 (MQ=255)
                                                                                                                           ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTtgatg                           <  2:1587257/150‑1 (MQ=255)
                                                                                                                           ggTACATAAATCGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGACTCTATGGGCAACAAAGTCATTATCGGGGTTGATTATTTTTTTGTAAATTTATGAGCGATATGTCAATCATTGACTtgatg                           >  2:142134/1‑150 (MQ=255)
                                                                                                                            gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTtgatga                          <  1:397913/150‑1 (MQ=255)
                                                                                                                             tACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGAc                         >  1:161730/1‑150 (MQ=255)
                                                                                                                              aCATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGAcc                        >  1:1488615/1‑150 (MQ=255)
                                                                                                                               cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCg                       <  2:1517117/150‑1 (MQ=255)
                                                                                                                               cATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTAAGGGCTACAAAGGCATAATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCg                       >  1:977639/1‑150 (MQ=255)
                                                                                                                                 tatatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGTAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCTcc                     >  1:458700/1‑150 (MQ=255)
                                                                                                                                   tatGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCATCTATGGGCTACAAAGTCATTATCGGGGATGGTTAGTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGccat                   >  1:938787/1‑148 (MQ=255)
                                                                                                                                     tGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTAAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCctct                 <  1:1502290/150‑1 (MQ=255)
                                                                                                                                     cctgcgaaGCCCGGATTTTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATTTGGCATTGATTGACTTGATGACCGCCctct                 >  2:1348540/7‑150 (MQ=255)
                                                                                                                                     cctgcgaaGCCCGGATTTTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCctct                 >  1:1537047/7‑150 (MQ=255)
                                                                                                                                      ctgcgaaGCCCGGATTTTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCTCCCTCTg                >  1:259664/6‑150 (MQ=255)
                                                                                                                                        gcgaaGCCCGGATTTTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTa              <  2:259664/147‑1 (MQ=255)
                                                                                                                                        gcgaaGCCCGGATTTTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTa              >  2:1138123/4‑150 (MQ=255)
                                                                                                                                        cAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACCAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGCGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTa              >  2:255567/1‑150 (MQ=255)
                                                                                                                                        cAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCCTTATCGGGGATGCTTATGTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTa              >  1:890253/1‑150 (MQ=255)
                                                                                                                                        cAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTa              >  1:1307830/1‑150 (MQ=255)
                                                                                                                                        cAAAAGCCCGGAAAATCCGGGCTACGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATCGGCTACAAAGCCATGATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCGCTGtt              >  1:75343/1‑149 (MQ=255)
                                                                                                                                            aGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTc          <  2:1162310/150‑1 (MQ=255)
                                                                                                                                              cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGGCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTg        <  2:59798/150‑1 (MQ=255)
                                                                                                                                              cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTTTATGTCTg        >  1:1429176/1‑150 (MQ=255)
                                                                                                                                              cccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTg        >  2:804100/1‑150 (MQ=255)
                                                                                                                                               ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGc       >  2:1632761/1‑150 (MQ=255)
                                                                                                                                               ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGc       <  2:1456189/150‑1 (MQ=255)
                                                                                                                                               ccGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGc       >  2:1183786/1‑150 (MQ=255)
                                                                                                                                                 ggAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGcaa     >  2:646775/1‑150 (MQ=255)
                                                                                                                                                  gAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGcaac    <  2:1712056/150‑1 (MQ=255)
                                                                                                                                                    aaaTCCGGGCTTCGCAGGTCCCTCAAGGTGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGcaacca  <  2:1104264/150‑1 (MQ=255)
                                                                                                                                                    aaaTCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGcaacca  <  2:299918/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTACAAGGATATGCAGTGCTGCTGCTGTTTTTGCCATGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATTATCGGGGATGGTTATTTTATTGTCAATTTATGAGCGATATGGCATTGATTGACTTGATGACCGCCCTCTGTATGTCTGCAACCA  >  NC_000913/3959776‑3960073

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