Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 287 Δ1 bp 5.9% intergenic (+32/‑50) 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_0009132870C. 5.9% 374.6 / 17.0 119intergenic (+32/‑50)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (64/47);  new base . (3/4);  total (68/51)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 123 total aligned reads displayed.

CAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGG  >  NC_000913/138‑436
                                                                                                                                                      |                                                                                                                                                       
cAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCAATACCACAGGTAACGGTGCGGGCTGACCCGTACAGGAAACACAG‑AAAAAATCACGCAc                                                                                                                                                         >  2:794852/1‑150 (MQ=255)
 aaTATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCAcc                                                                                                                                                        <  1:1780670/150‑1 (MQ=255)
 aaTATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCAcc                                                                                                                                                        >  1:875022/1‑150 (MQ=255)
     taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCTTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGAc                                                                                                                                                    <  2:1282877/150‑1 (MQ=255)
     taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGAc                                                                                                                                                    <  1:127395/150‑1 (MQ=255)
      aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACa                                                                                                                                                   >  2:889000/1‑150 (MQ=255)
            aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                            <  1:1711300/150‑1 (MQ=255)
             gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCCCCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                           >  1:221546/1‑150 (MQ=255)
               gcACAGACAGATAAATATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTAGCACCACCATCACCATTTCCTCAGGTAACGGTGCGGGCTGACGCGTAAAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑t                                                                                                                                        <  1:1291041/150‑1 (MQ=255)
               gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCCCCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑t                                                                                                                                        <  1:1349953/150‑1 (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)
                          tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGGGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGa                                                                                                                             <  2:1674985/150‑1 (MQ=255)
                           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)
                                                  ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGTAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCa                                                                                                     >  2:68383/1‑150 (MQ=255)
                                                  ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCCCCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                     >  2:629272/1‑150 (MQ=255)
                                                   cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                    >  2:794226/1‑150 (MQ=255)
                                                   cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                    >  2:1042905/1‑150 (MQ=255)
                                                   cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                    >  1:1295098/1‑150 (MQ=255)
                                                    aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAACCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                   <  2:263193/150‑1 (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‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                           <  2:1524718/150‑1 (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)
                                                             attaccaccaccattaccaccaccATCACCATAACCACCGGTAACGGTGCGGGCTGACGCGAACAGGAAACACAG‑AACAACGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                          <  2:267325/150‑1 (MQ=255)
                                                                aGCACCACGATTACCACCACTATCACCATTACCACAGGTAAGGGTGCGGGCTGACGCGTACTGGAAACACTG‑AAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                       <  2:986249/150‑1 (MQ=255)
                                                                aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                       <  1:661403/150‑1 (MQ=255)
                                                                 gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGt                                                                                      >  1:468653/1‑150 (MQ=255)
                                                                    ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                   >  2:1525110/1‑150 (MQ=255)
                                                                    ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                   >  1:52621/1‑150 (MQ=255)
                                                                     caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                  >  1:1189211/1‑150 (MQ=255)
                                                                        caTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcgg                                                                               <  2:1851151/150‑1 (MQ=255)
                                                                           tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTAc                                                                            <  1:1304255/150‑1 (MQ=255)
                                                                            accaccaccATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTGTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGAGGTAcg                                                                           >  1:309212/1‑149 (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‑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)
                                                                                    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)
                                                                                         ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATAAGTGGCAAATg                                                              >  1:1206029/1‑150 (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‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                 >  1:1469422/1‑150 (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)
                                                                                                                    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‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                           >  1:486253/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)
                                                                                                                                         aaaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg               >  2:453233/1‑150 (MQ=255)
                                                                                                                                          aaaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGATTTCTTCGTGTTGCCCATATTCTGGAAAGCAATGCcaggc              >  2:532674/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)
                                                                                                                                                      |                                                                                                                                                       
CAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAG‑AAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGG  >  NC_000913/138‑436

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: ‑