Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,797 T→G 6.8% N581H (AAC→CAC)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,7970TG 6.8% 1270.4 / 73.1 369N581H (AAC→CAC) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (169/174);  new base G (11/14);  total (181/188)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

Only 100 of 379 total aligned reads displayed.

AGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAA  >  NC_000913/3832649‑3832945
                                                                                                                                                    |                                                                                                                                                    
aGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGtt                                                                                                                                                      <  3:430434/149‑1 (MQ=255)
 ggTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGACTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTc                                                                                                                                                     <  1:18974/149‑1 (MQ=255)
      aaCTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcg                                                                                                                                                >  1:294633/1‑149 (MQ=255)
          gCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCa                                                                                                                                            >  4:213487/1‑149 (MQ=255)
          gCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCa                                                                                                                                            >  1:452758/1‑149 (MQ=255)
          gCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcg                                                                                                                                            <  2:317509/149‑2 (MQ=255)
           cacaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAg                                                                                                                                           <  4:371206/149‑1 (MQ=255)
           cacaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCACCATCGGCGTACCCCGCGCGTGCGAACgcgcgg                                                                                                                                           <  4:441892/149‑3 (MQ=255)
            acaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcggg                                                                                                                                          <  4:480071/149‑4 (MQ=255)
              aTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCtt                                                                                                                                        >  3:296495/1‑149 (MQ=255)
               tCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTc                                                                                                                                       >  1:386913/1‑149 (MQ=255)
               tCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTc                                                                                                                                       >  1:383095/1‑149 (MQ=255)
                cGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCa                                                                                                                                      <  3:560999/149‑1 (MQ=255)
                  cccGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCa                                                                                                                                      <  1:100497/147‑1 (MQ=255)
                     gCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATa                                                                                                                                 <  2:58655/149‑1 (MQ=255)
                       ttCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACa                                                                                                                               <  2:99510/149‑1 (MQ=255)
                         cAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGa                                                                                                                             <  2:486138/149‑1 (MQ=255)
                           gTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATa                                                                                                                                 >  2:269329/1‑143 (MQ=255)
                            tGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATAc                                                                                                                          <  4:495590/149‑1 (MQ=255)
                            tGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATAc                                                                                                                          <  4:508397/149‑1 (MQ=255)
                            tGAACACCGGCGCAACCATCACGTTATCACCTAACATGTATTGACGGTGAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATAc                                                                                                                          <  4:53023/149‑1 (MQ=255)
                             gAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACg                                                                                                                         <  2:544627/149‑1 (MQ=255)
                                cacCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGTTCTTCCGGGAACTCCATCATCTTGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCGCGCGCAGCTTCACGATACAGATACGGca                                                                                                                      <  2:427051/149‑1 (MQ=255)
                                    ggCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcatcat                                                                                                                  <  2:377940/149‑1 (MQ=255)
                                      cgcAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGc                                                                                                                <  4:328555/149‑1 (MQ=255)
                                       gcAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCg                                                                                                               >  3:143372/1‑149 (MQ=255)
                                         aaCCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGc                                                                                                             >  3:50464/1‑149 (MQ=255)
                                           ccATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCAt                                                                                                           <  3:526259/149‑1 (MQ=255)
                                            cATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCAtt                                                                                                          <  2:246434/149‑1 (MQ=255)
                                             aTCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATTCGGCAttt                                                                                                         <  3:412035/149‑1 (MQ=255)
                                               cACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCa                                                                                                       <  3:266176/149‑1 (MQ=255)
                                                aCGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAg                                                                                                      <  4:384020/149‑1 (MQ=255)
                                                aCGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGCCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAg                                                                                                      <  2:110328/149‑1 (MQ=255)
                                                  gTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                                                                    <  3:312590/149‑1 (MQ=255)
                                                  gTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                                                                    <  1:280307/149‑1 (MQ=255)
                                                    tGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGc                                                                                                  <  4:246646/149‑1 (MQ=255)
                                                      tCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGTATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGt                                                                                                <  3:235016/149‑1 (MQ=255)
                                                        gCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgg                                                                                                                                           >  4:558726/1‑102 (MQ=255)
                                                          cTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaag                                                                                            <  4:541274/149‑1 (MQ=255)
                                                           tAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaaga                                                                                           >  2:81624/1‑149 (MQ=255)
                                                            aaCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                                                                          <  2:440245/149‑1 (MQ=255)
                                                            aaCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                                                                          <  4:157758/149‑1 (MQ=255)
                                                             aCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                                                                        <  4:310267/134‑1 (MQ=255)
                                                             aCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGaagaag                                                                                         <  1:551737/149‑1 (MQ=255)
                                                             aCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGaagaag                                                                                         >  4:341363/1‑149 (MQ=255)
                                                                 gTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                                                                     <  1:583935/149‑1 (MQ=255)
                                                                     tGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgc                                                                                                                                               >  1:359593/1‑87 (MQ=255)
                                                                      gACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACcaca                                                                                >  2:525856/1‑149 (MQ=255)
                                                                          gTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                                                                    <  1:510608/125‑1 (MQ=255)
                                                                          gTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                                                                    >  2:510608/1‑125 (MQ=255)
                                                                          gTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                                                                            <  4:60128/149‑1 (MQ=255)
                                                                            cAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAg                                                                          <  3:91698/149‑1 (MQ=255)
                                                                              aGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGa                                                                        <  1:132070/149‑1 (MQ=255)
                                                                                   aTCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcat                                                                   <  1:107522/149‑1 (MQ=255)
                                                                                     cacaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                                                                          <  1:489736/124‑1 (MQ=255)
                                                                                      acaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCAT                                                                <  1:259014‑M1/149‑1 (MQ=255)
                                                                                      acaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCAT                                                                <  4:146919‑M1/149‑1 (MQ=255)
                                                                                        aaGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCG                                                              >  3:440046‑M1/1‑149 (MQ=255)
                                                                                          gCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa                                                            <  1:94555/149‑1 (MQ=255)
                                                                                            cGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAgg                                                          <  4:459023/149‑1 (MQ=255)
                                                                                                   tCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACgg                                                   >  1:99796/1‑149 (MQ=255)
                                                                                                       ggAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGAGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACa                                               >  2:279980‑M2/1‑148 (MQ=255)
                                                                                                          aCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACacga                                            <  1:182157‑M2/149‑5 (MQ=255)
                                                                                                            tCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATa                                          <  1:2358/149‑1 (MQ=255)
                                                                                                            tCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATa                                          >  4:114789/1‑149 (MQ=255)
                                                                                                            tCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATa                                          >  3:282937/1‑149 (MQ=255)
                                                                                                               atcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGa                                       <  1:61659/149‑1 (MQ=255)
                                                                                                               atcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGa                                       <  2:92451/149‑1 (MQ=255)
                                                                                                                tcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAt                                      >  3:461226/1‑149 (MQ=255)
                                                                                                                 catcatGCCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAtt                                     <  1:330544/149‑1 (MQ=255)
                                                                                                                   tcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTg                                   >  2:175485/1‑149 (MQ=255)
                                                                                                                   tcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTg                                   >  2:169971/1‑149 (MQ=255)
                                                                                                                     atGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                                                                     >  1:461321/1‑97 (MQ=255)
                                                                                                                     atGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                                                                     <  2:461321/97‑1 (MQ=255)
                                                                                                                      tGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTa                                >  3:556383/1‑149 (MQ=255)
                                                                                                                        gCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTAcc                              <  1:90399/149‑1 (MQ=255)
                                                                                                                          ccGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgt                            <  3:20143/149‑1 (MQ=255)
                                                                                                                           cGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgtg                           <  4:519693/149‑1 (MQ=255)
                                                                                                                            gCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgtgt                          >  4:357237/1‑149 (MQ=255)
                                                                                                                                 atcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAgc                                                                                        <  3:543161/82‑1 (MQ=255)
                                                                                                                                    ggCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCt                  <  3:37554/149‑1 (MQ=255)
                                                                                                                                    ggCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCt                  >  3:339499/1‑149 (MQ=255)
                                                                                                                                     gCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTa                 >  3:188075/1‑149 (MQ=255)
                                                                                                                                     gCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTa                 >  1:210795/1‑149 (MQ=255)
                                                                                                                                       gTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATg               >  1:15679/1‑149 (MQ=255)
                                                                                                                                           gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTg           >  3:573029/4‑149 (MQ=255)
                                                                                                                                           gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTg           >  2:188000/4‑149 (MQ=255)
                                                                                                                                           cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTg           <  1:420976/149‑1 (MQ=255)
                                                                                                                                            ctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTgg          >  3:120782/3‑149 (MQ=255)
                                                                                                                                            ccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTgg          >  4:418588/1‑149 (MQ=255)
                                                                                                                                            ccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTgg          >  3:168028/1‑149 (MQ=255)
                                                                                                                                             tgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACt                                                                      <  2:215290/87‑1 (MQ=255)
                                                                                                                                             cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgt                            <  4:300044/130‑1 (MQ=255)
                                                                                                                                             cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGga         <  4:285755/149‑1 (MQ=255)
                                                                                                                                             cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGga         >  2:129069/1‑149 (MQ=255)
                                                                                                                                             cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGga         >  2:405913/1‑149 (MQ=255)
                                                                                                                                              gcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGgag        <  3:481697/149‑1 (MQ=255)
                                                                                                                                                gcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGagag      <  4:271083/149‑1 (MQ=255)
                                                                                                                                                  gTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGca    >  4:131781/1‑149 (MQ=255)
                                                                                                                                                  gTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGca    >  3:241830/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGTAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAA  >  NC_000913/3832649‑3832945

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: