Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,880,963 C→A R29I (AGA→ATA)  casD ← type I‑E CRISPR system Cascade subunit CasD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,880,9630CA99.1% 398.1 / ‑5.7 113R29I (AGA→ATA) casDtype I‑E CRISPR system Cascade subunit CasD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (63/49);  minor base T (1/0);  total (64/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 115 total aligned reads displayed.

TCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAA  >  NC_000913/2880831‑2881096
                                                                                                                                    |                                                                                                                                     
tCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGt                                                                                                                                >  2:800518/1‑140 (MQ=255)
  aaGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAg                                                                                                                              >  2:541970/1‑140 (MQ=255)
    gAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGt                                                                                                                            >  1:129485/1‑140 (MQ=255)
        gAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTcgcg                                                                                                                        <  1:1169426/140‑1 (MQ=255)
         aGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGt                                                                                                                                <  2:277021/131‑1 (MQ=255)
         aGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGt                                                                                                                                >  1:277021/1‑131 (MQ=255)
         aGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGt                                                                                                                       >  1:235042/1‑140 (MQ=255)
           ttCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTc                                                                                                                     >  2:1118062/1‑140 (MQ=255)
               tCGCAGCACACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTc                                                                                                                 >  1:681059/1‑140 (MQ=255)
                        aCTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCgg                                                                                                        <  1:918368/140‑1 (MQ=255)
                        aCTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCgg                                                                                                        >  2:945094/1‑140 (MQ=255)
                         cTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGc                                                                                                       <  2:61663/140‑1 (MQ=255)
                         cTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGc                                                                                                       <  2:1332413/140‑1 (MQ=255)
                          tGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCt                                                                                                      <  1:199480/140‑1 (MQ=255)
                          tGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCt                                                                                                      <  2:1404234/140‑1 (MQ=255)
                          tGCAAATTGCACA‑‑CTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGc                                                                                                    >  1:493841/1‑140 (MQ=255)
                            cAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGc                                                                                                    >  2:918420/1‑140 (MQ=255)
                               aTTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGcccc                                                                                                 <  2:1026974/140‑1 (MQ=255)
                                    acacTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGc                                                                                            >  1:1117452/1‑140 (MQ=255)
                                        tctctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGcccc                                                                                                 <  2:795936/131‑1 (MQ=255)
                                        tctctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGcccc                                                                                                 >  1:795936/1‑131 (MQ=255)
                                         ctctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCa                                                                                       <  1:1242444/140‑1 (MQ=255)
                                         ctctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCa                                                                                       <  2:235042/140‑1 (MQ=255)
                                         ctctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCa                                                                                       >  1:434743/1‑140 (MQ=255)
                                          tctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCAt                                                                                      <  2:788956/140‑1 (MQ=255)
                                          tctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCAt                                                                                      <  1:32816/140‑1 (MQ=255)
                                                 aCGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCa                                                                               >  1:1151485/1‑140 (MQ=255)
                                                  cGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAg                                                                              >  2:307325/1‑140 (MQ=255)
                                                  cGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAg                                                                              <  1:425559/140‑1 (MQ=255)
                                                  cGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAg                                                                              >  2:28033/1‑140 (MQ=255)
                                                  cGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAg                                                                              >  1:450921/1‑140 (MQ=255)
                                                  cGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAg                                                                              >  1:574471/1‑140 (MQ=255)
                                                  cGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAg                                                                              >  2:425559/1‑140 (MQ=255)
                                                  cGCCTGTAATGAAGAAGAATCATCACGTTGGATCCCAAGACAAGCCCCGAGGAGCCCTAATAACCAGCTTCGCGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGACGGCTGCCCCCAGGCTTGCATTGGCCCAg                                                                              >  2:1002266/1‑140 (MQ=255)
                                                    ccTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCa                                                                            >  1:416675/1‑140 (MQ=255)
                                                    ccTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCa                                                                            >  1:1172923/1‑140 (MQ=255)
                                                    ccTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCa                                                                            <  2:434743/140‑1 (MQ=255)
                                                    ccTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCa                                                                            >  1:656141/1‑140 (MQ=255)
                                                    ccTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCa                                                                            >  1:906675/1‑140 (MQ=255)
                                                      tGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAg                                                                          >  2:727375/1‑140 (MQ=255)
                                                        tAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGcc                                                                        >  1:397777/1‑140 (MQ=255)
                                                        tAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGcc                                                                        >  1:734639/1‑140 (MQ=255)
                                                         aaTGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATACCCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCg                                                                       <  2:1385844/140‑1 (MQ=255)
                                                         aaTGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCg                                                                       <  1:800518/140‑1 (MQ=255)
                                                          aTGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGc                                                                      <  2:416675/140‑1 (MQ=255)
                                                          aTGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGc                                                                      <  1:1118062/140‑1 (MQ=255)
                                                          aTGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGc                                                                      <  2:1453762/140‑1 (MQ=255)
                                                           tGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCa                                                                     <  2:906675/140‑1 (MQ=255)
                                                            gaagaaGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCaa                                                                    <  1:429011/140‑1 (MQ=255)
                                                              agaagTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGa                                                                  >  1:905013/1‑140 (MQ=255)
                                                              agaagTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGa                                                                  <  2:1166998/140‑1 (MQ=255)
                                                                aagTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATc                                                                <  2:129485/140‑1 (MQ=255)
                                                                  gTATCATCACGTTGGATTCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCaa                                                              >  2:1235176/1‑140 (MQ=255)
                                                                  gTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCaa                                                              >  1:177976/1‑140 (MQ=255)
                                                                  gTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCaa                                                              >  1:933385/1‑140 (MQ=255)
                                                                  gTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCaa                                                              >  1:140753/1‑140 (MQ=255)
                                                                    atcatcACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATACCCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAAt                                                            <  2:695966/140‑1 (MQ=255)
                                                                    atcatcACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAAt                                                            <  2:734639/140‑1 (MQ=255)
                                                                    atcatcACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAAt                                                            <  2:1160500/140‑1 (MQ=255)
                                                                    atcatcACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAAt                                                            >  2:90621/1‑140 (MQ=255)
                                                                     tcatcaCGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATa                                                           >  1:84693/1‑140 (MQ=255)
                                                                       atcaCGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAg                                                         >  1:631246/1‑140 (MQ=255)
                                                                            gTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGAtctc                                                    >  1:112785/1‑140 (MQ=255)
                                                                             ttGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCa                                                   >  2:563933/1‑140 (MQ=255)
                                                                             ttGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGGAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCa                                                   >  1:122933/1‑140 (MQ=255)
                                                                              tGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCAt                                                  >  1:338596/1‑140 (MQ=255)
                                                                                 aTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATTTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGtt                                               >  1:585756/1‑140 (MQ=255)
                                                                                   cccAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCa                                             <  2:656141/140‑1 (MQ=255)
                                                                                        gACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCt                                        <  2:112785/140‑1 (MQ=255)
                                                                                        gACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCt                                        <  2:84693/140‑1 (MQ=255)
                                                                                        gACAAGCCCCGAGTAGCCATAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCt                                        <  1:658761/140‑1 (MQ=255)
                                                                                         aCAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTc                                       <  1:727375/140‑1 (MQ=255)
                                                                                          cAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCg                                      >  1:1164623/1‑140 (MQ=255)
                                                                                             gtcccGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCa                                   >  2:470216/3‑140 (MQ=255)
                                                                                             gCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCa                                   >  1:347373/1‑140 (MQ=255)
                                                                                             gCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCa                                   >  2:1055695/1‑140 (MQ=255)
                                                                                             gCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCa                                   >  1:331086/1‑140 (MQ=255)
                                                                                             gCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGCAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCa                                   >  2:307849/1‑140 (MQ=255)
                                                                                              ccccGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCat                                  >  1:498616/1‑140 (MQ=255)
                                                                                              ccccGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCat                                  <  2:450921/140‑1 (MQ=255)
                                                                                                   aGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCAttatta                             >  2:1291996/1‑140 (MQ=255)
                                                                                                   aGTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCAttatta                             >  2:641871/1‑140 (MQ=255)
                                                                                                    gTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTAc                            <  1:1144858/140‑1 (MQ=255)
                                                                                                    gTAGCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTAc                            <  1:28033/140‑1 (MQ=255)
                                                                                                       gCCCTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGaa                         >  1:413011/1‑140 (MQ=255)
                                                                                                          cTAATAACCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAAccc                      <  2:493841/140‑1 (MQ=255)
                                                                                                           taataaCCCGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGTCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCa                     <  2:172669/140‑1 (MQ=255)
                                                                                                                 cccGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGAttt               >  1:853169/1‑140 (MQ=255)
                                                                                                                 cccGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGAttt               >  1:581713/1‑140 (MQ=255)
                                                                                                                   cGCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTa             >  2:417745/1‑140 (MQ=255)
                                                                                                                    gCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTaa            <  2:498616/140‑1 (MQ=255)
                                                                                                                    gCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTaa            <  2:678508/140‑1 (MQ=255)
                                                                                                                    gCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTaa            <  2:933385/140‑1 (MQ=255)
                                                                                                                    gCTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTaa            <  1:1291996/140‑1 (MQ=255)
                                                                                                                     tttCGGGTCGGAAATATTCCGGTAGGTCGCGTCCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCTCAAGATCAAATAAGACCTCATTTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAAc           <  2:1399795/139‑1 (MQ=255)
                                                                                                                     cTTCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAAc           <  1:1055695/140‑1 (MQ=255)
                                                                                                                       tCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTg         >  1:273884/1‑140 (MQ=255)
                                                                                                                       tCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTg         >  1:88675/1‑140 (MQ=255)
                                                                                                                       tCGGGTCGGAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATAAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTg         >  2:300329/1‑140 (MQ=255)
                                                                                                                              ggAAATATTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTaaa  <  1:470216/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
TCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAA  >  NC_000913/2880831‑2881096

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