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% 1132.9 / 42.5 324intergenic (+36/‑46)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (132/171);  new base T (10/9);  total (143/181)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 327 total aligned reads displayed.

ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  NC_000913/142‑440
                                                                                                                                                     |                                                                                                                                                        
ataGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGa                                                                                                                                                          >  1:5273921/1‑150 (MQ=255)
 taGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                                                                         <  2:711492/150‑1 (MQ=255)
  aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACa                                                                                                                                                        >  1:3418809/1‑150 (MQ=255)
    gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑Gt                                                                                                                                                     >  2:1762443/1‑150 (MQ=255)
        aGCGCACAGACATATAAAAATTACAGACTACACAACACCCATGAAACGCATTAGCACCCCCATTACCACCACCATCACCATTCCCCCCGGTAACGGTGCGGGCTGACGCTTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                 <  1:1181882/150‑1 (MQ=255)
        aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCTTTAGCACCACCATTACCACCACCATCACCATTACCACCGTTATGGGTGGGGGCTGACGCGTACAGTAAACACATAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                 <  1:1758814/150‑1 (MQ=255)
        aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                 <  2:4542820/150‑1 (MQ=255)
         gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                                <  1:3179241/150‑1 (MQ=255)
          cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGc                                                                                                                                               >  1:6413255/1‑150 (MQ=255)
          cgcACACACAGATAAAAATTACAGAGAACACAACAGCCATGAAACGCATTCGCACCACCATGACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGc                                                                                                                                               <  1:1283126/150‑1 (MQ=255)
               agacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑ttttt                                                                                                                                        >  2:5852835/1‑150 (MQ=255)
                  cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATAGCCCGCACCTGACA‑GTGCGGGC‑‑tttttttt                                                                                                                                     >  1:2563341/1‑150 (MQ=255)
                   agaTAAAAATTACAGAGTACACCACATCCCTGACCCGCATTAGCACCACCATTACCACCACCATCACCCTTCCCACAGGTAACGGTGCGGGGTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑ttttttttt                                                                                                                                    <  2:1504142/150‑1 (MQ=255)
                           aTTCCAGAGTACACAACTTCCAAGAAACGCTTTCGCCCCACCATTACCACCACCCTCACCATTCCCCCAGGTAACTGTGCGGGCGGACGCTTACAGGAAACAAAGAAAAACGCCGCAACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCaaa                                                                                                                            <  2:125383/150‑1 (MQ=255)
                           aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCaaa                                                                                                                            <  2:2552130/150‑1 (MQ=255)
                            ttACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCT‑TTTTTTTTTCGACCaaa                                                                                                                            <  2:2898690/150‑1 (MQ=255)
                            ttACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAg                                                                                                                           >  1:1413030/1‑150 (MQ=255)
                               cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTa                                                                                                                        >  2:6185219/1‑150 (MQ=255)
                                agagTACACAAAATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTaa                                                                                                                       <  1:2135066/150‑1 (MQ=255)
                                    tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAg                                                                                                                   >  1:6003531/1‑150 (MQ=255)
                                    tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAg                                                                                                                   >  2:4504943/1‑150 (MQ=255)
                                      cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                                 <  1:2429644/150‑1 (MQ=255)
                                       acaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTa                                                                                                                <  1:52162/150‑1 (MQ=255)
                                         aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAGACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                                                                              >  2:4234910/1‑150 (MQ=255)
                                         aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                                                                              >  2:3877404/1‑150 (MQ=255)
                                           caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGGCA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                            >  1:2840364/1‑150 (MQ=255)
                                           caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                            <  1:5482246/150‑1 (MQ=255)
                                                aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCT‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                        >  1:4041170/1‑150 (MQ=255)
                                                aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                       >  1:6364923/1‑150 (MQ=255)
                                                 tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                      >  1:2114119/1‑150 (MQ=255)
                                                  gAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCT‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                      <  2:2105085/150‑1 (MQ=255)
                                                   aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACC‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCATGCGa                                                                                                    >  2:2383231/1‑150 (MQ=255)
                                                     aCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                                                                                  >  1:4532660/1‑150 (MQ=255)
                                                       gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAATGGTAACGAGGTAACAACCATGCGAgtgt                                                                                                >  2:2020431/1‑150 (MQ=255)
                                                       gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtgt                                                                                                >  2:2603118/1‑150 (MQ=255)
                                                           tAGCACCCCCATTACCACCACCATCACCATTACCCCCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:1961299/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:386959/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:3719797/150‑1 (MQ=255)
                                                              caccacCATTACCACCACCATCACGATTAACTCATGTAACGGTGCGGGCTGACGCGTACTGGAAACACAGAAAAAAGACCGCACCTGACA‑TTGCGTGA‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCATGCGAGTGTTGAAGtt                                                                                         <  2:5614671/150‑1 (MQ=255)
                                                                ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                       >  1:1515101/1‑150 (MQ=255)
                                                                 caccaTTACCACCACCATCACCATTACCACAGGTCACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTACCAACCATGCGCGGGTTGAAGTTcgg                                                                                      >  2:4455047/1‑150 (MQ=255)
                                                                 caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                      >  1:4102956/1‑150 (MQ=255)
                                                                   ccaTTACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg                                                                                    <  1:6123191/150‑1 (MQ=255)
                                                                     atgaccacccccATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                                                                                  <  2:5447887/147‑1 (MQ=255)
                                                                      ttACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                 <  1:4894640/150‑1 (MQ=255)
                                                                      ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                 <  2:2318351/150‑1 (MQ=255)
                                                                      ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                 >  2:4791794/1‑150 (MQ=255)
                                                                      ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                 >  1:3538904/1‑150 (MQ=255)
                                                                          caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                             <  1:3240635/150‑1 (MQ=255)
                                                                           accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAAAACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTAGACCAAAGGTAAAGAGGTAAAAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                                                                            >  2:4993726/1‑150 (MQ=255)
                                                                             caccaacaccatTACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCCGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                          <  2:2472249/150‑1 (MQ=255)
                                                                                 atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                                                                      >  1:6571546/1‑150 (MQ=255)
                                                                                  tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa                                                                     >  2:4282626/1‑150 (MQ=255)
                                                                                    accatTACCACAGGTACCGGGGCGGGCTGACGCGTACAGGACACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTCACGAGGTAACACCCCGGCGAGTCTTGAAGTTCGGCGCGACAGCACTGGcccag                                                                   >  2:5045441/1‑146 (MQ=255)
                                                                                      catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  1:1851175/150‑1 (MQ=255)
                                                                                      cagtACCACAGGTAACGGTGCGGGCTGACGCGCACAGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  2:5708808/147‑1 (MQ=255)
                                                                                       agtACCACAGGTAACGGTGCGGGCTGACGAGTACAGGAAACCCCCCAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa                                                                <  2:6281226/148‑1 (MQ=255)
                                                                                        ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg                                                               >  2:3918877/1‑150 (MQ=255)
                                                                                          aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                             <  1:1632451/150‑1 (MQ=255)
                                                                                           ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCT‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGAGGCAAATGCAGaa                                                             >  1:3053830/1‑150 (MQ=255)
                                                                                            cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                             <  1:3756086/150‑1 (MQ=255)
                                                                                            cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACg                                                           >  1:509909/1‑150 (MQ=255)
                                                                                                ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGtttt                                                       >  1:4329141/1‑150 (MQ=255)
                                                                                                 gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      >  1:758274/1‑150 (MQ=255)
                                                                                                 gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      >  1:1006051/1‑150 (MQ=255)
                                                                                                  tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                     <  1:3978324/150‑1 (MQ=255)
                                                                                                    aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGc                                                   <  2:5199756/150‑1 (MQ=255)
                                                                                                      ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                 <  1:5018015/150‑1 (MQ=255)
                                                                                                        tGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtgt                                               >  1:2363319/1‑150 (MQ=255)
                                                                                                         gCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGtt                                              <  1:365931/150‑1 (MQ=255)
                                                                                                           gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               >  1:1469457/1‑99 (MQ=255)
                                                                                                            ggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                           >  2:3162759/1‑150 (MQ=255)
                                                                                                              cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                         <  2:1413030/150‑1 (MQ=255)
                                                                                                              cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                         <  1:623465/150‑1 (MQ=255)
                                                                                                                gACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGata                                       >  2:2675875/1‑150 (MQ=255)
                                                                                                                 aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                      >  2:798810/1‑150 (MQ=255)
                                                                                                                 aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                      >  2:5671714/1‑150 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                  <  1:4281659/150‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                 <  2:365280/150‑1 (MQ=255)
                                                                                                                        aaGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                               <  2:408665/149‑1 (MQ=255)
                                                                                                                          ggAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCTGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                             <  1:2792416/150‑1 (MQ=255)
                                                                                                                          ggAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTACCAACCATGCGCGTGTTGAAGTTCGGCGGTCCATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGACAg                             >  2:5999107/1‑150 (MQ=255)
                                                                                                                          ggAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                             <  1:958326/150‑1 (MQ=255)
                                                                                                                             aaCACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCaa                          <  1:2603118/150‑1 (MQ=255)
                                                                                                                             aaCACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCaa                          <  1:570250/150‑1 (MQ=255)
                                                                                                                                gaaaaaaaaaGCCCGCA‑CGGTCAGGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCCGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                       >  2:4640070/5‑150 (MQ=255)
                                                                                                                                acaGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                       <  1:4791794/150‑1 (MQ=255)
                                                                                                                                  aGAAAAAAGCCCGCACCTGACA‑GTGCGGTC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                     <  2:3406159/150‑1 (MQ=255)
                                                                                                                                  aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                     <  1:4141586/150‑1 (MQ=255)
                                                                                                                                        aaGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg               <  2:3051667/150‑1 (MQ=255)
                                                                                                                                           cccGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggc            >  2:2223312/1‑150 (MQ=255)
                                                                                                                                            ccGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTTAAGTTCGGCGGTACAACAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca           >  1:5996202/1‑150 (MQ=255)
                                                                                                                                            ccGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcc           >  2:2782047/1‑149 (MQ=255)
                                                                                                                                            ccGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcc           >  1:2682923/1‑149 (MQ=255)
                                                                                                                                            ccGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACCGTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca           >  1:6441777/1‑150 (MQ=255)
                                                                                                                                              gCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcagg         <  1:2979342/150‑1 (MQ=255)
                                                                                                                                                aCCTGACA‑GTGCGGGC‑‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg       >  1:672907/1‑150 (MQ=255)
                                                                                                                                                 ccTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg      <  1:5852624/150‑1 (MQ=255)
                                                                                                                                                  cTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc     <  2:6571546/150‑1 (MQ=255)
                                                                                                                                                    gACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa   <  1:891521/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  NC_000913/142‑440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑