Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 866,056 C→T 5.4% P103P (CCG→CCA moeA ← molybdopterin molybdotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913866,0560CT 5.4% 237.1 / 2.1 74P103P (CCG→CCAmoeAmolybdopterin molybdotransferase
Reads supporting (aligned to +/- strand):  ref base C (32/38);  new base T (2/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

CCCGCCGGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATC  >  NC_000913/865907‑866205
                                                                                                                                                     |                                                                                                                                                     
cccGCCGGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCAcc                                                                                                                                                       >  2:1149177/1‑150 (MQ=255)
 ccgccgGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCg                                                                                                                                                      <  1:58166/150‑1 (MQ=255)
  cgccgGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCgg                                                                                                                                                     >  1:969133/1‑150 (MQ=255)
   gccgGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGgc                                                                                                                                                    <  2:1268317/150‑1 (MQ=255)
       ggAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCAc                                                                                                                                                >  1:1794244/1‑150 (MQ=255)
           aacaacCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGgcg                                                                                                                                                   >  2:1168577/1‑143 (MQ=255)
           aacaacCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGgcg                                                                                                                                                   <  1:1168577/143‑1 (MQ=255)
             caacCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCa                                                                                                                                          <  2:1779080/150‑1 (MQ=255)
              aacCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCActagctgg                                                                                                                                         >  1:330703/1‑143 (MQ=255)
                 cGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCataat                                                                                                                                      <  1:480982/150‑1 (MQ=255)
                  gCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCataata                                                                                                                                     <  1:318107/150‑1 (MQ=255)
                   cACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATAc                                                                                                                                    >  2:1196850/1‑150 (MQ=255)
                   cACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATAc                                                                                                                                    >  2:176720/1‑150 (MQ=255)
                   cACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATAc                                                                                                                                    <  2:1932599/150‑1 (MQ=255)
                       tGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAAt                                                                                                                                >  1:1452870/1‑150 (MQ=255)
                        gCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATg                                                                                                                               >  2:551804/1‑150 (MQ=255)
                         cAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGc                                                                                                                              <  2:43063/150‑1 (MQ=255)
                          agagATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCa                                                                                                                             <  2:1411182/150‑1 (MQ=255)
                          agagATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCa                                                                                                                             <  1:854312/150‑1 (MQ=255)
                           gagaTATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAg                                                                                                                            >  1:360302/1‑150 (MQ=255)
                           gagaTATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCACTAAAACGCACGCCATTGTCCATTTTTTCAGTCTGCTCCTGCATCACCACCGCTTCGCTGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAg                                                                                                                            <  2:1064976/150‑1 (MQ=255)
                            agaTATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAgg                                                                                                                           <  1:1196850/150‑1 (MQ=255)
                              atatCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTa                                                                                                                         <  1:466759/150‑1 (MQ=255)
                                atCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACTGGCGCACCGGTCATAATACGAATGCAgg                                                                                                                           >  2:1825062/1‑146 (MQ=255)
                                atCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACTGGCGCACCGGTCATAATACGAATGCAgg                                                                                                                           >  2:1832540/1‑146 (MQ=255)
                                atCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACTGGCGCACCGGTCATAATACGAATGCAgg                                                                                                                           <  1:1832540/146‑1 (MQ=255)
                                atCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACTGGCGCACCGGTCATAATACGAATGCAgg                                                                                                                           <  1:1825062/146‑1 (MQ=255)
                                       ccGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGc                                                                                                                <  2:146271/150‑1 (MQ=255)
                                           gACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCAtt                                                                                                            <  1:457104/150‑1 (MQ=255)
                                                cGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCa                                                                                                       >  1:790342/1‑150 (MQ=255)
                                                cGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCa                                                                                                       >  1:1385692/1‑150 (MQ=255)
                                                   atatTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATgg                                                                                                    <  2:599002/150‑1 (MQ=255)
                                                      ttttGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTAt                                                                                                 >  1:536082/1‑150 (MQ=255)
                                                         tGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGc                                                                                              <  1:1712513/150‑1 (MQ=255)
                                                           cccGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGAACCCGCAGGCCATTCACCATGGTATGGCTg                                                                                            <  2:1221652/150‑1 (MQ=255)
                                                             cGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGAc                                                                                          >  2:748125/1‑150 (MQ=255)
                                                              gCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGAcc                                                                                         >  2:1712427/1‑150 (MQ=255)
                                                              gCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGAcc                                                                                         >  2:1862758/1‑150 (MQ=255)
                                                              gCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGGTCAGTCTGCTCCTGCATCACCACCGCTTCGCCGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCCCCATGGTATGGCTGAcc                                                                                         >  2:798973/1‑150 (MQ=255)
                                                                tACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCgg                                                                                       <  2:444740/150‑1 (MQ=255)
                                                                         cagcagTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAgg                                                                                 <  1:443407/147‑1 (MQ=255)
                                                                         cagcagTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAgg                                                                                 >  2:443407/1‑147 (MQ=255)
                                                                                   acgcacgcCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAAc                                                                    >  2:1597568/1‑150 (MQ=255)
                                                                                      cacgcCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACggg                                                                 >  2:981510/1‑150 (MQ=255)
                                                                                      cacgcCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACggg                                                                 <  1:304733/150‑1 (MQ=255)
                                                                                       acgcCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGc                                                                <  1:1265323/150‑1 (MQ=255)
                                                                                          ccATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGc                                                             <  2:790342/150‑1 (MQ=255)
                                                                                            aTTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCgg                                                           <  2:969133/150‑1 (MQ=255)
                                                                                                tCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGc                                                       <  2:1859892/150‑1 (MQ=255)
                                                                                                 ccATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGcc                                                      >  1:1906124/1‑150 (MQ=255)
                                                                                                   aTTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCg                                                    >  1:363395/1‑150 (MQ=255)
                                                                                                     ttGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGa                                                  <  1:1540568/150‑1 (MQ=255)
                                                                                                      tGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTAcc                                                                          >  2:491785/1‑125 (MQ=255)
                                                                                                      tGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTAcc                                                                          <  1:491785/125‑1 (MQ=255)
                                                                                                        ttCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGc                                               <  1:1712427/150‑1 (MQ=255)
                                                                                                           aGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAat                                            <  2:536082/150‑1 (MQ=255)
                                                                                                                 cTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGc                                      <  1:748125/150‑1 (MQ=255)
                                                                                                                   ccTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTa                                    >  2:1640336/1‑150 (MQ=255)
                                                                                                                      gCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAAc                                 <  1:1150543/150‑1 (MQ=255)
                                                                                                                      gCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAAc                                 >  1:1546898/1‑150 (MQ=255)
                                                                                                                          caccacCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCAc                             >  1:1037849/1‑150 (MQ=255)
                                                                                                                           accaccGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCAcc                            >  2:902733/1‑150 (MQ=255)
                                                                                                                           accaccGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCAcc                            <  2:1145442/150‑1 (MQ=255)
                                                                                                                               ccGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGt                        <  1:1862758/150‑1 (MQ=255)
                                                                                                                                 gCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCTGGAGGCAATATCGGCTAAACGCACCGCGTAg                      <  2:1452870/150‑1 (MQ=255)
                                                                                                                                 gCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAg                      >  1:1406118/1‑150 (MQ=255)
                                                                                                                                    tCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCg                   >  1:559452/1‑150 (MQ=255)
                                                                                                                                        aGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGc                     <  1:1480100/144‑1 (MQ=255)
                                                                                                                                        aGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGc                     >  2:1480100/1‑144 (MQ=255)
                                                                                                                                           cTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTg            >  1:1172424/1‑150 (MQ=255)
                                                                                                                                              ccGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCgg         >  1:1319562/1‑150 (MQ=255)
                                                                                                                                              ccGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGACCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCgg         <  2:626132/150‑1 (MQ=255)
                                                                                                                                                  cACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGtt     <  2:60412/150‑1 (MQ=255)
                                                                                                                                                  cACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGtt     >  1:60412/1‑150 (MQ=255)
                                                                                                                                                   aCCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTa    <  1:551804/150‑1 (MQ=255)
                                                                                                                                                     cGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATc  <  2:1906124/150‑1 (MQ=255)
                                                                                                                                                     cGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATc  >  2:1343150/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCGCCGGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCACGCCATTGTCCATTTGTTCAGTCTGCTCCTGCATCACCACCGCTTCGCAGCCTTCCGGCACCGGCGCACCGGTCATAATACGAATGCAGGTACCCGCAGGCCATTCACCATGGTATGGCTGACCGGCAAAGGATTTACCGGCAACGGGCAGCGGTTGCCCGGAGGCAATATCGGCTAAACGCACCGCGTAGCCGTCCATTGCGGAGTTATC  >  NC_000913/865907‑866205

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