Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,881,944 G→A H67Y (CAT→TAT)  casC ← type I‑E CRISPR system Cascade subunit CasC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,881,9440GA98.5% 235.9 / ‑4.3 67H67Y (CAT→TAT) casCtype I‑E CRISPR system Cascade subunit CasC
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (31/35);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GTAACCGCATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATGAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCATGTTCATAT  >  NC_000913/2881808‑2882078
                                                                                                                                        |                                                                                                                                      
gTAACCGCATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAAt                                                                                                                                     <  2:1133390/140‑1 (MQ=255)
gTAACCGCATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAAt                                                                                                                                     <  2:193724/140‑1 (MQ=255)
gTAACCGCATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAAt                                                                                                                                     >  2:44856/1‑140 (MQ=255)
  aaCCGCATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATgg                                                                                                                                   >  1:944088/1‑140 (MQ=255)
  aaCCGCATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATgg                                                                                                                                   >  1:96688/1‑140 (MQ=255)
  aaCCGCATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATgg                                                                                                                                   >  1:7550/1‑140 (MQ=255)
      gCATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCt                                                                                                                               <  1:340798/140‑1 (MQ=255)
       cATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTg                                                                                                                              >  2:33772/1‑140 (MQ=255)
          cGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTgaga                                                                                                                           <  2:72612/140‑1 (MQ=255)
          cGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTgaga                                                                                                                           <  2:1807/140‑1 (MQ=255)
               gAAGTCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGa                                                                                                                      <  2:946/140‑1 (MQ=255)
               gAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGa                                                                                                                      <  2:7550/140‑1 (MQ=255)
               gAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGa                                                                                                                      <  2:466507/140‑1 (MQ=255)
               gAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGa                                                                                                                      <  2:918992/140‑1 (MQ=255)
               gAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGa                                                                                                                      <  2:203314/140‑1 (MQ=255)
               gAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGa                                                                                                                      <  2:944088/140‑1 (MQ=255)
               gAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGa                                                                                                                      <  2:96688/140‑1 (MQ=255)
                   tcctttCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCa                                                                                                                  <  1:227946/137‑1 (MQ=255)
                     ttttCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCAcc                                                                                                                <  2:343415/140‑1 (MQ=255)
                                 tcaACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGc                                                                                                    <  1:515686/140‑1 (MQ=255)
                                       gATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGtaata                                                                                              >  2:1182984/1‑140 (MQ=255)
                                        aTTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGtaataa                                                                                             >  2:467284/1‑140 (MQ=255)
                                         tttACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAAc                                                                                            >  1:974174/1‑140 (MQ=255)
                                         tttACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAAc                                                                                            <  1:33772/140‑1 (MQ=255)
                                          ttACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAAcc                                                                                           >  2:1103720/1‑140 (MQ=255)
                                           tACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCa                                                                                          >  2:67814/1‑140 (MQ=255)
                                                   agcagcGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACg                                                                                  >  2:550830/1‑140 (MQ=255)
                                                     cagcGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATGAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCa                                                                                >  2:418049/1‑140 (MQ=255)
                                                        cgcTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCg                                                                             <  2:569990/140‑1 (MQ=255)
                                                        cgcTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCg                                                                             >  2:1062722/1‑140 (MQ=255)
                                                         gcTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGc                                                                            <  2:732658/140‑1 (MQ=255)
                                                         gcTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGc                                                                            <  1:643336/140‑1 (MQ=255)
                                                             aTGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGt                                                                        >  2:1101691/1‑140 (MQ=255)
                                                                  ttATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAg                                                                   >  1:854022/1‑140 (MQ=255)
                                                                     tCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCt                                                                >  1:1200779/1‑140 (MQ=255)
                                                                        atgatTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTg                                                             >  2:814141/1‑140 (MQ=255)
                                                                          gatTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGAc                                                           >  2:874544/1‑140 (MQ=255)
                                                                          gatTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGAc                                                           >  1:794730/1‑140 (MQ=255)
                                                                             tttttGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTg                                                        <  1:814141/140‑1 (MQ=255)
                                                                                      aaaaCGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTa                                               >  2:521772/1‑140 (MQ=255)
                                                                                       aaaCGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTAc                                              <  1:689503/140‑1 (MQ=255)
                                                                                                 cAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTg                                    <  2:794730/140‑1 (MQ=255)
                                                                                                 cAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTg                                    <  2:854022/140‑1 (MQ=255)
                                                                                                   aGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTgcc                                  <  1:1014409/140‑1 (MQ=255)
                                                                                                         tGCCGAAGAACATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGaa                            <  1:521772/140‑1 (MQ=255)
                                                                                                              aagaaCATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAg                       <  1:78018/140‑1 (MQ=255)
                                                                                                                 aaCATCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGt                    <  1:1103720/140‑1 (MQ=255)
                                                                                                                    aTCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCtt                 <  2:573971/140‑1 (MQ=255)
                                                                                                                    aTCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCtt                 <  1:432672/140‑1 (MQ=255)
                                                                                                                     tCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCttt                <  1:874544/140‑1 (MQ=255)
                                                                                                                     tCACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCttt                >  2:789623/1‑140 (MQ=255)
                                                                                                                      cACGTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTGTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTc               >  2:1121018/1‑140 (MQ=255)
                                                                                                                         gTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGc            <  1:538545/140‑1 (MQ=255)
                                                                                                                         gTAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGc            >  2:210150/1‑140 (MQ=255)
                                                                                                                          tAATTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCa           <  1:390881/140‑1 (MQ=255)
                                                                                                                           aaTTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCAt          >  1:295238/1‑140 (MQ=255)
                                                                                                                           aaTTGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCAt          <  2:188778/140‑1 (MQ=255)
                                                                                                                             ttGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCATGt        >  2:614280/1‑140 (MQ=255)
                                                                                                                             ttGTGCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCATGt        >  1:1185624/1‑140 (MQ=255)
                                                                                                                                tgCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCATGTTCa     >  2:653398/1‑140 (MQ=255)
                                                                                                                                tgCAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCATGTTCa     >  1:378811/1‑140 (MQ=255)
                                                                                                                                  cAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCAt          <  2:235445/133‑1 (MQ=255)
                                                                                                                                  cAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCAt          >  1:235445/1‑133 (MQ=255)
                                                                                                                                  cAAGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCATGTTCata   >  2:392351/1‑140 (MQ=255)
                                                                                                                                   aaGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTg             >  2:815561/1‑129 (MQ=255)
                                                                                                                                   aaGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTg             <  1:815561/129‑1 (MQ=255)
                                                                                                                                   aaGATAAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCATGTTCatat  >  2:1033795/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GTAACCGCATCGGCAGAAATCTTTTCGGCTTCATCAACTGATTTACCGGAGAGCAGCGCTAATGTCTTATCGATGATTTTTTGGTCAAAACGTTCACCAAGTTTTTGCCGAAGAACATCACGTAATTGTGCAAGATGAATGGTTCTGAGACTGGATTCACCAATATTTTGTGCGTAATAACCACTTTTACGCATCGCACGTTTAAGGCTTTGACTTGAAATTCTTACTCGTCTTTTGCCGCCGAAAATAGCGTCTTTCTGCATGTTCATAT  >  NC_000913/2881808‑2882078

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