Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,945,020 C→T 8.9% intergenic (+35/+16) csdE → / ← tcdA sulfur acceptor protein CsdE/tRNA threonylcarbamoyladenosine dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,945,0200CT 8.9% 405.4 / 18.7 137intergenic (+35/+16)csdE/tcdAsulfur acceptor protein CsdE/tRNA threonylcarbamoyladenosine dehydratase
Reads supporting (aligned to +/- strand):  ref base C (68/55);  new base T (5/7);  total (73/64)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 141 total aligned reads displayed.

GTCACCACTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCG  >  NC_000913/2944871‑2945169
                                                                                                                                                     |                                                                                                                                                     
       cTGGCATTGTTTGATGAGCTGGGATTACGTGCGCTGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATGAGACGCGAACGTCTATCCGAcc                                                                                                                                               <  2:2194236/150‑1 (MQ=255)
       cTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATGAGACGCGAACGTCTATCCGAcc                                                                                                                                               <  2:2203101/150‑1 (MQ=255)
        tGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATGAGACGCGAACGTCTTATCCGAcc                                                                                                                                               <  2:1493996/150‑1 (MQ=255)
        tGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATGAGACGCGAACGTCTTATCCGAcc                                                                                                                                               <  1:1206439/150‑1 (MQ=255)
         ggCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATGAGACGCGAACGTCTTATCCGAccc                                                                                                                                              <  1:2409087/150‑2 (MQ=255)
         ggCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATGAGACGCGAACGTCTTATCCAAccc                                                                                                                                              <  2:2202493/150‑2 (MQ=255)
         ggCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCt                                                                                                                                              <  2:1768125/150‑1 (MQ=255)
                  tgatgaGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTa                                                                                                                                     <  2:888611/150‑1 (MQ=255)
                  tgatgaGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTa                                                                                                                                     <  2:1553669/150‑1 (MQ=255)
                  tgatgaGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTa                                                                                                                                     <  2:868927/150‑1 (MQ=255)
                  tgatgaGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTa                                                                                                                                     >  2:1899409/1‑150 (MQ=255)
                    atgaGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAATCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAc                                                                                                                                   <  2:1806882/150‑1 (MQ=255)
                    atgaGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAc                                                                                                                                   >  2:2459663/1‑150 (MQ=255)
                     tgaGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAAcc                                                                                                                                  <  2:20686/150‑1 (MQ=255)
                       aGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCt                                                                                                                                >  1:1413474/1‑150 (MQ=255)
                       aGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCt                                                                                                                                >  1:1781298/1‑150 (MQ=255)
                        gCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGTGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTg                                                                                                                               >  2:729273/1‑150 (MQ=255)
                        gCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTg                                                                                                                               >  2:348727/1‑150 (MQ=255)
                             gATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGAcgcg                                                                                                                          >  1:2080338/1‑150 (MQ=255)
                                tACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGcgccg                                                                                                                       <  2:2316280/150‑1 (MQ=255)
                                 aCGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGcgccgc                                                                                                                      >  2:23760/1‑150 (MQ=255)
                                  cGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGGTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGCACCGGGTCGGATAAGACGTTCGCGTCTCATCCGCCCTGATTGTTTAACCCTGACGCGCCGCt                                                                                                                     >  1:2350828/1‑150 (MQ=255)
                                     gcgcAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTc                                                                                                                  >  2:1497178/1‑150 (MQ=255)
                                     gcgcAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTc                                                                                                                  >  1:2056213/1‑150 (MQ=255)
                                       gcAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGc                                                                                                                <  1:1328189/150‑1 (MQ=255)
                                        cAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGcc                                                                                                               <  1:1122481/150‑1 (MQ=255)
                                        cAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGcc                                                                                                               <  1:1503261/150‑1 (MQ=255)
                                        cAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGGTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGcc                                                                                                               <  1:787838/150‑1 (MQ=255)
                                          gCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcat                                                                                                             >  1:1309785/1‑150 (MQ=255)
                                           cTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCATGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcatc                                                                                                            >  2:154590/1‑150 (MQ=255)
                                           cTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcatc                                                                                                            >  1:2314335/1‑150 (MQ=255)
                                           cTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcatc                                                                                                            <  1:1615742/150‑1 (MQ=255)
                                            ttAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcatca                                                                                                           <  2:781337/150‑1 (MQ=255)
                                            ttAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcatca                                                                                                           <  1:1208024/150‑1 (MQ=255)
                                            ttAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcatca                                                                                                           <  1:730581/150‑1 (MQ=255)
                                              aGCGCCTCACGCGGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCcatcata                                                                                                         >  1:945299/1‑149 (MQ=255)
                                              aGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCatcatc                                                                                                         <  1:1313762/150‑1 (MQ=255)
                                               gcgcCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCAtct                                                                                                        <  1:880814/150‑1 (MQ=255)
                                               gcgcCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCAtct                                                                                                        >  2:1791794/1‑150 (MQ=255)
                                                 gcCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCAtcttc                                                                                                      >  1:1435291/1‑150 (MQ=255)
                                                 gcCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGGCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCAtcttc                                                                                                      <  2:1345591/150‑1 (MQ=255)
                                                   cTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATcttctt                                                                                                    >  1:227376/1‑150 (MQ=255)
                                                    tCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTATCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCttcttc                                                                                                   <  2:9245/150‑1 (MQ=255)
                                                     cACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCa                                                                                                  >  1:2299349/1‑150 (MQ=255)
                                                       cGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAgc                                                                                                >  1:1512997/1‑150 (MQ=255)
                                                        gCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAgcg                                                                                               <  1:2229796/150‑1 (MQ=255)
                                                         cAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAgcgc                                                                                              >  1:324333/1‑150 (MQ=255)
                                                           gCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCAt                                                                                            >  1:1496430/1‑150 (MQ=255)
                                                              aGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATgag                                                                                         <  1:545283/150‑1 (MQ=255)
                                                                   ttAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAAcc                                                                                    <  2:138458/150‑1 (MQ=255)
                                                                     aaaTGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGCCCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGc                                                                                  >  2:538393/1‑150 (MQ=255)
                                                                     aaaTGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGc                                                                                  >  1:549904/1‑150 (MQ=255)
                                                                      aaTGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGca                                                                                 >  1:378740/1‑150 (MQ=255)
                                                                      aaTGCATTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGca                                                                                 <  2:2111511/150‑1 (MQ=255)
                                                                          cGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGcaacaa                                                                             >  1:1746996/1‑150 (MQ=255)
                                                                          cGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGcaacaa                                                                             >  1:2265420/1‑150 (MQ=255)
                                                                           gTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACaaa                                                                            >  2:1159585/1‑150 (MQ=255)
                                                                               aGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCa                                                                        >  2:810767/1‑150 (MQ=255)
                                                                                 cGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGCTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCaaa                                                                      <  1:2184633/150‑1 (MQ=255)
                                                                                   aGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAgg                                                                    <  2:1874679/150‑1 (MQ=255)
                                                                                   aGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAgg                                                                    <  1:2291565/150‑1 (MQ=255)
                                                                                    ggCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGt                                                                   >  2:2377403/1‑150 (MQ=255)
                                                                                      cgaaaaTCGATGCGACGAAGCAGCTTTAAGCACAGAAACGCGTCGCATAATACCATCTCGTAACATCCGACCTGATTGTTTAACCCTTACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAACGTgg                                                                 >  2:2325147/6‑150 (MQ=255)
                                                                                      cgaaaaTCGATGAGACGAATCAGGTTTAAGCACAGAAACGCGTCGCATAATACCATCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTgg                                                                 >  2:2333494/6‑150 (MQ=255)
                                                                                         ttATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCgg                                                              <  1:1173354/150‑1 (MQ=255)
                                                                                          tATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGt                                                             <  2:1878606/150‑1 (MQ=255)
                                                                                              gCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCAcc                                                         >  2:412814/1‑150 (MQ=255)
                                                                                                 gCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGCCCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATc                                                      >  2:1575702/1‑150 (MQ=255)
                                                                                                  cgacgaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCg                                                     <  1:23760/150‑1 (MQ=255)
                                                                                                    acgaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGtt                                                   >  1:758472/1‑150 (MQ=255)
                                                                                                     cgaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTg                                                                  <  2:658356/134‑1 (MQ=255)
                                                                                                     cgaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTg                                                                  >  1:658356/1‑134 (MQ=255)
                                                                                                     cgaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCTGTCACCATCGTTg                                                  <  2:88052/150‑1 (MQ=255)
                                                                                                     cgaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTg                                                  <  1:2027124/150‑1 (MQ=255)
                                                                                                     cgaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTg                                                  >  2:334109/1‑150 (MQ=255)
                                                                                                      gaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGc                                                 <  1:1791794/150‑1 (MQ=255)
                                                                                                      gaAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGc                                                 >  1:2211137/1‑150 (MQ=255)
                                                                                                       aaGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGCGGTCACCATCGTTGCcg                                               <  1:633260/150‑1 (MQ=255)
                                                                                                          cAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCcgcg                                             <  1:3577/150‑1 (MQ=255)
                                                                                                                  aGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATcc                                     <  1:1672100/150‑1 (MQ=255)
                                                                                                                  aGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATcc                                     <  2:324333/150‑1 (MQ=255)
                                                                                                                  aGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATcc                                     <  1:1515915/150‑1 (MQ=255)
                                                                                                                  aGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCACCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATcc                                     >  2:1985641/1‑150 (MQ=255)
                                                                                                                    cacaGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCac                                   >  2:240977/1‑149 (MQ=255)
                                                                                                                    cacaGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAg                                   >  2:1188350/1‑150 (MQ=255)
                                                                                                                         aatcaggTCGGATGAGACGCGAACGTCTTATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGcc                              >  2:230185/6‑149 (MQ=38)
                                                                                                                         aaCCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGca                              >  2:1456344/1‑150 (MQ=255)
                                                                                                                           ccGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGcaca                            <  1:1077010/150‑1 (MQ=255)
                                                                                                                            cGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACaa                           >  1:1714237/1‑150 (MQ=255)
                                                                                                                             gggTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAat                          >  2:1635635/1‑150 (MQ=255)
                                                                                                                             gggTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAat                          <  1:937390/150‑1 (MQ=255)
                                                                                                                              ggTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAatc                         <  2:651057/150‑1 (MQ=255)
                                                                                                                               gTCGGATGAGACGCGAACGTCTTATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAatcc                        >  2:627835/1‑150 (MQ=255)
                                                                                                                                tCGGATGAGACGCGAACGTCTTATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAatcca                       >  2:2435438/1‑150 (MQ=255)
                                                                                                                                   gATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAatccatcc                    >  1:1476782/1‑150 (MQ=255)
                                                                                                                                        gACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTc               >  1:1853060/1‑150 (MQ=255)
                                                                                                                                        gACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTc               >  2:313261/1‑150 (MQ=255)
                                                                                                                                          cGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCgg             <  2:2299349/150‑1 (MQ=255)
                                                                                                                                                   ctcATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGc    <  1:2457188/150‑1 (MQ=255)
                                                                                                                                                   ctcATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGGGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGc    <  2:1925244/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCACCACTGGCATTGTTTGATGAGCTGGGATTACGTGCGCAGCTTAGCGCCTCACGCAGCCAGGGGTTAAATGCGTTAAGCGAGGCGATTATCGCTGCGACGAAGCAGGTTTAAGCACAGAACCGGGTCGGATAAGACGTTCGCGTCTCATCCGACCTGATTGTTTAACCCTGACGCGCCGCTTTCGCCATCATCTTCTTCAGCGCATGAGAAACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCG  >  NC_000913/2944871‑2945169

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