Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 291 A→T 5.8% intergenic (+36/‑46) thrL → / → thrA thr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132910AT 5.8% 306.7 / 4.7 122intergenic (+36/‑46)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (65/46);  new base T (3/5);  total (70/52)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

TAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  NC_000913/143‑440
                                                                                                                                                     |                                                                                                                                                       
taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCTTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                                                                        <  2:1282877/150‑1 (MQ=255)
taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                                                                        <  1:127395/150‑1 (MQ=255)
taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                                                                        >  1:185838/1‑150 (MQ=255)
 aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACa                                                                                                                                                       >  2:889000/1‑150 (MQ=255)
        gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCCCCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                               >  1:221546/1‑150 (MQ=255)
        gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                               <  2:322803/150‑1 (MQ=255)
          gcACAGACAGATAAATATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTAGCACCACCATCACCATTTCCTCAGGTAACGGTGCGGGCTGACGCGTAAAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCt                                                                                                                                            <  1:1291041/150‑1 (MQ=255)
          gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCCCCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCt                                                                                                                                            <  1:1349953/150‑1 (MQ=255)
           cacaGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtt                                                                                                                                           >  1:120608/1‑150 (MQ=255)
            acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttt                                                                                                                                          <  1:882304/150‑1 (MQ=255)
            acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCttt                                                                                                                                          <  1:400825/150‑1 (MQ=255)
              agacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttt                                                                                                                                        <  2:800754/150‑1 (MQ=255)
                 cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCCTTAGCACCACCATTACCACCACCATCACCATTACCACCCGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCtttttttt                                                                                                                                     <  1:1255852/150‑1 (MQ=255)
                   taTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAGACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTc                                                                                                                                   <  1:1199767/149‑1 (MQ=255)
                   gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTc                                                                                                                                   >  1:1088351/1‑150 (MQ=255)
                    aTAAAAATTACAGAGTACACAACCTCCATGAAACGCCTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCg                                                                                                                                  <  2:1648594/150‑1 (MQ=255)
                    aTAAAAATTACAGAGTACACAACATCCATGAACCGCATTATCACCACCATTACCCCCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCttttttctgtg                                                                                                                                  <  2:916601/150‑6 (MQ=38)
                      aaaaaTTACATAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAc                                                                                                                                <  2:931587/150‑1 (MQ=255)
                       aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAcc                                                                                                                               >  2:337736/1‑150 (MQ=255)
                       aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGAAGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAcc                                                                                                                               >  2:1212455/1‑150 (MQ=255)
                        aaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCa                                                                                                                              >  1:764344/1‑150 (MQ=255)
                             aCAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGt                                                                                                                         >  2:1304255/1‑150 (MQ=255)
                                 agTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTaaaa                                                                                                                     >  1:940211/1‑148 (MQ=255)
                                     cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTAACACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                                 >  2:1127119/1‑150 (MQ=255)
                                         acaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaaca                                                                                                             >  2:173051/1‑150 (MQ=255)
                                             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGTAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCa                                                                                                         >  2:68383/1‑150 (MQ=255)
                                             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                         >  2:629272/1‑150 (MQ=255)
                                              cATGAAACGCATTAGCCCCACCATTACCACCACCATCACCAATACCACAGGTAACGGTGCGGGCTGACGCGGACAGGAAACACAGAAAAAAGCCCGCACCTGACCGGCCGGGCTTGTGTTTGCCACCAAAGGTAACGAGGTAACACCCAt                                                                                                        >  1:345543/1‑150 (MQ=255)
                                              cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                        >  2:1042905/1‑150 (MQ=255)
                                              cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                        >  2:794226/1‑150 (MQ=255)
                                               aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAACCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                       <  2:263193/150‑1 (MQ=255)
                                                tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGCCCAAAGGTAACGCGGTAACACCCATGc                                                                                                      >  1:498655/1‑150 (MQ=255)
                                                   aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                                                                                   <  1:17970/150‑1 (MQ=255)
                                                     cGCATTAGCACCACCATTACCACCCCGATCACCATTACCACAGGTATCTGTGCGGGCTGACGCGTACAGGAAACCCAGAATAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                                                                                 <  2:1046502/150‑1 (MQ=255)
                                                      gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCGGTTTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAgtgt                                                                                                >  2:605782/1‑150 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAAGGAGGTAACAACCATGCGAGTGtt                                                                                               >  1:1302419/1‑150 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               <  2:1524718/150‑1 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               >  2:1540879/1‑150 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               <  1:794226/150‑1 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAGGCGAGTGtt                                                                                               >  1:646182/1‑150 (MQ=255)
                                                       atttcccaccacCAATAACACCACCTTGACCATTACCTCTGGTAAGGGTGCGGTCTGACGCGTACAGGAAACTCAGAAAAACCCCCTGACCTGACATTGCGGTGTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               <  2:1747370/144‑1 (MQ=255)
                                                           aGCACCACGATTACCACCACTATCACCATTACCACAGGTAAGGGTGCGGGCTGACGCGTACTGGAAACACTGAAAAATGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                           <  2:986249/150‑1 (MQ=255)
                                                            gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGt                                                                                          >  1:468653/1‑150 (MQ=255)
                                                               ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                       >  1:52621/1‑150 (MQ=255)
                                                                caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                      >  1:1189211/1‑150 (MQ=255)
                                                                      tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                                                                                <  1:1304255/150‑1 (MQ=255)
                                                                         caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                             >  1:1661706/1‑150 (MQ=255)
                                                                         caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                             >  2:1858619/1‑150 (MQ=255)
                                                                          accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTTACGCGGTAACAACCCTGCGAGTGTTGCAGTTCGGCGGttactca                                                                            >  2:1665685/1‑144 (MQ=255)
                                                                             gcccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACAGTGCGGGCTGTTTTTTTCGCCCAAAGGTAACGAGGTAACACCCATGCGAGTGTTGAAGTTCGGCGGTACCTCAGTg                                                                         >  1:1561072/6‑150 (MQ=255)
                                                                             gcccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                                                                         >  1:1096718/6‑150 (MQ=255)
                                                                             accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                                                                         <  1:385351/150‑1 (MQ=255)
                                                                             accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAACAGTg                                                                         >  2:1789685/1‑150 (MQ=255)
                                                                              cccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGTGTTTTGAAGTGCGGCGGTACATCAGTgg                                                                        >  2:647151/5‑150 (MQ=255)
                                                                              ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                                                                        <  2:185838/150‑1 (MQ=255)
                                                                               ccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                                                                       >  2:778383/4‑150 (MQ=255)
                                                                               catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                                                                       >  2:875795/1‑150 (MQ=255)
                                                                                atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAAGGAGGTAACAACCATGCGTGTGTTGAAGTTCGGCGGTACATCAGTGGCa                                                                      >  2:68171/1‑150 (MQ=255)
                                                                                atcaccatAACCACAGGTAACGGGGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                                                                      <  2:976253/150‑1 (MQ=255)
                                                                                     catTACGACAGGGTACGGTGCGGGCTGACGGGTTGAGGAAACACAGAAAAAAGCCCGCAGCTGCCAGTGAGGGCTTTTTTTTTCGACATAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGCTGCCTCAGTGGCAAATGc                                                                 >  2:103694/1‑150 (MQ=255)
                                                                                     catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCGGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  2:903786/150‑1 (MQ=255)
                                                                                          ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGTGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAATTGCAGAAc                                                            >  2:162898/1‑150 (MQ=255)
                                                                                             caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGtt                                                         >  1:102143/1‑150 (MQ=255)
                                                                                                gTAACGGTTCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      <  1:1224505/150‑1 (MQ=255)
                                                                                                 tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGTTAAAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                     >  1:1841539/1‑150 (MQ=255)
                                                                                                  aaCGGTGCGGGCTGACGCGTAAAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTg                                                    <  1:227582/150‑1 (MQ=255)
                                                                                                  aaCGGGGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTg                                                    <  2:746515/150‑1 (MQ=255)
                                                                                                   aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGc                                                   <  1:540336/150‑1 (MQ=255)
                                                                                                    cGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCg                                                  <  2:1509070/150‑1 (MQ=255)
                                                                                                      gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCTAATGCAGAACGTTTTCTGCgtg                                                <  1:831456/150‑1 (MQ=255)
                                                                                                      gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtg                                                <  1:230020/150‑1 (MQ=255)
                                                                                                          gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                            <  2:1245147/150‑1 (MQ=255)
                                                                                                             cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGGCGa                                         >  1:1566791/1‑150 (MQ=255)
                                                                                                             cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                         <  1:173051/150‑1 (MQ=255)
                                                                                                               gACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                        >  1:1023890/1‑150 (MQ=255)
                                                                                                                  gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTc                                    <  1:1326494/150‑1 (MQ=255)
                                                                                                                    gTACAGGAAACACAGAAAAAATCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACTAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                  <  1:1789309/150‑1 (MQ=255)
                                                                                                                     tACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                 <  1:1854253/150‑1 (MQ=255)
                                                                                                                       cAGGAAACACAGAAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                >  2:1227994/1‑150 (MQ=255)
                                                                                                                       cAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAAAGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATAAGTGGCAAATGCAGCACGTTTTCTGCGTGTTGCCGATATTCTGGaa                               >  1:939994/1‑150 (MQ=255)
                                                                                                                       cAGCAAACACATAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                               <  2:1824584/150‑1 (MQ=255)
                                                                                                                          gAAACACAGAAAAAAGCCCGCACCTGTCAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                            <  1:1653642/150‑1 (MQ=255)
                                                                                                                              cacaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATg                        >  2:663383/1‑150 (MQ=255)
                                                                                                                              cacaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAAAGAGGTAACACCCATGCGTGTGTTGAAGTTCGGCGTGACTTCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATTTTGTGGAAAGCAATg                        >  2:336141/1‑150 (MQ=255)
                                                                                                                               acaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGc                       >  2:929723/1‑150 (MQ=255)
                                                                                                                                gaaaaaaaaaGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                       >  1:717114/5‑150 (MQ=39)
                                                                                                                                    aaaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                   >  2:453233/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTCGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggca                 >  2:363427/1‑150 (MQ=255)
                                                                                                                                        aaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg               <  2:102143/150‑1 (MQ=255)
                                                                                                                                         aGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggg              >  2:1625000/1‑150 (MQ=255)
                                                                                                                                           cccGCACCTGACAGTGCGGGCTTTTTTTTTCGCCCAAAGGTAACGAGGTAACAACCATGCGTGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggc            >  2:708835/1‑150 (MQ=255)
                                                                                                                                           cccGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggccggggc            >  2:214490/1‑150 (MQ=255)
                                                                                                                                           cccGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggc            >  1:1151506/1‑150 (MQ=255)
                                                                                                                                             cGCACCTGGCAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                        <  2:580333/104‑1 (MQ=255)
                                                                                                                                             cGCACCTGGCAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                        >  1:580333/1‑104 (MQ=255)
                                                                                                                                             cGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATAAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag          >  1:986301/1‑150 (MQ=255)
                                                                                                                                                 ccTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg      >  1:1005940/1‑150 (MQ=255)
                                                                                                                                                   tGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc    >  2:1424862/1‑150 (MQ=255)
                                                                                                                                                    gACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa   >  2:74465/1‑150 (MQ=255)
                                                                                                                                                     aCAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGAGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAc  <  1:1127119/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  NC_000913/143‑440

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