Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,181,988 G→A W298* (TGG→TAG)  yegV → putative kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,181,9880GA94.6% 251.4 / ‑3.2 74W298* (TGG→TAG) yegVputative kinase
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (24/46);  total (27/47)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTGGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCACTTGCTTGATCAC  >  NC_000913/2181840‑2182133
                                                                                                                                                    |                                                                                                                                                 
gacgCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTagg                                                                                                                                                 <  2:841452/151‑1 (MQ=255)
 acgCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTagg                                                                                                                                                 >  2:716240/1‑150 (MQ=255)
      tCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTGAAACCCTTTGGGGCGGCCGACAGCAATGCCGGTGGCGTACTTGCCGGGCTGCCCTCTGGCCTGCCATTGGGGGATGCCGATTTTTTGGGCAATGAAGGGGCGCCGTAGGTTGTCg                                                                                                                                           >  2:939653/1‑151 (MQ=255)
         ggCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGc                                                                                                                                        <  1:1307114/151‑1 (MQ=255)
               aTTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGGTTTACTGGGCAATGCAGTGGCGTCGTGGGTGGGCggggctccgg                                                                                                                                  >  1:1149466/1‑144 (MQ=255)
                   cGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCgg                                                                                                                              <  1:603473/151‑1 (MQ=255)
                       aTTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTg                                                                                                                            >  2:1095201/1‑149 (MQ=255)
                            cAACGCAAGTTTTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTgc                                                                                                                     <  1:1018812/151‑1 (MQ=255)
                             aaCGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGGTGGCGGCTTACTCGGCTAGGCCGTGGCGTCGTGGGTGTTCGGGCGTCGGGGCGGGGAtt                                                                                                                         >  1:91364/1‑146 (MQ=255)
                                 cAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTGGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCaa                                                                                                                <  2:833002/151‑1 (MQ=255)
                                 cAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCaa                                                                                                                <  2:264170/151‑1 (MQ=255)
                                      tGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAAcgcg                                                                                                            <  2:482639/150‑1 (MQ=255)
                                         cGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAg                                                                                                        <  2:1042417/151‑1 (MQ=255)
                                         cGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGGTGTCGGGGAACGGGGCGGGGATTGTGCGCCCAcccgcg                                                                                                          >  1:146271/1‑149 (MQ=255)
                                          gACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGGGATTGTGCGCCAACGCGCGAgg                                                                                                       >  2:1409496/1‑151 (MQ=255)
                                             acCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAac                                                                                                    <  2:1407088/151‑1 (MQ=255)
                                              ccATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAact                                                                                                   <  2:243063/151‑1 (MQ=255)
                                              ccATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGGAATACTGGGGAATGCAGGGGCGTCGTAGGTTGTCGGGCCTCGGGGGGGGGGGTGTGGGCCCAcacccgggggaa                                                                                                    >  1:910489/1‑139 (MQ=255)
                                               cATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAacta                                                                                                  <  2:997972/151‑1 (MQ=255)
                                                aTTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGGGATTGTGCGCCCACGCGCGGGGGactac                                                                                                 >  2:1059454/1‑151 (MQ=255)
                                                aTTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGGGATTGTGCGCCAACGCGCGAGGaa                                                                                                     >  2:172028/1‑147 (MQ=255)
                                                     ggcgggcgACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTActcctc                                                                                            <  1:1167411/151‑1 (MQ=255)
                                                           cgACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGcacac                                                                                      <  1:768815/151‑1 (MQ=255)
                                                           cgACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAAAGCAGGGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGGTTGTGGGCCAACCCGCGAGGGACTACTCCCCGcac                                                                                        >  1:1189367/1‑149 (MQ=255)
                                                             aCAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGAAGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACaa                                                                                    <  2:838544/151‑1 (MQ=255)
                                                              cAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACaaa                                                                                   <  2:132635/151‑1 (MQ=255)
                                                               agccATGCCGGTGGCGTTCTTGCCGGGCTGGCCTCTGGACTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGACGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCGACGCGCGAGGACCTTCTCCTCGCGCACaaa                                                                                   <  1:108082/147‑1 (MQ=255)
                                                               aGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACaaaa                                                                                  <  1:732720/151‑1 (MQ=255)
                                                               aGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACaaaa                                                                                  <  1:989456/151‑1 (MQ=255)
                                                                gTCATGCCTGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGGATTACTGGGCAATGCAGTGGCGTCGTAAGTTGTCGGGGATCGGGGCGGTGATTTTGCGCCAACCCGCGGGGAACCACTTCTCGCACACaaaa                                                                                  >  2:928044/1‑150 (MQ=255)
                                                                  cATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAAc                                                                                <  2:665676/150‑1 (MQ=255)
                                                                    tGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAAc                                                                                <  1:1167247/148‑1 (MQ=255)
                                                                    tGCCGGTGGCGTACTTGCCGGGCGGGCGTCTGGTCTGGCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGt                                                                              <  1:532713/150‑1 (MQ=255)
                                                                     gCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGtat                                                                            <  2:544336/151‑1 (MQ=255)
                                                                        ggTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTTTTAATGGGCAATGAATGGGCGTCGGGGGTTGTCGGGCCTCCGGGCGGGGGTTGTGCGCCCACGCGCGGGGGACTACACCTCGCCCACaaaa                                                                                  >  1:1382045/1‑142 (MQ=255)
                                                                         gTGGCGTACTTGCCGGGCTGGCCTCTGGTCGGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGAt                                                                        <  2:777788/151‑1 (MQ=255)
                                                                           ggCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGCACCCGTAAACCTGGGTGACGCAGTTGCGCCGTAGGTTGTCGGGCATTGGGGCGGCGTGTGTGGGGCAACGCCTGAGGAACAAACCCCCCCCCCCAAAAACGta                                                                             >  1:221976/1‑144 (MQ=255)
                                                                           ggCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATc                                                                       <  1:892565/150‑1 (MQ=255)
                                                                           ggCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCg                                                                      <  2:564793/151‑1 (MQ=255)
                                                                              gTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTg                                                                   <  1:1371419/151‑1 (MQ=255)
                                                                              gTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTg                                                                   <  2:1173679/151‑1 (MQ=255)
                                                                                aCTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCg                                                                 <  1:1311071/151‑1 (MQ=255)
                                                                                  ttGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGGGATTGTGCGCCCACGCGCGAGGAAATACTCCTCGCACACAAAAAAAGATAGATCGCTGc                                                                  >  2:1160862/1‑148 (MQ=255)
                                                                                    gCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACa                                                              <  2:187878/150‑1 (MQ=255)
                                                                                    gCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGGATTACTGGGCAATGCAGTGGCGTCGTAGGGTGTCGGGCATCGGGGCGGGGGTTGTGCGCCCACCCGCGGGGGAATACTCCTCGCCCACAAAAAAGTATAGATCGCTGCg                                                                 >  1:1064636/1‑147 (MQ=255)
                                                                                     ccGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGt                                                            <  2:645599/151‑1 (MQ=255)
                                                                                     ccGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGt                                                            <  1:1404525/151‑1 (MQ=255)
                                                                                       gggCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTgg                                                          <  2:589220/151‑1 (MQ=255)
                                                                                          cTGGCCTCTGGTCTGCCACTGGCGGATGTCGTCTTTATCGGCAATGTAGTGGCGCCCGAGGTTGTCGGGCATCGGGGTGGTGAATGCTGGCCCACACGCGAGGGACTACTCCTCCCACCCAAAAAAGTATAAAGCCCCGCGACAGGGGCtc                                                       >  1:379586/1‑150 (MQ=255)
                                                                                             gCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATg                                                    <  2:911801/151‑1 (MQ=255)
                                                                                                 ctGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCt                                                  <  2:137088/149‑1 (MQ=255)
                                                                                                  tGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGGGGCGTAGGCTGTCGGGCATCGGGGCGGGGGTTGTGCGCCAACCCGCGAGGAACTAATCCTCGCCCACAAAAACGGATAGATCGCTGCGGCAGTGGGTAAAGCtt                                                 >  1:180088/1‑148 (MQ=255)
                                                                                                  tGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGGTTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGGTGCGGCAgggggcaaagct                                                  >  1:527367/1‑140 (MQ=255)
                                                                                                        gctaCTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCa                                           <  2:623131/146‑1 (MQ=255)
                                                                                                           aCTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAgg                                      <  2:1412136/151‑1 (MQ=255)
                                                                                                           aCTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAgg                                      <  2:416184/151‑1 (MQ=255)
                                                                                                              ggCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGc                                     <  1:538296/149‑1 (MQ=255)
                                                                                                               gCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGGGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGt                                  <  1:100913/151‑1 (MQ=255)
                                                                                                                cGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGtt                                 <  1:1294750/151‑1 (MQ=255)
                                                                                                                    tGCCGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTg                             >  2:463684/1‑151 (MQ=255)
                                                                                                                      ccGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGtt                           <  2:117968/151‑1 (MQ=255)
                                                                                                                       cGTATTACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGtt                           >  1:470948/1‑150 (MQ=255)
                                                                                                                           ttACTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTACTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCCAATCCTGTACTCAATAGGCCGTTGCTGTTGGTc                       >  1:1027819/1‑150 (MQ=255)
                                                                                                                            tACTGGGCAATGCTGTGGCGTCGTATGTTGGCGGGCATCGGGTCTGTGAGTGTGCGCCAACGCGCGATGAACTACTCCTCGCGCACCAAAACGTATAGATCGCGGCGACACTGGCTAATGCTGTACTCAATAGGGCGTTGCTGTTg                          <  1:1081999/146‑1 (MQ=255)
                                                                                                                             aCTGGGCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCaa                     <  2:1282305/150‑1 (MQ=255)
                                                                                                                                gggCAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAgcg                  <  1:1296114/150‑1 (MQ=255)
                                                                                                                                   cAATGCAGTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCAc              >  2:748243/1‑151 (MQ=255)
                                                                                                                                          gTGGCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCActtgcttg       <  1:1148011/151‑1 (MQ=255)
                                                                                                                                            ggCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGGCGTTGCTGTTGGTCAAGCGCCACTTGCTTGa      >  1:790418/1‑150 (MQ=255)
                                                                                                                                             gCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCACTTGCTTGAt     >  2:1217476/1‑150 (MQ=255)
                                                                                                                                             gCGTCGTAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCACTTGCTTGATc    <  2:241294/151‑1 (MQ=255)
                                                                                                                                              cgtcgtAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTTGTCAAGCGCCACTTGCTTGATc    >  1:1031037/1‑150 (MQ=255)
                                                                                                                                              cgtcgtAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCACTTGCTTGATCa   <  1:323054/151‑1 (MQ=255)
                                                                                                                                              cgtcgtAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACCCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCCGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCActtgcttg       >  2:82400/1‑147 (MQ=255)
                                                                                                                                              cgtcgtAGGTTGTCGGGCACCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCACTTGCTTGATCa   <  1:1387019/151‑1 (MQ=255)
                                                                                                                                               gtcgtAGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCACTTGCTTGATCAc  >  2:504142/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GACGCTTCTGGCTGCATTCCGGCATTTCCAACGCAAGTTGTCGACACCATTGGGGCGGGCGACAGTCATGCCGGTGGCGTACTTGCCGGGCTGGCCTCTGGTCTGCCACTGGCGGATGCCGTATTACTGGGCAATGCAGTGGCGTCGTGGGTTGTCGGGCATCGGGGCGGTGATTGTGCGCCAACGCGCGAGGAACTACTCCTCGCACACAAAAACGTATAGATCGCTGCGACAGTGGCTAATGCTGTACTCAATAGGCCGTTGCTGTTGGTCAAGCGCCACTTGCTTGATCAC  >  NC_000913/2181840‑2182133

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