Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 287 Δ1 bp 6.4% intergenic (+32/‑50) thrL → / → thrA thr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132870C. 6.4% 994.3 / 65.2 282intergenic (+32/‑50)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (136/127);  new base . (12/6);  total (149/133)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 289 total aligned reads displayed.

CAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGG  >  NC_000913/138‑436
                                                                                                                                                     |                                                                                                                                                       
cAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCAc                                                                                                                                                         >  2:645465/1‑150 (MQ=255)
 aaTATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCAcc                                                                                                                                                        <  1:1107192/150‑1 (MQ=255)
 aaTATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCAcc                                                                                                                                                        >  2:2004780/1‑150 (MQ=255)
 aaTATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCAcc                                                                                                                                                        >  1:40028/1‑150 (MQ=255)
  atatAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATAAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCt                                                                                                                                                       >  1:1960072/1‑150 (MQ=255)
    ataGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGa                                                                                                                                                     >  1:2724528/1‑150 (MQ=255)
     tcggCATAGCGCACAGACAGATAAAAATTACAGAGTACAAAACAACCATGAAACGCATTAGCACCACCATTACCACCACCAACACCATTACCACCGTTAACGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGAc                                                                                                                                                    <  1:1201568/148‑1 (MQ=255)
      aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACa                                                                                                                                                   >  1:4977591/1‑150 (MQ=255)
        gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGt                                                                                                                                                >  2:4928762/1‑150 (MQ=255)
        gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGt                                                                                                                                                >  2:4797524/1‑150 (MQ=255)
          aTAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGc                                                                                                                                              <  2:5979280/150‑1 (MQ=255)
          aTAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGc                                                                                                                                              <  1:5047316/150‑1 (MQ=255)
           tAGCGCACATACATATAAAAATTACAGAGAACACACCAACCATGAAACGCATTAGCACCACCATAACCACCACCACCACCATTACCACAGGTAACGGTGCGGGCTGACGCTTACAGTAAACACAGAAAAAAGCCCGCACCTGACAGTGCg                                                                                                                                             <  1:2101890/150‑1 (MQ=255)
           tAGCGCACAGACAGATAAAAATTACAGATTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAAGGGTGCGGGCTGACGCGTACAGTAAACAAATAAAAAAGCCCGCACCTGACAGTGCg                                                                                                                                             <  1:1672712/150‑1 (MQ=255)
            aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTCCCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGTAAACACAGAAAAAAGCCCGCACCTGACAGTGCgg                                                                                                                                            <  1:6154671/150‑1 (MQ=255)
             gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                           >  1:5742339/1‑150 (MQ=255)
             gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                           <  2:4474894/150‑1 (MQ=255)
             gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                           <  2:6365277/150‑1 (MQ=255)
             gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCggg                                                                                                                                           >  2:6154315/1‑150 (MQ=255)
                cacaGAAATATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCACCACCATTCCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAAAAAAAAGCCCGCACTGTCAGGTGCGGGCtt                                                                                                                                       <  1:3177784/150‑1 (MQ=255)
                    gacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACATAAAAAAGCCCGCACTGTCAGGTGCGGGCtttttt                                                                                                                                   <  1:5716947/150‑1 (MQ=255)
                       agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                                                                                >  2:5822756/1‑150 (MQ=255)
                       agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACACAAAAAAGCCCGCACCTGACAGTGCGGGCttttttttt                                                                                                                                <  1:728572/150‑1 (MQ=255)
                        gaTAAAAATTACAGCGTCCACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCCCAGGTAACGGTGCGGGCTGACGCGTACCGGAAACACACAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTc                                                                                                                               <  1:193706/150‑1 (MQ=255)
                        gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTc                                                                                                                               <  1:1011453/150‑1 (MQ=255)
                         aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATAACAACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACCGTGCGGGCTTTTTTTTTCg                                                                                                                              >  1:5386544/1‑150 (MQ=255)
                          tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGa                                                                                                                             <  1:4381445/150‑1 (MQ=255)
                          tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGa                                                                                                                             >  2:222286/1‑150 (MQ=255)
                            aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAcc                                                                                                                           >  1:1243933/1‑150 (MQ=255)
                               aTTACAGAGTACACAACATCCCTGAAACGCATTAGCACCCCCCTTACCACCCCCATCACCATTACCACAGCGAACGGGGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGGCCGCACCTGACTGTGCGGGCTTTTTTTTTCGACCaaa                                                                                                                        <  2:2965988/150‑1 (MQ=255)
                               aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCaaa                                                                                                                        <  2:111005/150‑1 (MQ=255)
                                 tACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACAACCACCATCACAATTACCACAGGTAACGGTGAGGGCTGACGCGTACAGGAAACACAGAAAAATGCCCGCACCTGACAGTGCGGGCTTTTTTATAAGGCCAAAgg                                                                                                                      >  1:2308992/1‑150 (MQ=255)
                                                aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                                                                       >  1:4859979/1‑150 (MQ=255)
                                                aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                                                                       >  2:6586278/1‑150 (MQ=255)
                                                aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                                                                       >  1:6098751/1‑150 (MQ=255)
                                                  ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                     >  2:15957/1‑150 (MQ=255)
                                                   cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                    <  1:3841669/150‑1 (MQ=255)
                                                   cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                    <  2:5917407/150‑1 (MQ=255)
                                                   cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                    >  1:3669718/1‑150 (MQ=255)
                                                     tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                  <  1:21985/150‑1 (MQ=255)
                                                       aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCgc                                                                                                >  1:6426319/1‑149 (MQ=255)
                                                           gCATTCGCACCACCATTACCCCCCCCATCACCCTTACCCCAGGTAACGGGGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCGGACAGTGCGGGCTTTTTTTTTCGACCACAGGTAACGAGGTAACAACCATGCGAgtgt                                                                                            <  2:1845550/150‑1 (MQ=255)
                                                             aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATTGCCCGCAGCTGACAGTGCGGGCTTTTTTTTTCGACCAAATGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                          >  1:1487937/1‑150 (MQ=255)
                                                              atAGCACCACCATTACCACCACCACCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACCGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                                                                         <  1:4090107/149‑1 (MQ=255)
                                                               tatcacaacCATTCCCCACACCATCACCATTTCCTCAGGTAACGTTGCGTGCTGTCGCGTACATGAAACACAGAAAAATGCCCGCACCTGACTGTGCTGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                        <  1:4433661/147‑1 (MQ=255)
                                                                   accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                    >  1:4248855/1‑150 (MQ=255)
                                                                   accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                    >  1:913033/1‑150 (MQ=255)
                                                                    ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCAGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTACCGAGGTAACAACCATGCGAGTGTTTAAGTTcg                                                                                   >  2:1172099/1‑150 (MQ=255)
                                                                    ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGCtttt                                                                                                                                     >  1:9956/1‑100 (MQ=38)
                                                                    ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCAGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGCGTGTTGAAGTTcg                                                                                   >  2:5672625/1‑150 (MQ=255)
                                                                     caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                  >  1:2482163/1‑150 (MQ=255)
                                                                     caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGc                                                                                                                                          <  2:6674037/95‑1 (MQ=38)
                                                                     caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACTGTCAGGTGCGGGc                                                                                                                                          >  1:6674037/1‑95 (MQ=38)
                                                                      accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                 >  2:744603/1‑150 (MQ=255)
                                                                             ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                                                                          <  2:2499417/150‑1 (MQ=255)
                                                                              caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                         >  1:5916094/1‑150 (MQ=255)
                                                                               accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                                                                        >  2:2170256/1‑150 (MQ=255)
                                                                               accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGGACATCa                                                                        >  1:2961239/1‑150 (MQ=255)
                                                                                 caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                      <  2:913033/150‑1 (MQ=255)
                                                                                 caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGGACATCAGt                                                                      >  1:1562459/1‑150 (MQ=255)
                                                                                 caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAGGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                      >  2:2833003/1‑150 (MQ=255)
                                                                                          catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                        >  1:2145349/1‑123 (MQ=255)
                                                                                          catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                             >  2:3354048/1‑150 (MQ=255)
                                                                                           atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                             >  2:5345985/1‑150 (MQ=255)
                                                                                               ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATTGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAATGGTAACGAGGTAACAACCATGGGAGTGTTGAAGTTCGGCGGTTCATCAGTGGCAAATGCAGaag                                                        >  1:5195068/1‑149 (MQ=255)
                                                                                               ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc                                                        <  2:3438905/150‑1 (MQ=255)
                                                                                                cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGAAGAACg                                                       >  1:4306085/1‑150 (MQ=255)
                                                                                                 acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGAc                                                                                                                            <  2:572700/80‑1 (MQ=255)
                                                                                                     gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                  <  1:188258/150‑1 (MQ=255)
                                                                                                     gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                  <  2:4582815/150‑1 (MQ=255)
                                                                                                     gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCAGACAGTGCGGGCATTTTTTTTCGACCAAAGGTAACGAGGTAACCACCAAGCGCGAGTTGAAGTTCGGCGGTACATCAGTGGCAAATTCCGAACGtttta                                                  >  1:3470938/1‑149 (MQ=255)
                                                                                                      tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                 <  2:3669718/150‑1 (MQ=255)
                                                                                                          ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                             <  1:3364690/150‑1 (MQ=255)
                                                                                                          ggTGCGGGCCGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                             <  2:5183306/150‑1 (MQ=255)
                                                                                                              cGGGCTGACGCGTACAGGAAACCCCGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCACAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTg                                         <  1:1779954/150‑1 (MQ=255)
                                                                                                               gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                        <  1:5825351/150‑1 (MQ=255)
                                                                                                               gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                        <  1:6370172/150‑1 (MQ=255)
                                                                                                                     aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGCGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAATCAGAACGTTTTCTGCGTGTTGCCGatat                                  >  1:3738145/1‑150 (MQ=255)
                                                                                                                      cgcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATAtt                                 <  1:2337356/150‑1 (MQ=255)
                                                                                                                          tACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                             <  2:2284072/150‑1 (MQ=255)
                                                                                                                           aCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                            <  1:6083390/150‑1 (MQ=255)
                                                                                                                           aCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                            >  2:6090071/1‑150 (MQ=255)
                                                                                                                             aGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTTGGACCAAAGGTAAAGGGGTTACAAACCTGCGGGTGTTGAAGGTCGGCGGTACAGCCGTGGCCCATGCAGAACGTTTTGTGCGTGTTGGCGATATTGTGGaa                           >  1:5865614/1‑150 (MQ=255)
                                                                                                                              ggAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                         <  2:1243933/150‑1 (MQ=255)
                                                                                                                                  acacAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAt                     <  1:1422153/150‑1 (MQ=255)
                                                                                                                                   cacaGAACAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATg                    <  2:40028/150‑1 (MQ=255)
                                                                                                                                    acaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                   >  2:3099943/1‑150 (MQ=255)
                                                                                                                                     aaaaaaaaaGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTAGACCAAAGGTAACCAGGTAAAAACCAATCGAGTGTTGAAGTTAGGCGGTAAATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGACGATATTATTGAAAGCAATGcc                  >  1:6508976/4‑150 (MQ=25)
                                                                                                                                      aGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                 <  1:5583098/150‑1 (MQ=255)
                                                                                                                                         aaaaaGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggc              >  2:1430145/1‑150 (MQ=255)
                                                                                                                                          aaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggca             >  2:2858154/1‑150 (MQ=255)
                                                                                                                                            aaGCCCGCTCCTGACAGTGGGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg           <  2:4366796/150‑1 (MQ=255)
                                                                                                                                            aaGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg           >  2:5325906/1‑150 (MQ=255)
                                                                                                                                             aGCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                              <  1:4548903/130‑1 (MQ=255)
                                                                                                                                              gCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGTCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg         <  1:2638948/150‑1 (MQ=255)
                                                                                                                                              gCCCGCACTGTCAGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg         >  2:1597404/1‑150 (MQ=255)
                                                                                                                                               cccGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                             >  2:1758527/1‑113 (MQ=255)
                                                                                                                                                 cGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag      <  1:6090087/150‑1 (MQ=255)
                                                                                                                                                   cACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt    >  2:3909164/1‑150 (MQ=255)
                                                                                                                                                     ccTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg  <  2:4859979/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGG  >  NC_000913/138‑436

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: