Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 291 A→T 6.3% 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 6.3% 1004.8 / 42.7 287intergenic (+36/‑46)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (140/129);  new base T (12/6);  total (152/135)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 288 total aligned reads displayed.

ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCA  >  NC_000913/142‑439
                                                                                                                                                     |                                                                                                                                                       
      aaaGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGc                                                                                                                                                  >  1:2498114/3‑150 (MQ=255)
       tcgcgcACAGACATATAAAAATTACATAGTACACAACATCCATGAAACGCATTAGCACCACCCTTCCCACCACCATCACCATTACCACAGGTAACGGTGCGGGCGGACGCTTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCg                                                                                                                                                 <  1:6379439/148‑1 (MQ=255)
        aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                <  1:5488552/150‑1 (MQ=255)
        aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                <  1:4421942/150‑1 (MQ=255)
        aGCGCACAGAAATATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCTTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                                <  1:5878531/150‑1 (MQ=255)
         gcgcACAGACATATAAAAATTACAGACTACACAACATCCATGAAACGCATTAGCACCACCATTCCCACCACCATCACCATTACCACAGGTAAGGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                               <  1:2130463/150‑1 (MQ=255)
         gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                               <  2:5932608/150‑1 (MQ=255)
            cacaGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGG‑C‑tt                                                                                                                                          <  1:1640864/150‑1 (MQ=255)
                   agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑ttttttttt                                                                                                                                   >  1:4536865/1‑150 (MQ=255)
                   agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑ttttttttt                                                                                                                                   >  1:4289894/1‑150 (MQ=255)
                   agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑ttttttttt                                                                                                                                   >  2:5218710/1‑150 (MQ=255)
                   agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑ttttttttt                                                                                                                                   >  2:4989019/1‑150 (MQ=255)
                    gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTc                                                                                                                                  <  1:1011453/150‑1 (MQ=255)
                      tAAAAATTCCAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGa                                                                                                                                <  1:4381403/150‑1 (MQ=255)
                           aTTACAGAGTACACAACATCCCTGAAACGCATTAGCACCCCCCTTACCACCCCCATCACCATTACCACAGCGAACGGGGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGGCCGCACCTGACT‑GTGCGGG‑C‑TTTTTTTTTCGACCaaa                                                                                                                           <  2:2965988/150‑1 (MQ=255)
                           aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCaaa                                                                                                                           <  2:1073545/150‑1 (MQ=255)
                                 gagTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCCAAGGTAAc                                                                                                                     >  1:2422679/1‑150 (MQ=255)
                                   ttACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGa                                                                                                                   <  1:5786263/149‑1 (MQ=255)
                                      aacaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCTTTACCACCGGTAAGGGTGCGGGCTGACGCGTACTGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                                <  1:1227789/149‑1 (MQ=255)
                                       acaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTa                                                                                                               >  2:6626959/1‑150 (MQ=255)
                                            aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                                                                          >  1:6098679/1‑150 (MQ=255)
                                              ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                        >  2:15957/1‑150 (MQ=255)
                                               cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                       <  1:3841669/150‑1 (MQ=255)
                                                aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                      <  1:5655875/150‑1 (MQ=255)
                                                 tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                     >  1:91635/1‑150 (MQ=255)
                                                 tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                     >  1:3282694/1‑150 (MQ=255)
                                                   aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCgc                                                                                                   >  2:5331141/1‑149 (MQ=255)
                                                    accGCATTACCACCACCTTTACACCCCCCTCCACCATTACCCCAGGTAAGGGTGCGGGCTTACGCGTACCGTAAACACAGAAAAAAGCCGGAACCTTACC‑GTGCGGGAT‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                                                                   <  2:3368037/148‑1 (MQ=255)
                                                       gCATTCGCACCACCATTACCCCCCCCATCACCCTTACCCCAGGTAACGGGGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCGGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCACAGGTAACGAGGTAACAACCATGCGAgtgt                                                                                               <  2:1845550/150‑1 (MQ=255)
                                                       gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑ttttttttt                                                                                                                                   <  2:3753402/114‑1 (MQ=255)
                                                        cATTAGCACCCCCATTACCACCACCCTCCCCATTACCCCAGGTAACGGTGCGGGCTGACGCGGACAGGAAACACCGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                              <  2:5352481/150‑1 (MQ=255)
                                                        cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                              <  1:1437123/150‑1 (MQ=255)
                                                         aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                             <  1:3402985/150‑1 (MQ=255)
                                                         aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                             <  1:1490503/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                           <  1:4951542/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTACCAACCATGCGAGTGTTGaa                                                                                           >  1:5295842/1‑150 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                           <  1:4819855/150‑1 (MQ=255)
                                                            aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                          >  1:3018042/1‑150 (MQ=255)
                                                             gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                                                                         >  1:3709483/1‑150 (MQ=255)
                                                              caccacCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAATGTAAAGAGGTAACAACCATGCGAGTGTTGAAttt                                                                                        >  1:1301947/1‑150 (MQ=255)
                                                                ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                      >  1:2021907/1‑150 (MQ=255)
                                                                 caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                     >  1:6738843/1‑150 (MQ=255)
                                                                  accaTTCCCACCACCATCACCATTACCACAGTTAACGGTGCGGGATGACCCTTACAGTAAACACACAAAAAACCCCCAACCTTACA‑TTGCGGG‑A‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                    <  1:2848004/150‑1 (MQ=255)
                                                                  accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                    >  2:744603/1‑150 (MQ=255)
                                                                         ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                                                                             <  2:2499417/150‑1 (MQ=255)
                                                                         ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                                                                             <  1:1939337/150‑1 (MQ=255)
                                                                             caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                         <  2:913033/150‑1 (MQ=255)
                                                                             caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAGGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                         >  2:2833003/1‑150 (MQ=255)
                                                                               ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                                                                       >  2:4846576/1‑150 (MQ=255)
                                                                                catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑ttttttttttc                                                                                                                                 <  2:2116983/91‑3 (MQ=255)
                                                                                catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                                                                      >  2:5825351/1‑150 (MQ=255)
                                                                                  tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa                                                                    >  1:5922196/1‑150 (MQ=255)
                                                                                  tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa                                                                    >  1:3739455/1‑150 (MQ=255)
                                                                                     ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTGGCAAATg                                                                 >  1:2150868/1‑150 (MQ=255)
                                                                                        ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg                                                              >  2:6209750/1‑150 (MQ=255)
                                                                                             acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGAc                                                                                                                               >  1:572700/1‑80 (MQ=255)
                                                                                               aGGTAACGGTGCGGGCTGACGCGAACAGGAAACACAGAAAAAAGCCCGCACCAGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAAAGAGGTAACAACCAAGCGAGTGGTGAAGTTCGGCGGTACATCAGTGGCAAATTCAGAACGttt                                                       >  1:2306221/1‑150 (MQ=255)
                                                                                                 gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                     <  1:188258/150‑1 (MQ=255)
                                                                                                 gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                     <  2:1351075/150‑1 (MQ=255)
                                                                                                       gTGCGGGCTGACGCGTACAGGAAACACAGAATATAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCCGTGGCAAATGCAGAACGTTTTCTGCgtg                                               >  1:3536086/1‑150 (MQ=255)
                                                                                                           gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                           <  1:5825379/150‑1 (MQ=255)
                                                                                                                 aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                     <  1:744603/150‑1 (MQ=255)
                                                                                                                    cgTACAGGAAACACAGAAAAAATCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAAAAACCAAGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAAATCAGAAAGTTTTGTGCGTGTTTCCGATAttat                                  >  1:6511984/1‑148 (MQ=255)
                                                                                                                     gTACAGGAAACCCAGACCACAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAACGGTAACGAGGTAACACCCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAACTGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  2:5767510/150‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAATAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGGCCAAAGGTAAACAGGTAACAACCATGCGAGGGTTTAAGTTCGGCGGTACATCCGTTGCAAAAGCAGAAAGTTTTGTGCGGTTTTCCGAGATTCTg                                 >  1:4991310/1‑150 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  2:6006607/150‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  1:1198229/150‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  1:6154315/150‑1 (MQ=255)
                                                                                                                      tACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                <  2:2284072/150‑1 (MQ=255)
                                                                                                                       aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                               >  2:6090087/1‑150 (MQ=255)
                                                                                                                       aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                               <  1:6286380/150‑1 (MQ=255)
                                                                                                                         aGGAAACACAGCAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                             <  2:4426181/150‑1 (MQ=255)
                                                                                                                           gAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                           >  2:2181571/1‑150 (MQ=255)
                                                                                                                           gAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTATGGAAAGc                           >  1:4132860/1‑150 (MQ=255)
                                                                                                                            aaaCACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑A‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCa                          <  2:5686481/150‑1 (MQ=255)
                                                                                                                               cacaGAACAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATg                       <  2:40028/150‑1 (MQ=255)
                                                                                                                                acaGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCCTGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCacggc                      >  2:4757632/1‑146 (MQ=255)
                                                                                                                                acaGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                      >  2:3099943/1‑150 (MQ=255)
                                                                                                                                 aaaaaaaaaGCCCGCA‑CTGTCAGGTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCCAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGcc                     >  1:6508957/4‑150 (MQ=39)
                                                                                                                                  aGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAGGCca                    >  2:6535678/1‑150 (MQ=255)
                                                                                                                                    aaaaaaGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                  >  2:4875209/1‑150 (MQ=255)
                                                                                                                                    aaaaaaGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                  >  2:3569083/1‑150 (MQ=255)
                                                                                                                                     aaaaaGCCCGCA‑CTGTCAGGTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggc                 >  2:1430145/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggca                >  2:2858154/1‑150 (MQ=255)
                                                                                                                                        aaGCCCGCACCTGACA‑GTGCGGG‑CTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcag               >  2:1227077/1‑150 (MQ=255)
                                                                                                                                        aaGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg              <  2:11343/150‑1 (MQ=255)
                                                                                                                                         aGCCCGCA‑CTGTCAGGTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAAAGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 >  2:4548903/1‑130 (MQ=255)
                                                                                                                                          gCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGTCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg            <  1:2638948/150‑1 (MQ=255)
                                                                                                                                          gCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTAGACCAAAGGTAACCAGGTAACAAACAAGCGAGTGTTTAATTTCGGCGGTACAAAACTGGCAAAATCAGAACGGTTTCTTCGTGTTTACGATATTATGGGAAACAAAGCCAggcagggg            >  1:1804211/1‑150 (MQ=255)
                                                                                                                                          gCCCGCA‑CTGTCAGGTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg            >  2:1597404/1‑150 (MQ=255)
                                                                                                                                           cccGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                <  1:1758511/113‑1 (MQ=255)
                                                                                                                                           cccGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                >  2:1758527/1‑113 (MQ=255)
                                                                                                                                            ccGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGAGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca          >  2:6391427/1‑150 (MQ=255)
                                                                                                                                             cGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag         <  1:6090071/150‑1 (MQ=255)
                                                                                                                                             cGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAGGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGGAAATGCAGAACGTTTTCTTCGTGTTTCCGAGAGTATGGAAATCAAAGCCAcgcaggggccg         >  1:2319838/1‑148 (MQ=255)
                                                                                                                                               cACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt       >  2:3909164/1‑150 (MQ=255)
                                                                                                                                               cACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt       >  2:2227334/1‑150 (MQ=255)
                                                                                                                                               cACCTGACA‑GTGCGGG‑C‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAAACATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGCCAGGCAGGGGCAGGt       >  1:51843/1‑150 (MQ=255)
                                                                                                                                                  cTGACA‑GTGCGGG‑C‑TTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAAATCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGCCAGGCAGGGGCAGGTGGc    >  1:4021297/1‑150 (MQ=255)
                                                                                                                                                    gACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa  >  2:4262832/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGG‑C‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCA  >  NC_000913/142‑439

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