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.

ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCA  >  NC_000913/142‑439
                                                                                                                                                     |                                                                                                                                                       
      aaaGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGc                                                                                                                                                  >  1:2498114/3‑150 (MQ=255)
       tcgcgcACAGACATATAAAAATTACATAGTACACAACATCCATGAAACGCATTAGCACCACCCTTCCCACCACCATCACCATTACCACAGGTAACGGTGCGGGCGGACGCTTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCg                                                                                                                                                 <  1:6379439/148‑1 (MQ=255)
        aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                                                                                <  1:5488552/150‑1 (MQ=255)
        aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                                                                                <  1:4421942/150‑1 (MQ=255)
        aGCGCACAGAAATATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCTTACAGTAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                                                                                <  1:5878531/150‑1 (MQ=255)
         gcgcACAGACATATAAAAATTACAGACTACACAACATCCATGAAACGCATTAGCACCACCATTCCCACCACCATCACCATTACCACAGGTAAGGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                               <  1:2130463/150‑1 (MQ=255)
         gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                               <  2:5932608/150‑1 (MQ=255)
            cacaGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCtt                                                                                                                                          <  1:1640864/150‑1 (MQ=255)
                   agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                                                                                   >  1:4536865/1‑150 (MQ=255)
                   agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                                                                                   >  1:4289894/1‑150 (MQ=255)
                   agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                                                                                   >  2:5218710/1‑150 (MQ=255)
                   agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                                                                                   >  2:4989019/1‑150 (MQ=255)
                    gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTc                                                                                                                                  <  1:1011453/150‑1 (MQ=255)
                      tAAAAATTCCAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGa                                                                                                                                <  1:4381403/150‑1 (MQ=255)
                           aTTACAGAGTACACAACATCCCTGAAACGCATTAGCACCCCCCTTACCACCCCCATCACCATTACCACAGCGAACGGGGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGGCCGCACCTGACTGTGCGGGCTTTTTTTTTCGACCaaa                                                                                                                           <  2:2965988/150‑1 (MQ=255)
                           aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCaaa                                                                                                                           <  2:1073545/150‑1 (MQ=255)
                                 gagTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCCAAGGTAAc                                                                                                                     >  1:2422679/1‑150 (MQ=255)
                                   ttACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGa                                                                                                                   <  1:5786263/149‑1 (MQ=255)
                                      aacaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCTTTACCACCGGTAAGGGTGCGGGCTGACGCGTACTGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                                <  1:1227789/149‑1 (MQ=255)
                                       acaacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTa                                                                                                               >  2:6626959/1‑150 (MQ=255)
                                            aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                                                                          >  1:6098679/1‑150 (MQ=255)
                                              ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                        >  2:15957/1‑150 (MQ=255)
                                               cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                       <  1:3841669/150‑1 (MQ=255)
                                                aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                      <  1:5655875/150‑1 (MQ=255)
                                                 tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                     >  1:91635/1‑150 (MQ=255)
                                                 tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                     >  1:3282694/1‑150 (MQ=255)
                                                   aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCgc                                                                                                   >  2:5331141/1‑149 (MQ=255)
                                                    accGCATTACCACCACCTTTACACCCCCCTCCACCATTACCCCAGGTAAGGGTGCGGGCTTACGCGTACCGTAAACACAGAAAAAAGCCGGAACCTTACCGTGCGGGATTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                                                                   <  2:3368037/148‑1 (MQ=255)
                                                       gCATTCGCACCACCATTACCCCCCCCATCACCCTTACCCCAGGTAACGGGGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCGGACAGTGCGGGCTTTTTTTTTCGACCACAGGTAACGAGGTAACAACCATGCGAgtgt                                                                                               <  2:1845550/150‑1 (MQ=255)
                                                       gCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                                                                                   <  2:3753402/114‑1 (MQ=255)
                                                        cATTAGCACCCCCATTACCACCACCCTCCCCATTACCCCAGGTAACGGTGCGGGCTGACGCGGACAGGAAACACCGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                              <  2:5352481/150‑1 (MQ=255)
                                                        cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                                                                              <  1:1437123/150‑1 (MQ=255)
                                                         aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                             <  1:3402985/150‑1 (MQ=255)
                                                         aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                             <  1:1490503/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                           <  1:4951542/150‑1 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTACCAACCATGCGAGTGTTGaa                                                                                           >  1:5295842/1‑150 (MQ=255)
                                                           tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                           <  1:4819855/150‑1 (MQ=255)
                                                            aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                                                                          >  1:3018042/1‑150 (MQ=255)
                                                             gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                                                                         >  1:3709483/1‑150 (MQ=255)
                                                              caccacCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAATGTAAAGAGGTAACAACCATGCGAGTGTTGAAttt                                                                                        >  1:1301947/1‑150 (MQ=255)
                                                                ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                      >  1:2021907/1‑150 (MQ=255)
                                                                 caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                     >  1:6738843/1‑150 (MQ=255)
                                                                  accaTTCCCACCACCATCACCATTACCACAGTTAACGGTGCGGGATGACCCTTACAGTAAACACACAAAAAACCCCCAACCTTACATTGCGGGATTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                    <  1:2848004/150‑1 (MQ=255)
                                                                  accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                    >  2:744603/1‑150 (MQ=255)
                                                                         ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                                                                             <  2:2499417/150‑1 (MQ=255)
                                                                         ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                                                                             <  1:1939337/150‑1 (MQ=255)
                                                                             caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                         <  2:913033/150‑1 (MQ=255)
                                                                             caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAGGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                         >  2:2833003/1‑150 (MQ=255)
                                                                               ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                                                                       >  2:4846576/1‑150 (MQ=255)
                                                                                catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttttc                                                                                                                                 <  2:2116983/91‑3 (MQ=255)
                                                                                catcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGc                                                                      >  2:5825351/1‑150 (MQ=255)
                                                                                  tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa                                                                    >  1:5922196/1‑150 (MQ=255)
                                                                                  tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaa                                                                    >  1:3739455/1‑150 (MQ=255)
                                                                                     ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACATCAGTGGCAAATg                                                                 >  1:2150868/1‑150 (MQ=255)
                                                                                        ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCCAAGGTAACGAGGTAACAACCATGCGAGTGGTGAAGTTCGGCGGTACATCAGTGGCAAATGCAg                                                              >  2:6209750/1‑150 (MQ=255)
                                                                                             acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAc                                                                                                                               >  1:572700/1‑80 (MQ=255)
                                                                                               aGGTAACGGTGCGGGCTGACGCGAACAGGAAACACAGAAAAAAGCCCGCACCAGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAAAGAGGTAACAACCAAGCGAGTGGTGAAGTTCGGCGGTACATCAGTGGCAAATTCAGAACGttt                                                       >  1:2306221/1‑150 (MQ=255)
                                                                                                 gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                     <  1:188258/150‑1 (MQ=255)
                                                                                                 gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                     <  2:1351075/150‑1 (MQ=255)
                                                                                                       gTGCGGGCTGACGCGTACAGGAAACACAGAATATAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCCGTGGCAAATGCAGAACGTTTTCTGCgtg                                               >  1:3536086/1‑150 (MQ=255)
                                                                                                           gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                           <  1:5825379/150‑1 (MQ=255)
                                                                                                                 aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                     <  1:744603/150‑1 (MQ=255)
                                                                                                                    cgTACAGGAAACACAGAAAAAATCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAAAAACCAAGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAAATCAGAAAGTTTTGTGCGTGTTTCCGATAttat                                  >  1:6511984/1‑148 (MQ=255)
                                                                                                                     gTACAGGAAACCCAGACCACAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAACGGTAACGAGGTAACACCCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAACTGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  2:5767510/150‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAATAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGGCCAAAGGTAAACAGGTAACAACCATGCGAGGGTTTAAGTTCGGCGGTACATCCGTTGCAAAAGCAGAAAGTTTTGTGCGGTTTTCCGAGATTCTg                                 >  1:4991310/1‑150 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  2:6006607/150‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  1:1198229/150‑1 (MQ=255)
                                                                                                                     gTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 <  1:6154315/150‑1 (MQ=255)
                                                                                                                      tACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                <  2:2284072/150‑1 (MQ=255)
                                                                                                                       aCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                               >  2:6090087/1‑150 (MQ=255)
                                                                                                                       aCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                               <  1:6286380/150‑1 (MQ=255)
                                                                                                                         aGGAAACACAGCAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaaa                             <  2:4426181/150‑1 (MQ=255)
                                                                                                                           gAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                           >  2:2181571/1‑150 (MQ=255)
                                                                                                                           gAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTATGGAAAGc                           >  1:4132860/1‑150 (MQ=255)
                                                                                                                            aaaCACAGAAAAAAGCCCGCACCTGACAGTGCGGGATTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCa                          <  2:5686481/150‑1 (MQ=255)
                                                                                                                               cacaGAACAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATg                       <  2:40028/150‑1 (MQ=255)
                                                                                                                                acaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCCTGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCacggc                      >  2:4757632/1‑146 (MQ=255)
                                                                                                                                acaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                      >  2:3099943/1‑150 (MQ=255)
                                                                                                                                 aaaaaaaaaGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCCAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGcc                     >  1:6508957/4‑150 (MQ=39)
                                                                                                                                  aGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAGGCca                    >  2:6535678/1‑150 (MQ=255)
                                                                                                                                    aaaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                  >  2:4875209/1‑150 (MQ=255)
                                                                                                                                    aaaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcagg                  >  2:3569083/1‑150 (MQ=255)
                                                                                                                                     aaaaaGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggc                 >  2:1430145/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggca                >  2:2858154/1‑150 (MQ=255)
                                                                                                                                        aaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcag               >  2:1227077/1‑150 (MQ=255)
                                                                                                                                        aaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg              <  2:11343/150‑1 (MQ=255)
                                                                                                                                         aGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAAAGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                                 >  2:4548903/1‑130 (MQ=255)
                                                                                                                                          gCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGTCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg            <  1:2638948/150‑1 (MQ=255)
                                                                                                                                          gCCCGCACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACCAGGTAACAAACAAGCGAGTGTTTAATTTCGGCGGTACAAAACTGGCAAAATCAGAACGGTTTCTTCGTGTTTACGATATTATGGGAAACAAAGCCAggcagggg            >  1:1804211/1‑150 (MQ=255)
                                                                                                                                          gCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg            >  2:1597404/1‑150 (MQ=255)
                                                                                                                                           cccGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                <  1:1758511/113‑1 (MQ=255)
                                                                                                                                           cccGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                >  2:1758527/1‑113 (MQ=255)
                                                                                                                                            ccGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGAGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca          >  2:6391427/1‑150 (MQ=255)
                                                                                                                                             cGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag         <  1:6090071/150‑1 (MQ=255)
                                                                                                                                             cGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAGGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGGAAATGCAGAACGTTTTCTTCGTGTTTCCGAGAGTATGGAAATCAAAGCCAcgcaggggccg         >  1:2319838/1‑148 (MQ=255)
                                                                                                                                               cACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt       >  2:3909164/1‑150 (MQ=255)
                                                                                                                                               cACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt       >  2:2227334/1‑150 (MQ=255)
                                                                                                                                               cACCTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGTAACAAACATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGCCAGGCAGGGGCAGGt       >  1:51843/1‑150 (MQ=255)
                                                                                                                                                  cTGACAGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTAGGCGGTACATCAGTGGCAAAATCAGAACGTTTTATGCGTGTTGCCGATATTATGGAAAGCAATGCCAGGCAGGGGCAGGTGGc    >  1:4021297/1‑150 (MQ=255)
                                                                                                                                                    gACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCa  >  2:4262832/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
ATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCA  >  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: