Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,307,298 G→A 9.7% P168S (CCC→TCC)  alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,307,2980GA 9.7% 232.2 / 11.6 83P168S (CCC→TCC) alsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (40/34);  new base A (4/4);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

AGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCA  >  NC_000913/4307149‑4307447
                                                                                                                                                     |                                                                                                                                                     
aGGGGCGTTTTCCGCATGTACGAAAAGATCGCGCTATGGGTAATTTCGGGGCTGTTGTTCGTAGCAGCGTCTTCGCGCCATTCCAGAGCATTTTGTTTCCAGGCACCCAGGATTGCCACACGCGCTGTGTTGGGTCATATCTCCCTgggg                                                                                                                                                       <  2:241164/150‑1 (MQ=255)
aGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                       >  2:960140/1‑150 (MQ=255)
 ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                      <  2:256071/150‑1 (MQ=255)
      gTTTTCCGCATTGACCAAAAGATCGCGCAAGTGGTAATTTCGGGTCTGTTGTGCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTCGTGTCCAGGCACCCCGGCTCGCCACACGCGCAGTGTTGGGTCATAGCTCCCAGGGGGATATg                                                                                                                                                 <  2:1247148/150‑1 (MQ=255)
      gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATg                                                                                                                                                 <  2:1285688/150‑1 (MQ=255)
        tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                               <  1:780618/150‑1 (MQ=255)
         ttCCGCATGGACGAAAAGAGCGCGCAATGGGTAATTTCGGGGCGGTTGTTCGGACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGAGTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                                                                                              <  1:581430/150‑1 (MQ=255)
          tCCGCATGGACGAAAAGATCGCGCAATGGGTAATATCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAccc                                                                                                                                             <  2:123447/150‑1 (MQ=255)
          gccGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGGCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAccc                                                                                                                                             <  1:1340339/149‑1 (MQ=255)
           ccGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGACCCa                                                                                                                                            <  2:668542/150‑1 (MQ=255)
           ccGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGACCCa                                                                                                                                            <  2:1290108/150‑1 (MQ=255)
            cGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTCGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGAGATGCCCCAg                                                                                                                                           >  2:66265/1‑150 (MQ=255)
               atgtgcGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTc                                                                                                                                        <  2:1444327/145‑1 (MQ=255)
                 ggACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                                      <  1:363412/150‑1 (MQ=255)
                 ggACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                                      <  2:1447402/150‑1 (MQ=255)
                 ggACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGATCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTTTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGc                                                                                                                                      >  2:1065715/1‑150 (MQ=255)
                  gACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                     >  1:534054/1‑150 (MQ=255)
                  gACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                               >  2:343622/1‑140 (MQ=255)
                  gACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                               <  1:343622/140‑1 (MQ=255)
                    cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                                   >  2:1592220/1‑150 (MQ=255)
                    cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                                   >  2:1057844/1‑150 (MQ=255)
                    cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACACAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                                   >  2:865038/1‑150 (MQ=255)
                    cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGGTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                                   <  1:1358014/150‑1 (MQ=255)
                     gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGc                                                                                                                                  <  2:87443/150‑1 (MQ=255)
                     gAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATATACCAGGGGGATATGAACCAGATAGACTGc                                                                                                                                  >  1:1669488/1‑150 (MQ=255)
                     caaaaGACCGCGCAATTGGTAATTTCGGGGCGGGTCTGCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGc                                                                                                                                  <  2:111742/149‑1 (MQ=255)
                      aaaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGcc                                                                                                                                 >  1:1294959/1‑150 (MQ=255)
                         aGATCGCGAAATGGGTAATTTCGGTGCTGTTGTTCGTACAAGCGTCTTAGCGACATTCCAGAGCTATTGCTTGCCCGGCCCCCAGGATTGCCACACGCGCAGTGTTGGGTCGTATCTCCCAGGGGGATATGACCCCGTTCGCCTGCaacc                                                                                                                              >  2:954231/1‑149 (MQ=255)
                          gATCGCGCAATGGGTAATTTCGGTGCTGTTGTTCGTACCAGCTGCTTATCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCacac                                                                                                                             <  2:1059745/150‑1 (MQ=255)
                              gcgcAATGGGTAATTTCGGGGCTGGTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGCCCCCGTTCGCCTGCCACAccttt                                                                                                                         >  2:76192/1‑147 (MQ=255)
                                 cAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGtgtga                                                                                                                      >  2:960334/1‑149 (MQ=255)
                                      ggTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCg                                                                                                                 >  1:532862/1‑150 (MQ=255)
                                            ttCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACg                                                                                                           >  1:234440/1‑150 (MQ=255)
                                               ggggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGGAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgcg                                                                                                        >  2:113200/1‑150 (MQ=255)
                                                gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGgcgc                                                                                                       <  2:234440/150‑1 (MQ=255)
                                                  gCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGTTATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCCCAc                                                                                                     >  1:808045/1‑150 (MQ=255)
                                                   cTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGGTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAcc                                                                                                    <  1:244140/150‑1 (MQ=255)
                                                    tgttgtTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGACTGCCACACCGTGTGCACCCGTCCACGGCGCACCg                                                                                                   >  2:1073250/1‑150 (MQ=255)
                                                       tgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                                <  1:1399646/150‑1 (MQ=255)
                                                         ttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                                              <  2:1240598/150‑1 (MQ=255)
                                                               ccAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                        >  1:1382880/1‑150 (MQ=255)
                                                               ccAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCACAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCacac                                                                                        >  1:1376225/1‑150 (MQ=255)
                                                                      cTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGCTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGCAc                                                                                 >  1:65813/1‑150 (MQ=255)
                                                                      cTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAc                                                                                 <  1:1592220/150‑1 (MQ=255)
                                                                        tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                                               <  1:1103028/150‑1 (MQ=255)
                                                                          gcgcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCa                                                                             >  1:166185/1‑150 (MQ=255)
                                                                             ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcc                                                                          >  2:313777/1‑150 (MQ=255)
                                                                                   cAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCGCCCCGTTCATCCACACTGCGACCCCCATCCCCGTAc                                                                    >  1:700510/1‑150 (MQ=255)
                                                                                        cAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                               <  1:960140/150‑1 (MQ=255)
                                                                                        cAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGGACCgaga                                                               >  1:292120/1‑150 (MQ=255)
                                                                                              gTTTCCAGGCACCCAGGATTTCCACACGCGCAGTGTTGGGTCATATCTCCCAGTGTGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCTAACCCCATCCCCGTACCGAGATAGgcc                                                         <  1:571175/150‑1 (MQ=255)
                                                                                               tttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccg                                                        <  2:1382880/150‑1 (MQ=255)
                                                                                               tttCCAGGAACCCAGGATTTCCACACGCGCAGTGTTGGGTCATATCTCCAAGTGGGATATGACCAAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccg                                                        <  2:1376225/150‑1 (MQ=255)
                                                                                                   cAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAg                                                    >  2:606762/1‑150 (MQ=255)
                                                                                                     ggCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGaa                                                  >  1:729669/1‑150 (MQ=255)
                                                                                                      gCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAAc                                                 <  1:1561004/150‑1 (MQ=255)
                                                                                                        aCCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCAGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCa                                               >  1:852083/1‑150 (MQ=255)
                                                                                                          ccAGGATTGCCACACGCGCAGTTTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                             <  1:1065715/150‑1 (MQ=255)
                                                                                                           cAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCAACGGCGCACCGTTCATCCACACTGCGACCCCCCTCCCCGTACCGAGATAGGCCGCCAGAACCAgtt                                            >  2:93124/1‑150 (MQ=255)
                                                                                                              gATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                         >  1:79287/1‑150 (MQ=255)
                                                                                                               aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGttgttg                                        <  1:1487252/150‑1 (MQ=255)
                                                                                                                  gCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACCCTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCAGGTGTTGCgg                                     >  2:456377/1‑149 (MQ=255)
                                                                                                                  gCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGt                                     >  2:427481/1‑150 (MQ=255)
                                                                                                                  gCCACACGCGCAGTGTTGGGTAATATCTCCCAGGGGGATATTACCCAGTTAGACTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTTAGAACCACAACCCCGTTGCGAGAAAGGCCTCCAGAAACAGTTGTTGCga                                     >  2:767134/1‑149 (MQ=255)
                                                                                                                   ccACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCCCCAGAACCAGTTGTTGCGTa                                    <  1:825604/150‑1 (MQ=255)
                                                                                                                   ccACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTACGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTa                                    <  1:1125206/150‑1 (MQ=255)
                                                                                                                    cacaCGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTaa                                   >  2:1647173/1‑150 (MQ=255)
                                                                                                                        cgcgCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCCGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGCACCCCCTCCCCGGACCGAGATTGGCCGCCAGCACCAGTTGGTGCGTAAGGCg                               >  1:1347491/1‑150 (MQ=255)
                                                                                                                        cgcgCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCg                               <  1:867863/150‑1 (MQ=255)
                                                                                                                        cgcgCAGTGTTGGGGCATATCTCCCAGGGGTATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCg                               <  1:1063490/150‑1 (MQ=255)
                                                                                                                               tgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCTTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTct                        <  2:1349139/150‑1 (MQ=255)
                                                                                                                                gtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACTTTTTTTTTTAATGTTAAGGGGTCAACGGTGATCTACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTcta                       <  2:3019/150‑1 (MQ=255)
                                                                                                                                 ttGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCAGCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctac                      <  1:93124/150‑1 (MQ=255)
                                                                                                                                  tGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctact                     <  2:166185/150‑1 (MQ=255)
                                                                                                                                   gggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctacta                    >  1:641647/1‑150 (MQ=255)
                                                                                                                                   gggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCACATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTgttcta                    >  1:221336/1‑150 (MQ=255)
                                                                                                                                    ggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCGCTCCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCAGTTGTTGCGTAAGGCGGTTTTctactcc                   >  1:797507/1‑148 (MQ=255)
                                                                                                                                     gTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACGAGTTGTTGCGTAAGGCGGGTTTCTACTACg                  >  1:895298/1‑150 (MQ=255)
                                                                                                                                      tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGt                 >  1:767192/1‑150 (MQ=255)
                                                                                                                                         tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCCCGGCGCACCGTTCATCCACCCTGCGCACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGGTTTCTACTCCGTccc              >  1:704767/1‑150 (MQ=255)
                                                                                                                                         tatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc              <  1:570408/150‑1 (MQ=255)
                                                                                                                                         tatCTCCCAGGGGGATATGACACAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGAGAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTAGTACGTccc              >  2:1404786/1‑150 (MQ=255)
                                                                                                                                            cccccTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCAACCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTATTCTACAACGTCCCAgg           >  2:1478748/3‑150 (MQ=255)
                                                                                                                                             tcCCAGGGGGATTTGACCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCGCACCGTTCATCCACCCTGCGACCCCCATCCCCGTACCGAGATAGGCCGCCGGAACCCGGTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAggg          >  2:1007411/1‑149 (MQ=255)
                                                                                                                                               ccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGGTTTCTACTACGTCCCAGgaga        >  1:364737/1‑150 (MQ=255)
                                                                                                                                                cAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGagag       <  2:532862/150‑1 (MQ=255)
                                                                                                                                                   ggggATATGACCCAGTGAGCCTGCCACACCGTGTGCACCCTTCCACGGCGCACAGGTCCTCCACTCTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCTGAACCTGTTGTGGCGTAAGGCGGTTTTATACTACGTCCCAGGAGAGTTg    <  1:76192/150‑1 (MQ=255)
                                                                                                                                                    gggATATGACCCCGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCCTCCACACTGCGAACCCCCTCCCCGTACCGAGAGAGGCCGCCAGACGCAGTTGTTGCGTAAGGAGTTTTTCTACTACGTCCCAGGAGAGTTGc   >  2:1133249/1‑150 (MQ=255)
                                                                                                                                                    gggATATGACCCAGTTCGCCAGCCACACAGTTTGTACCCGTTCCAGGAGCAGCGTCGATCTACACTGCGAACCCCATCCCCGTACCGAGTTAGGCGGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGc   <  2:1063072/150‑1 (MQ=255)
                                                                                                                                                     ggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTCCCGAGAGAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCa  >  2:104036/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCA  >  NC_000913/4307149‑4307447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑