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.

ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑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)
    gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGt                                                                                                                                                     >  2:1762443/1‑150 (MQ=255)
        aGCGCACAGACATATAAAAATTACAGACTACACAACACCCATGAAACGCATTAGCACCCCCATTACCACCACCATCACCATTCCCCCCGGTAACGGTGCGGGCTGACGCTTACAGTAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                                                                                 <  1:1181882/150‑1 (MQ=255)
        aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCTTTAGCACCACCATTACCACCACCATCACCATTACCACCGTTATGGGTGGGGGCTGACGCGTACAGTAAACACATAAAAAAGCCCGCACCTGACAGTGCgg                                                                                                                                                 <  1:1758814/150‑1 (MQ=255)
        aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                                                                                 <  2:4542820/150‑1 (MQ=255)
         gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                                <  1:3179241/150‑1 (MQ=255)
          cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGc                                                                                                                                               >  1:6413255/1‑150 (MQ=255)
          cgcACACACAGATAAAAATTACAGAGAACACAACAGCCATGAAACGCATTCGCACCACCATGACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGGAAACACACAAAAAAGCCCGCACCTGACAGTGCGGGc                                                                                                                                               <  1:1283126/150‑1 (MQ=255)
               agacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑ttttt                                                                                                                                        >  2:5852835/1‑150 (MQ=255)
                  cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATAGCCCGCACCTGACAGTGCGGGC‑‑tttttttt                                                                                                                                     >  1:2563341/1‑150 (MQ=255)
                   agaTAAAAATTACAGAGTACACCACATCCCTGACCCGCATTAGCACCACCATTACCACCACCATCACCCTTCCCACAGGTAACGGTGCGGGGTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑ttttttttt                                                                                                                                    <  2:1504142/150‑1 (MQ=255)
                           aTTCCAGAGTACACAACTTCCAAGAAACGCTTTCGCCCCACCATTACCACCACCCTCACCATTCCCCCAGGTAACTGTGCGGGCGGACGCTTACAGGAAACAAAGAAAAACGCCGCAACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCaaa                                                                                                                            <  2:125383/150‑1 (MQ=255)
                           aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCaaa                                                                                                                            <  2:2552130/150‑1 (MQ=255)
                            ttACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCaaa                                                                                                                            <  2:2898690/150‑1 (MQ=255)
                            ttACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAg                                                                                                                           >  1:1413030/1‑150 (MQ=255)
                               cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTa                                                                                                                        >  2:6185219/1‑150 (MQ=255)
                                agagTACACAAAATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTaa                                                                                                                       <  1:2135066/150‑1 (MQ=255)
                                    tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAg                                                                                                                   >  1:6003531/1‑150 (MQ=255)
                                    tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAg                                                                                                                   >  2:4504943/1‑150 (MQ=255)
                                      cacaACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                                 <  1:2429644/150‑1 (MQ=255)
                                       acaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTa                                                                                                                <  1:52162/150‑1 (MQ=255)
                                         aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAGACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                                                                              >  2:4234910/1‑150 (MQ=255)
                                         aacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                                                                              >  2:3877404/1‑150 (MQ=255)
                                           caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGGCAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                            >  1:2840364/1‑150 (MQ=255)
                                           caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                            <  1:5482246/150‑1 (MQ=255)
                                                aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                        >  1:4041170/1‑150 (MQ=255)
                                                aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                       >  1:6364923/1‑150 (MQ=255)
                                                 tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                      >  1:2114119/1‑150 (MQ=255)
                                                  gAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                      <  2:2105085/150‑1 (MQ=255)
                                                   aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACCGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCATGCGa                                                                                                    >  2:2383231/1‑150 (MQ=255)
                                                     aCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgt                                                                                                  >  1:4532660/1‑150 (MQ=255)
                                                       gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAATGGTAACGAGGTAACAACCATGCGAgtgt                                                                                                >  2:2020431/1‑150 (MQ=255)
                                                       gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAgtgt                                                                                                >  2:2603118/1‑150 (MQ=255)
                                                           tAGCACCCCCATTACCACCACCATCACCATTACCCCCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:1961299/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:386959/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:3719797/150‑1 (MQ=255)
                                                              caccacCATTACCACCACCATCACGATTAACTCATGTAACGGTGCGGGCTGACGCGTACTGGAAACACAGAAAAAAGACCGCACCTGACATTGCGTGA‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCATGCGAGTGTTGAAGtt                                                                                         <  2:5614671/150‑1 (MQ=255)
                                                                ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                       >  1:1515101/1‑150 (MQ=255)
                                                                 caccaTTACCACCACCATCACCATTACCACAGGTCACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTACCAACCATGCGCGGGTTGAAGTTcgg                                                                                      >  2:4455047/1‑150 (MQ=255)
                                                                 caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                      >  1:4102956/1‑150 (MQ=255)
                                                                   ccaTTACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg                                                                                    <  1:6123191/150‑1 (MQ=255)
                                                                     atgaccacccccATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGt                                                                                  <  2:5447887/147‑1 (MQ=255)
                                                                      ttACCACCACCATCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                 <  1:4894640/150‑1 (MQ=255)
                                                                      ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                 <  2:2318351/150‑1 (MQ=255)
                                                                      ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                 >  2:4791794/1‑150 (MQ=255)
                                                                      ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTa                                                                                 >  1:3538904/1‑150 (MQ=255)
                                                                          caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                             <  1:3240635/150‑1 (MQ=255)
                                                                           accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAAAACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTAGACCAAAGGTAAAGAGGTAAAAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                                                                            >  2:4993726/1‑150 (MQ=255)
                                                                             caccaacaccatTACCACCGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCCGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                          <  2:2472249/150‑1 (MQ=255)
                                                                                 atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCa                                                                      >  1:6571546/1‑150 (MQ=255)
                                                                                  tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa                                                                     >  2:4282626/1‑150 (MQ=255)
                                                                                    accatTACCACAGGTACCGGGGCGGGCTGACGCGTACAGGACACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTCACGAGGTAACACCCCGGCGAGTCTTGAAGTTCGGCGCGACAGCACTGGcccag                                                                   >  2:5045441/1‑146 (MQ=255)
                                                                                      catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  1:1851175/150‑1 (MQ=255)
                                                                                      cagtACCACAGGTAACGGTGCGGGCTGACGCGCACAGGAAACACACAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  2:5708808/147‑1 (MQ=255)
                                                                                       agtACCACAGGTAACGGTGCGGGCTGACGAGTACAGGAAACCCCCCAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa                                                                <  2:6281226/148‑1 (MQ=255)
                                                                                        ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg                                                               >  2:3918877/1‑150 (MQ=255)
                                                                                          aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                             <  1:1632451/150‑1 (MQ=255)
                                                                                           ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGAGGCAAATGCAGaa                                                             >  1:3053830/1‑150 (MQ=255)
                                                                                            cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                             <  1:3756086/150‑1 (MQ=255)
                                                                                            cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACg                                                           >  1:509909/1‑150 (MQ=255)
                                                                                                ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGtttt                                                       >  1:4329141/1‑150 (MQ=255)
                                                                                                 gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      >  1:758274/1‑150 (MQ=255)
                                                                                                 gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      >  1:1006051/1‑150 (MQ=255)
                                                                                                  tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                     <  1:3978324/150‑1 (MQ=255)
                                                                                                    aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGc                                                   <  2:5199756/150‑1 (MQ=255)
                                                                                                      ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                 <  1:5018015/150‑1 (MQ=255)
                                                                                                        tGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtgt                                               >  1:2363319/1‑150 (MQ=255)
                                                                                                         gCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGtt                                              <  1:365931/150‑1 (MQ=255)
                                                                                                           gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                               >  1:1469457/1‑99 (MQ=255)
                                                                                                            ggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                           >  2:3162759/1‑150 (MQ=255)
                                                                                                              cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                         <  2:1413030/150‑1 (MQ=255)
                                                                                                              cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGa                                         <  1:623465/150‑1 (MQ=255)
                                                                                                                gACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGata                                       >  2:2675875/1‑150 (MQ=255)
                                                                                                                 aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                      >  2:798810/1‑150 (MQ=255)
                                                                                                                 aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                      >  2:5671714/1‑150 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                  <  1:4281659/150‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                 <  2:365280/150‑1 (MQ=255)
                                                                                                                        aaGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                               <  2:408665/149‑1 (MQ=255)
                                                                                                                          ggAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCTGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                             <  1:2792416/150‑1 (MQ=255)
                                                                                                                          ggAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTACCAACCATGCGCGTGTTGAAGTTCGGCGGTCCATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGACAg                             >  2:5999107/1‑150 (MQ=255)
                                                                                                                          ggAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                             <  1:958326/150‑1 (MQ=255)
                                                                                                                             aaCACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCaa                          <  1:2603118/150‑1 (MQ=255)
                                                                                                                             aaCACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCaa                          <  1:570250/150‑1 (MQ=255)
                                                                                                                                gaaaaaaaaaGCCCGCACGGTCAGGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCCGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                       >  2:4640070/5‑150 (MQ=255)
                                                                                                                                acaGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                       <  1:4791794/150‑1 (MQ=255)
                                                                                                                                  aGAAAAAAGCCCGCACCTGACAGTGCGGTC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                     <  2:3406159/150‑1 (MQ=255)
                                                                                                                                  aGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                     <  1:4141586/150‑1 (MQ=255)
                                                                                                                                        aaGCCCGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg               <  2:3051667/150‑1 (MQ=255)
                                                                                                                                           cccGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggc            >  2:2223312/1‑150 (MQ=255)
                                                                                                                                            ccGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTTAAGTTCGGCGGTACAACAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca           >  1:5996202/1‑150 (MQ=255)
                                                                                                                                            ccGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcc           >  2:2782047/1‑149 (MQ=255)
                                                                                                                                            ccGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcc           >  1:2682923/1‑149 (MQ=255)
                                                                                                                                            ccGCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACCGTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca           >  1:6441777/1‑150 (MQ=255)
                                                                                                                                              gCACCTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcagg         <  1:2979342/150‑1 (MQ=255)
                                                                                                                                                aCCTGACAGTGCGGGC‑‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg       >  1:672907/1‑150 (MQ=255)
                                                                                                                                                 ccTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg      <  1:5852624/150‑1 (MQ=255)
                                                                                                                                                  cTGACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGc     <  2:6571546/150‑1 (MQ=255)
                                                                                                                                                    gACAGTGCGGGC‑‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa   <  1:891521/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑‑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: