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.

CAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGG  >  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)
        gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑Gt                                                                                                                                                >  2:4928762/1‑150 (MQ=255)
        gCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑Gt                                                                                                                                                >  2:4797524/1‑150 (MQ=255)
          aTAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGc                                                                                                                                              <  2:5979280/150‑1 (MQ=255)
          aTAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGc                                                                                                                                              <  1:5047316/150‑1 (MQ=255)
           tAGCGCACATACATATAAAAATTACAGAGAACACACCAACCATGAAACGCATTAGCACCACCATAACCACCACCACCACCATTACCACAGGTAACGGTGCGGGCTGACGCTTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCg                                                                                                                                             <  1:2101890/150‑1 (MQ=255)
           tAGCGCACAGACAGATAAAAATTACAGATTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACCGGTAAGGGTGCGGGCTGACGCGTACAGTAAACAAATAAAAAAGCCCGCACCTGACA‑GTGCg                                                                                                                                             <  1:1672712/150‑1 (MQ=255)
            aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTCCCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACTGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                            <  1:6154671/150‑1 (MQ=255)
             gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                           >  1:5742339/1‑150 (MQ=255)
             gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                           <  2:4474894/150‑1 (MQ=255)
             gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                           <  2:6365277/150‑1 (MQ=255)
             gcgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCggg                                                                                                                                           >  2:6154315/1‑150 (MQ=255)
                cacaGAAATATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCACCACCATTCCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAAAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑tt                                                                                                                                       <  1:3177784/150‑1 (MQ=255)
                    gacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACATAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑tttttt                                                                                                                                   <  1:5716947/150‑1 (MQ=255)
                       agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttttt                                                                                                                                >  2:5822756/1‑150 (MQ=255)
                       agaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑ttttttttt                                                                                                                                <  1:728572/150‑1 (MQ=255)
                        gaTAAAAATTACAGCGTCCACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCCCAGGTAACGGTGCGGGCTGACGCGTACCGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTc                                                                                                                               <  1:193706/150‑1 (MQ=255)
                        gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTc                                                                                                                               <  1:1011453/150‑1 (MQ=255)
                         aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATAACAACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACC‑GTGCGGGC‑TTTTTTTTTCg                                                                                                                              >  1:5386544/1‑150 (MQ=255)
                          tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCCCCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGa                                                                                                                             <  1:4381445/150‑1 (MQ=255)
                          tAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGa                                                                                                                             >  2:222286/1‑150 (MQ=255)
                            aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAcc                                                                                                                           >  1:1243933/1‑150 (MQ=255)
                               aTTACAGAGTACACAACATCCCTGAAACGCATTAGCACCCCCCTTACCACCCCCATCACCATTACCACAGCGAACGGGGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGGCCGCACCTGACT‑GTGCGGGC‑TTTTTTTTTCGACCaaa                                                                                                                        <  2:2965988/150‑1 (MQ=255)
                               aTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCaaa                                                                                                                        <  2:111005/150‑1 (MQ=255)
                                 tACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACAACCACCATCACAATTACCACAGGTAACGGTGAGGGCTGACGCGTACAGGAAACACAGAAAAATGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTATAAGGCCAAAgg                                                                                                                      >  1:2308992/1‑150 (MQ=255)
                                                aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                                                                       >  1:4859979/1‑150 (MQ=255)
                                                aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                                                                       >  2:6586278/1‑150 (MQ=255)
                                                aTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTaacaac                                                                                                       >  1:6098751/1‑150 (MQ=255)
                                                  ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                                                                                     >  2:15957/1‑150 (MQ=255)
                                                   cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                    <  1:3841669/150‑1 (MQ=255)
                                                   cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                    <  2:5917407/150‑1 (MQ=255)
                                                   cATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAt                                                                                                    >  1:3669718/1‑150 (MQ=255)
                                                     tGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                  <  1:21985/150‑1 (MQ=255)
                                                       aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCgc                                                                                                >  1:6426319/1‑149 (MQ=255)
                                                           gCATTCGCACCACCATTACCCCCCCCATCACCCTTACCCCAGGTAACGGGGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCGGACA‑GTGCGGGC‑TTTTTTTTTCGACCACAGGTAACGAGGTAACAACCATGCGAgtgt                                                                                            <  2:1845550/150‑1 (MQ=255)
                                                             aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATTGCCCGCAGCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAATGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                          >  1:1487937/1‑150 (MQ=255)
                                                              atAGCACCACCATTACCACCACCACCCCCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACCGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                                                                         <  1:4090107/149‑1 (MQ=255)
                                                               tatcacaacCATTCCCCACACCATCACCATTTCCTCAGGTAACGTTGCGTGCTGTCGCGTACATGAAACACAGAAAAATGCCCGCACCTGACT‑GTGCTGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                        <  1:4433661/147‑1 (MQ=255)
                                                                   accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                    >  1:4248855/1‑150 (MQ=255)
                                                                   accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                    >  1:913033/1‑150 (MQ=255)
                                                                    ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCAGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTACCGAGGTAACAACCATGCGAGTGTTTAAGTTcg                                                                                   >  2:1172099/1‑150 (MQ=255)
                                                                    ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGC‑tttt                                                                                                                                     >  1:9956/1‑100 (MQ=38)
                                                                    ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCAGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGCGTGTTGAAGTTcg                                                                                   >  2:5672625/1‑150 (MQ=255)
                                                                     caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                  >  1:2482163/1‑150 (MQ=255)
                                                                     caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGc                                                                                                                                          <  2:6674037/95‑1 (MQ=38)
                                                                     caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGc                                                                                                                                          >  1:6674037/1‑95 (MQ=38)
                                                                      accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                 >  2:744603/1‑150 (MQ=255)
                                                                             ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACAt                                                                          <  2:2499417/150‑1 (MQ=255)
                                                                              caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                                                                         >  1:5916094/1‑150 (MQ=255)
                                                                               accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                                                                        >  2:2170256/1‑150 (MQ=255)
                                                                               accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGGACATCa                                                                        >  1:2961239/1‑150 (MQ=255)
                                                                                 caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                      <  2:913033/150‑1 (MQ=255)
                                                                                 caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGGACATCAGt                                                                      >  1:1562459/1‑150 (MQ=255)
                                                                                 caccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAGGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                      >  2:2833003/1‑150 (MQ=255)
                                                                                          catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                        >  1:2145349/1‑123 (MQ=255)
                                                                                          catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                             >  2:3354048/1‑150 (MQ=255)
                                                                                           atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                             >  2:5345985/1‑150 (MQ=255)
                                                                                               ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAATTGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAATGGTAACGAGGTAACAACCATGGGAGTGTTGAAGTTCGGCGGTTCATCAGTGGCAAATGCAGaag                                                        >  1:5195068/1‑149 (MQ=255)
                                                                                               ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc                                                        <  2:3438905/150‑1 (MQ=255)
                                                                                                cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGAAGAACg                                                       >  1:4306085/1‑150 (MQ=255)
                                                                                                 acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGAc                                                                                                                            <  2:572700/80‑1 (MQ=255)
                                                                                                     gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                  <  1:188258/150‑1 (MQ=255)
                                                                                                     gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                  <  2:4582815/150‑1 (MQ=255)
                                                                                                     gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCAGACA‑GTGCGGGC‑ATTTTTTTTCGACCAAAGGTAACGAGGTAACCACCAAGCGCGAGTTGAAGTTCGGCGGTACATCAGTGGCAAATTCCGAACGtttta                                                  >  1:3470938/1‑149 (MQ=255)
                                                                                                      tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                 <  2:3669718/150‑1 (MQ=255)
                                                                                                          ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                             <  1:3364690/150‑1 (MQ=255)
                                                                                                          ggTGCGGGCCGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                             <  2:5183306/150‑1 (MQ=255)
                                                                                                              cGGGCTGACGCGTACAGGAAACCCCGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCACAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTg                                         <  1:1779954/150‑1 (MQ=255)
                                                                                                               gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                        <  1:5825351/150‑1 (MQ=255)
                                                                                                               gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                        <  1:6370172/150‑1 (MQ=255)
                                                                                                                     aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGCGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAATCAGAACGTTTTCTGCGTGTTGCCGatat                                  >  1:3738145/1‑150 (MQ=255)
                                                                                                                      cgcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATAtt                                 <  1:2337356/150‑1 (MQ=255)
                                                                                                                          tACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                             <  2:2284072/150‑1 (MQ=255)
                                                                                                                           aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                            <  1:6083390/150‑1 (MQ=255)
                                                                                                                           aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                            >  2:6090071/1‑150 (MQ=255)
                                                                                                                             aGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTTGGACCAAAGGTAAAGGGGTTACAAACCTGCGGGTGTTGAAGGTCGGCGGTACAGCCGTGGCCCATGCAGAACGTTTTGTGCGTGTTGGCGATATTGTGGaa                           >  1:5865614/1‑150 (MQ=255)
                                                                                                                              ggAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                         <  2:1243933/150‑1 (MQ=255)
                                                                                                                                  acacAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAAt                     <  1:1422153/150‑1 (MQ=255)
                                                                                                                                   cacaGAACAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATg                    <  2:40028/150‑1 (MQ=255)
                                                                                                                                    acaGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                   >  2:3099943/1‑150 (MQ=255)
                                                                                                                                     aaaaaaaaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTAGACCAAAGGTAACCAGGTAAAAACCAATCGAGTGTTGAAGTTAGGCGGTAAATCAGTGGCAAATGCAGAACGTTTTATGCGTGTTGACGATATTATTGAAAGCAATGcc                  >  1:6508976/4‑150 (MQ=25)
                                                                                                                                      aGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                 <  1:5583098/150‑1 (MQ=255)
                                                                                                                                         aaaaaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggc              >  2:1430145/1‑150 (MQ=255)
                                                                                                                                          aaaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggca             >  2:2858154/1‑150 (MQ=255)
                                                                                                                                            aaGCCCGCTCCTGACA‑GTGGGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg           <  2:4366796/150‑1 (MQ=255)
                                                                                                                                            aaGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg           >  2:5325906/1‑150 (MQ=255)
                                                                                                                                             aGCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTg                              <  1:4548903/130‑1 (MQ=255)
                                                                                                                                              gCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAACGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGTCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg         <  1:2638948/150‑1 (MQ=255)
                                                                                                                                              gCCCGCA‑CTGTCAGGTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg         >  2:1597404/1‑150 (MQ=255)
                                                                                                                                               cccGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                             >  2:1758527/1‑113 (MQ=255)
                                                                                                                                                 cGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggcag      <  1:6090087/150‑1 (MQ=255)
                                                                                                                                                   cACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGt    >  2:3909164/1‑150 (MQ=255)
                                                                                                                                                     ccTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg  <  2:4859979/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGG  >  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: ‑