Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 680,038 C→A 5.1% C177* (TGC→TGA ybeU → DUF1266 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913680,0380CA 5.1% 569.6 / 14.9 158C177* (TGC→TGAybeUDUF1266 family protein
Reads supporting (aligned to +/- strand):  ref base C (75/75);  new base A (4/4);  total (79/79)
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.58e-01

Only 100 of 161 total aligned reads displayed.

CGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTGGCAGCCATTAACTGTCCAGTCTGAATGCCCGG  >  NC_000913/679891‑680186
                                                                                                                                                   |                                                                                                                                                    
cGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCg                                                                                                                                                     <  4:1210349/149‑1 (MQ=255)
 ggTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGa                                                                                                                                                    <  2:264911/149‑1 (MQ=255)
    gtgCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGa                                                                                                                                                 <  1:1180684/149‑1 (MQ=255)
        tCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGacaac                                                                                                                                             >  1:1072954/1‑149 (MQ=255)
        tCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTCAGGAGGTCTCTACTGGCAATCTTCTCAGTGCGAAGacaac                                                                                                                                             >  4:951531/1‑149 (MQ=255)
          aacaaGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAacac                                                                                                                                           <  2:579296/149‑1 (MQ=255)
          aacaaGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAacac                                                                                                                                           >  4:303225/1‑149 (MQ=255)
             aaGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTc                                                                                                                                        >  1:1092295/1‑149 (MQ=255)
                           agaagaGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGgcgc                                                                                                                          <  2:979390/149‑1 (MQ=255)
                            gaagaGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCt                                                                                                                         >  2:1269135/1‑149 (MQ=255)
                                    tttATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGt                                                                                                                 >  3:1132658/1‑149 (MQ=255)
                                    tttATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGt                                                                                                                 >  4:169324/1‑149 (MQ=255)
                                            tGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGac                                                                                                                                                >  2:267422/1‑110 (MQ=255)
                                            tGCAATCGCGCGTCTATGGACGAGCTCACCATTACTATCAAAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATAc                                                                                                         >  1:424213/1‑149 (MQ=255)
                                                   gcgcgTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAgcg                                                                                                  >  1:855966/1‑149 (MQ=255)
                                                         cTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGa                                                                                            <  3:303225/149‑1 (MQ=255)
                                                            tGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTc                                                                                         <  1:1269135/149‑1 (MQ=255)
                                                               aCGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTAGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATAc                                                                                                         <  3:972343/130‑1 (MQ=255)
                                                                   gCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTaa                                                                                                                       >  2:1045618/1‑112 (MQ=255)
                                                                     tCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCg                                                                                >  3:1282892/1‑149 (MQ=255)
                                                                       aCCTTTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGaa                                                                                                                                                   <  4:810175/80‑1 (MQ=255)
                                                                       aCCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGc                                                                                                                                                      >  3:1323216/1‑77 (MQ=255)
                                                                       aCCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGacaaca                                                                                                                                            <  4:448418/87‑1 (MQ=255)
                                                                       aCCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGacaaca                                                                                                                                            >  3:448418/1‑87 (MQ=255)
                                                                       aCCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGaa                                                                              >  1:15748/1‑149 (MQ=255)
                                                                        ccATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGaa                                                                              <  1:1055187/148‑1 (MQ=255)
                                                                        ccATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGaa                                                                              >  2:1055187/1‑148 (MQ=255)
                                                                         cATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGa                                                                                                       >  3:725717/1‑122 (MQ=255)
                                                                         cATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTa                                                                                      <  4:1046500/139‑1 (MQ=255)
                                                                          aTTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGt                                                                                                                 <  1:709923/111‑1 (MQ=255)
                                                                          aTTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGt                                                                                                                 >  2:709923/1‑111 (MQ=255)
                                                                           ttACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGc                                                                          >  3:20537/1‑149 (MQ=255)
                                                                              cTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTgg                                                                       >  2:1157663/1‑149 (MQ=255)
                                                                               tATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGt                                                                                                                 >  1:380315/1‑106 (MQ=255)
                                                                               tATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTa                                                                                                                <  3:770129/107‑1 (MQ=255)
                                                                               tATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTAGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGc                                                                      >  1:230009/1‑149 (MQ=255)
                                                                                aTCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCgg                                                                                              >  4:1341329/1‑124 (MQ=255)
                                                                                aTCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTgg                                                                       <  3:424242/147‑1 (MQ=255)
                                                                                aTCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCa                                                                     >  4:559615/1‑149 (MQ=255)
                                                                                   acaGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCgg                                                                                              <  2:403156/121‑1 (MQ=255)
                                                                                   acaGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCgg                                                                                              >  1:403156/1‑121 (MQ=255)
                                                                                    caGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCgg                                                                 <  4:1132658/149‑1 (MQ=255)
                                                                                    caGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCgg                                                                 <  1:683379/149‑1 (MQ=255)
                                                                                    caGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCgg                                                                 <  3:559615/149‑1 (MQ=255)
                                                                                       cTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGaa                                                              <  4:1282892/149‑1 (MQ=255)
                                                                                          gATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATa                                                                                                     <  3:192352/107‑1 (MQ=255)
                                                                                          gATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATa                                                                                                     >  4:192352/1‑107 (MQ=255)
                                                                                                cTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGa                                                                                                                                                 >  3:1463031/1‑57 (MQ=255)
                                                                                                cTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGa                                                                                                       >  1:86451/1‑99 (MQ=255)
                                                                                                cTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGa                                                                                                       <  2:86451/99‑1 (MQ=255)
                                                                                                 tATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAgg                                                                                                                             <  4:642926/76‑1 (MQ=255)
                                                                                                 tATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAgg                                                                                                                             >  3:642926/1‑76 (MQ=255)
                                                                                                 tATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGttt                                                    >  2:690524/1‑149 (MQ=255)
                                                                                                  aTTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAaga                                                                            >  3:1152439/1‑124 (MQ=255)
                                                                                                  aTTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAaga                                                                            <  4:1152439/124‑1 (MQ=255)
                                                                                                  aTTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTc                                                   <  1:103930/149‑1 (MQ=255)
                                                                                                  aTTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTATCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAg                                                                                           <  1:288576/109‑1 (MQ=255)
                                                                                                   tttttCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGa                                                                                                                               <  2:1414016/72‑1 (MQ=255)
                                                                                                   tttttCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGa                                                                                                                               >  1:1414016/1‑72 (MQ=255)
                                                                                                   tttttCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCt                                                  <  2:15748/149‑1 (MQ=255)
                                                                                                    ttttCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAg                                                             <  4:1187699/137‑1 (MQ=255)
                                                                                                    ttttCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTa                                                 <  3:404200/149‑1 (MQ=255)
                                                                                                     tttCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATATCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTAt                                                <  2:424213/149‑1 (MQ=255)
                                                                                                          gCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACGAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCAc                                           >  2:198740/1‑149 (MQ=255)
                                                                                                          gCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATAc                                                                                                         <  3:1063685/87‑1 (MQ=255)
                                                                                                            cTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGAt                                                         >  2:1239326/1‑133 (MQ=255)
                                                                                                            cTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGAt                                                         <  1:1239326/133‑1 (MQ=255)
                                                                                                             tATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACGAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCt                                        <  1:198740/149‑1 (MQ=255)
                                                                                                               ttCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCgg                                                                 >  4:775023/1‑122 (MQ=255)
                                                                                                               ttCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCgg                                                                 <  3:775023/122‑1 (MQ=255)
                                                                                                                tCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATa                                                                                                          >  4:1403589/1‑80 (MQ=255)
                                                                                                                tCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATa                                                                                                          <  3:1403589/80‑1 (MQ=255)
                                                                                                                tCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTccc                                    >  2:1086879/1‑149 (MQ=255)
                                                                                                                           tctcTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGa                                                                                                       >  2:293388/1‑72 (MQ=255)
                                                                                                                            ctctACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTa                                                 >  4:467781/1‑125 (MQ=255)
                                                                                                                            ctctACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTa                                                 <  3:467781/125‑1 (MQ=255)
                                                                                                                              ctACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCt                                                                         >  1:318257/1‑99 (MQ=255)
                                                                                                                               tACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTGGCAGCCATTAAc                      <  3:1195123/149‑1 (MQ=255)
                                                                                                                                 cgggCAATCTTCTCAGTGCGAAGAAAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGAAGCGATCGTTTCTATGc                                              <  1:297529/120‑1 (MQ=255)
                                                                                                                                 cTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGAAGCGATCGTTTCTATGc                                              >  2:297529/1‑122 (MQ=255)
                                                                                                                                  tGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCgg                                                                 >  2:34926/1‑103 (MQ=255)
                                                                                                                                   ggCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAg                                                             <  1:901403/106‑1 (MQ=255)
                                                                                                                                   ggCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAg                                                             >  2:901403/1‑106 (MQ=255)
                                                                                                                                    gCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAACTACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGtt                                                     <  4:732211/113‑1 (MQ=255)
                                                                                                                                    gCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTGGCAGCCAt                          >  1:1037320/1‑140 (MQ=255)
                                                                                                                                    gCAATCCTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTATCCCTGTACAAAGACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTGGCAGCCAt                          <  2:1037320/140‑1 (MQ=255)
                                                                                                                                     ttgtcttcgCACTGAGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCgg                                                                                              <  2:1294529/68‑1 (MQ=38)
                                                                                                                                     ttgtcttcgCACTGAGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCgg                                                                                              >  1:1294529/4‑71 (MQ=38)
                                                                                                                                     ttgtcttcgCACTGAGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAagag                                                                           >  2:999560/4‑90 (MQ=255)
                                                                                                                                     ttgtcttcgCACTGAGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAagag                                                                           >  1:614266/4‑90 (MQ=255)
                                                                                                                                     ttgtcttcgCACTGAGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAagag                                                                           <  1:999560/87‑1 (MQ=255)
                                                                                                                                     ttgtcttcgCACTGAGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTgg                                 >  1:1466634/4‑132 (MQ=255)
                                                                                                                                     ttgtcttcgCACTGAGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTgg                                 <  2:1466634/129‑1 (MQ=255)
                                                                                                                                     cAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGa                                                                                            >  4:132288/1‑73 (MQ=255)
                                                                                                                                      aaTCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGttt                                                    >  1:803257/1‑112 (MQ=255)
                                                                                                                                       aTCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTGGCAGCCATTaa                       >  3:929926/1‑140 (MQ=255)
                                                                                                                                        tcttctCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGaa                                                              >  1:579650/1‑100 (MQ=255)
                                                                                                                                        tcttctCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGaa                                                              <  2:579650/100‑1 (MQ=255)
                                                                                                                                                 tGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTGGCAGCCATTAACTGTCCAGTCTGAATGccc    <  2:619875/149‑1 (MQ=255)
                                                                                                                                                   cGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTGGCAGCCATTAACTGTCCAGTCTGAATGCCCgg  <  2:855966/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGACGTCTCTACTGGCAATCTTCTCAGTGCGAAGACAACACTTCGCTACGCGAGGCGCTAACCCTGTACAAATACGATAGCGCCGGGAGTCGTATGTTCGAAGAGCTGGCAGCGGGAAGCGATCGTTTCTATGCCACGCTTCCCTGGCAGCCATTAACTGTCCAGTCTGAATGCCCGG  >  NC_000913/679891‑680186

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