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.

TAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  NC_000913/143‑440
                                                                                                                                                     |                                                                                                                                                       
taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCTTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGAc                                                                                                                                                        <  2:1282877/150‑1 (MQ=255)
taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGAc                                                                                                                                                        <  1:127395/150‑1 (MQ=255)
taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGAc                                                                                                                                                        >  1:185838/1‑150 (MQ=255)
 aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACa                                                                                                                                                       >  2:889000/1‑150 (MQ=255)
        gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCCCCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                               >  1:221546/1‑150 (MQ=255)
        gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                               <  2:322803/150‑1 (MQ=255)
          gcACAGACAGATAAATATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTAGCACCACCATCACCATTTCCTCAGGTAACGGTGCGGGCTGACGCGTAAAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑t                                                                                                                                            <  1:1291041/150‑1 (MQ=255)
          gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCCCCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑t                                                                                                                                            <  1:1349953/150‑1 (MQ=255)
           cacaGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tt                                                                                                                                           >  1:120608/1‑150 (MQ=255)
            acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttt                                                                                                                                          <  1:882304/150‑1 (MQ=255)
            acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttt                                                                                                                                          <  1:400825/150‑1 (MQ=255)
              agacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttt                                                                                                                                        <  2:800754/150‑1 (MQ=255)
                 cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCCTTAGCACCACCATTACCACCACCATCACCATTACCACCCGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑tttttttt                                                                                                                                     <  1:1255852/150‑1 (MQ=255)
                   taTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAGACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTc                                                                                                                                   <  1:1199767/149‑1 (MQ=255)
                   gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTc                                                                                                                                   >  1:1088351/1‑150 (MQ=255)
                    aTAAAAATTACAGAGTACACAACCTCCATGAAACGCCTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCg                                                                                                                                  <  2:1648594/150‑1 (MQ=255)
                    aTAAAAATTACAGAGTACACAACATCCATGAACCGCATTATCACCACCATTACCCCCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑ttttttctgtg                                                                                                                                  <  2:916601/150‑6 (MQ=38)
                      aaaaaTTACATAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAc                                                                                                                                <  2:931587/150‑1 (MQ=255)
                       aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAcc                                                                                                                               >  2:337736/1‑150 (MQ=255)
                       aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGAAGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAcc                                                                                                                               >  2:1212455/1‑150 (MQ=255)
                        aaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCa                                                                                                                              >  1:764344/1‑150 (MQ=255)
                             aCAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGt                                                                                                                         >  2:1304255/1‑150 (MQ=255)
                                 agTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTaaaa                                                                                                                     >  1:940211/1‑148 (MQ=255)
                                     cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTAACACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                                 >  2:1127119/1‑150 (MQ=255)
                                         acaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaaca                                                                                                             >  2:173051/1‑150 (MQ=255)
                                             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGTAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCa                                                                                                         >  2:68383/1‑150 (MQ=255)
                                             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                         >  2:629272/1‑150 (MQ=255)
                                              cATGAAACGCATTAGCCCCACCATTACCACCACCATCACCAATACCACAGGTAACGGTGCGGGCTGACGCGGACAGGAAACACAG‑AAAAAAGCCCGCACCTGACC‑GGCCGGGC‑TTGTGTTTGCCACCAAAGGTAACGAGGTAACACCCAt                                                                                                        >  1:345543/1‑150 (MQ=255)
                                              cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                        >  2:1042905/1‑150 (MQ=255)
                                              cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                        >  2:794226/1‑150 (MQ=255)
                                               aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAACCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                       <  2:263193/150‑1 (MQ=255)
                                                tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGCCCAAAGGTAACGCGGTAACACCCATGc                                                                                                      >  1:498655/1‑150 (MQ=255)
                                                   aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                                                                                   <  1:17970/150‑1 (MQ=255)
                                                     cGCATTAGCACCACCATTACCACCCCGATCACCATTACCACAGGTATCTGTGCGGGCTGACGCGTACAGGAAACCCAG‑AATAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                                                                                 <  2:1046502/150‑1 (MQ=255)
                                                      gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑GGTTTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAgtgt                                                                                                >  2:605782/1‑150 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAAGGAGGTAACAACCATGCGAGTGtt                                                                                               >  1:1302419/1‑150 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               <  2:1524718/150‑1 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               >  2:1540879/1‑150 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               <  1:794226/150‑1 (MQ=255)
                                                       cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAGGCGAGTGtt                                                                                               >  1:646182/1‑150 (MQ=255)
                                                       atttcccaccacCAATAACACCACCTTGACCATTACCTCTGGTAAGGGTGCGGTCTGACGCGTACAGGAAACTCAG‑AAAAACCCCCTGACCTGACA‑TTGCGGTG‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               <  2:1747370/144‑1 (MQ=255)
                                                           aGCACCACGATTACCACCACTATCACCATTACCACAGGTAAGGGTGCGGGCTGACGCGTACTGGAAACACTG‑AAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                           <  2:986249/150‑1 (MQ=255)
                                                            gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGt                                                                                          >  1:468653/1‑150 (MQ=255)
                                                               ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                       >  1:52621/1‑150 (MQ=255)
                                                                caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                      >  1:1189211/1‑150 (MQ=255)
                                                                      tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                                                                                <  1:1304255/150‑1 (MQ=255)
                                                                         caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                             >  1:1661706/1‑150 (MQ=255)
                                                                         caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                             >  2:1858619/1‑150 (MQ=255)
                                                                          accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTTACGCGGTAACAACCCTGCGAGTGTTGCAGTTCGGCGGttactca                                                                            >  2:1665685/1‑144 (MQ=255)
                                                                             gcccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TGTTTTTTTCGCCCAAAGGTAACGAGGTAACACCCATGCGAGTGTTGAAGTTCGGCGGTACCTCAGTg                                                                         >  1:1561072/6‑150 (MQ=255)
                                                                             gcccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                                                                         >  1:1096718/6‑150 (MQ=255)
                                                                             accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg                                                                         <  1:385351/150‑1 (MQ=255)
                                                                             accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAACAGTg                                                                         >  2:1789685/1‑150 (MQ=255)
                                                                              cccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGTGTTTTGAAGTGCGGCGGTACATCAGTgg                                                                        >  2:647151/5‑150 (MQ=255)
                                                                              ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                                                                        <  2:185838/150‑1 (MQ=255)
                                                                               ccgcaccgtTACCTGTGGTAATGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                                                                       >  2:778383/4‑150 (MQ=255)
                                                                               catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                                                                       >  2:875795/1‑150 (MQ=255)
                                                                                atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAAGGAGGTAACAACCATGCGTGTGTTGAAGTTCGGCGGTACATCAGTGGCa                                                                      >  2:68171/1‑150 (MQ=255)
                                                                                atcaccatAACCACAGGTAACGGGGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                                                                      <  2:976253/150‑1 (MQ=255)
                                                                                     catTACGACAGGGTACGGTGCGGGCTGACGGGTTGAGGAAACACAG‑AAAAAAGCCCGCAGCTGCCA‑GTGAGGGC‑TTTTTTTTTCGACATAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGCTGCCTCAGTGGCAAATGc                                                                 >  2:103694/1‑150 (MQ=255)
                                                                                     catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCGGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  2:903786/150‑1 (MQ=255)
                                                                                          ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGTGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAATTGCAGAAc                                                            >  2:162898/1‑150 (MQ=255)
                                                                                             caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGtt                                                         >  1:102143/1‑150 (MQ=255)
                                                                                                gTAACGGTTCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      <  1:1224505/150‑1 (MQ=255)
                                                                                                 tAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACGAGGTTAAAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                     >  1:1841539/1‑150 (MQ=255)
                                                                                                  aaCGGTGCGGGCTGACGCGTAAAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTg                                                    <  1:227582/150‑1 (MQ=255)
                                                                                                  aaCGGGGCGGGCTGCCGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTg                                                    <  2:746515/150‑1 (MQ=255)
                                                                                                   aCGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGc                                                   <  1:540336/150‑1 (MQ=255)
                                                                                                    cGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCg                                                  <  2:1509070/150‑1 (MQ=255)
                                                                                                      gTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCTAATGCAGAACGTTTTCTGCgtg                                                <  1:831456/150‑1 (MQ=255)
                                                                                                      gTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtg                                                <  1:230020/150‑1 (MQ=255)
                                                                                                          gggCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                            <  2:1245147/150‑1 (MQ=255)
                                                                                                             cTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGGCGa                                         >  1:1566791/1‑150 (MQ=255)
                                                                                                             cTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                         <  1:173051/150‑1 (MQ=255)
                                                                                                               gACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGat                                        >  1:1023890/1‑150 (MQ=255)
                                                                                                                  gcgTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTc                                    <  1:1326494/150‑1 (MQ=255)
                                                                                                                    gTACAGGAAACACAG‑AAAAAATCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACTAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                  <  1:1789309/150‑1 (MQ=255)
                                                                                                                     tACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                 <  1:1854253/150‑1 (MQ=255)
                                                                                                                       cAGGAAACACAGAAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                >  2:1227994/1‑150 (MQ=255)
                                                                                                                       cAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTAGACCAAAGGTAAAGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATAAGTGGCAAATGCAGCACGTTTTCTGCGTGTTGCCGATATTCTGGaa                               >  1:939994/1‑150 (MQ=255)
                                                                                                                       cAGCAAACACAT‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                               <  2:1824584/150‑1 (MQ=255)
                                                                                                                          gAAACACAG‑AAAAAAGCCCGCACCTGTCA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                            <  1:1653642/150‑1 (MQ=255)
                                                                                                                              cacaG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATg                        >  2:663383/1‑150 (MQ=255)
                                                                                                                              cacaG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAAAGAGGTAACACCCATGCGTGTGTTGAAGTTCGGCGTGACTTCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATTTTGTGGAAAGCAATg                        >  2:336141/1‑150 (MQ=255)
                                                                                                                               acaG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGc                       >  2:929723/1‑150 (MQ=255)
                                                                                                                                gaaaaaaaaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                       >  1:717114/5‑150 (MQ=39)
                                                                                                                                    aaaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                   >  2:453233/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTCGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggca                 >  2:363427/1‑150 (MQ=255)
                                                                                                                                        aaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg               <  2:102143/150‑1 (MQ=255)
                                                                                                                                         aGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggg              >  2:1625000/1‑150 (MQ=255)
                                                                                                                                           cccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGCCCAAAGGTAACGAGGTAACAACCATGCGTGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggc            >  2:708835/1‑150 (MQ=255)
                                                                                                                                           cccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggccggggc            >  2:214490/1‑150 (MQ=255)
                                                                                                                                           cccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggc            >  1:1151506/1‑150 (MQ=255)
                                                                                                                                             cGCACCTGGCA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                        <  2:580333/104‑1 (MQ=255)
                                                                                                                                             cGCACCTGGCA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                        >  1:580333/1‑104 (MQ=255)
                                                                                                                                             cGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATAAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag          >  1:986301/1‑150 (MQ=255)
                                                                                                                                                 ccTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg      >  1:1005940/1‑150 (MQ=255)
                                                                                                                                                   tGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc    >  2:1424862/1‑150 (MQ=255)
                                                                                                                                                    gACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa   >  2:74465/1‑150 (MQ=255)
                                                                                                                                                     aCA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGAGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAc  <  1:1127119/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  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: ‑