Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,182,939 T→C D238G (GAC→GGC)  yegX ← putative family 25 glycosyl hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,182,9390TC98.5% 215.6 / ‑3.8 65D238G (GAC→GGC) yegXputative family 25 glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (32/32);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTT  >  NC_000913/2182797‑2183083
                                                                                                                                              |                                                                                                                                                
cGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAAtgc                                                                                                                                          >  1:245317/1‑151 (MQ=255)
cGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAAtgc                                                                                                                                          <  1:884267/151‑1 (MQ=255)
 ggggAGTCTGGTGGGAGTATATCACTGCCTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAAtgct                                                                                                                                         >  1:909744/1‑151 (MQ=255)
   ggAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAAtgctgc                                                                                                                                       >  1:192750/1‑151 (MQ=255)
            tggGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAgcgc                                                                                                                              >  1:262696/1‑151 (MQ=255)
             gggAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAgcgc                                                                                                                              <  2:357340/150‑1 (MQ=255)
                             ttcttATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCAtt                                                                                                             >  2:428480/1‑151 (MQ=255)
                               cttATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGt                                                                                                           >  2:126819/1‑151 (MQ=255)
                               cttATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGt                                                                                                           >  2:806984/1‑151 (MQ=255)
                                  aTTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTcc                                                                                                         >  2:60874/1‑150 (MQ=255)
                                  aTTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCg                                                                                                        <  1:244870/151‑1 (MQ=255)
                                   tttAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCTTCCACCGTGCCATTAAATACATTAAAATCTACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCgg                                                                                                       >  2:171492/1‑151 (MQ=255)
                                        ggCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGgacgac                                                                                                  <  2:844781/151‑1 (MQ=255)
                                        ggCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGgacgac                                                                                                  <  1:456378/151‑1 (MQ=255)
                                          cGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGt                                                                                                <  1:720832/151‑1 (MQ=255)
                                            tttCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTg                                                                                              >  1:6614/1‑151 (MQ=255)
                                              tCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAt                                                                                            >  2:469703/1‑151 (MQ=255)
                                              gcgttAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACATGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGa                                                                                             >  2:863166/4‑150 (MQ=255)
                                               cTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAt                                                                                            >  1:307261/1‑150 (MQ=255)
                                                 ttAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAt                                                                                            <  1:292234/148‑1 (MQ=255)
                                                   aaTCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAtagtag                                                                                       <  2:902506/151‑1 (MQ=255)
                                                   aaTCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAtagtag                                                                                       <  1:812840/151‑1 (MQ=255)
                                                   aaTCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGAtagta                                                                                        <  2:414738/150‑1 (MQ=255)
                                                      cccATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATagtagt                                                                                      >  2:437199/1‑149 (MQ=255)
                                                       ccATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAg                                                                                   <  2:285555/151‑1 (MQ=255)
                                                           cAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGccac                                                                               >  1:151843/1‑151 (MQ=255)
                                                               gAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGccacccac                                                                           >  1:507686/1‑151 (MQ=255)
                                                                 aTGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACccacca                                                                         <  2:811575/151‑1 (MQ=255)
                                                                       gCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGata                                                                   <  2:686160/151‑1 (MQ=255)
                                                                       gCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGat                                                                    <  1:103612/150‑1 (MQ=255)
                                                                       gCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGat                                                                    <  2:627420/150‑1 (MQ=255)
                                                                             cctccACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCAtt                                                             >  1:354055/1‑151 (MQ=255)
                                                                             cctccACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCAt                                                              <  1:58007/150‑1 (MQ=255)
                                                                              ctccACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTg                                                            >  2:705102/1‑151 (MQ=255)
                                                                                ccACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGa                                                           <  1:730966/150‑1 (MQ=255)
                                                                                 cACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGaaa                                                         <  1:496681/151‑1 (MQ=255)
                                                                                  aCCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGaaa                                                         >  2:868763/1‑150 (MQ=255)
                                                                                   ccGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAAt                                                        >  1:413006/1‑150 (MQ=255)
                                                                                    cGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAAt                                                        >  2:70357/1‑149 (MQ=255)
                                                                                      tGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGcc                                                    >  2:648084/1‑151 (MQ=255)
                                                                                       gCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGccc                                                   <  2:293158/151‑1 (MQ=255)
                                                                                        ccATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCg                                                  <  1:687013/151‑1 (MQ=255)
                                                                                          aTTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGcc                                                <  1:807049/151‑1 (MQ=255)
                                                                                            tAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGcc                                                >  1:422046/1‑149 (MQ=255)
                                                                                            tAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCa                                               <  1:405720/150‑1 (MQ=255)
                                                                                               aTACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGAtt                                           <  2:682692/151‑1 (MQ=255)
                                                                                                       aaaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGAt                                   >  1:132832/1‑151 (MQ=255)
                                                                                                        aaTCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATa                                  <  2:758864/151‑1 (MQ=255)
                                                                                                          tCCACCGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaa                                 >  2:220775/1‑150 (MQ=255)
                                                                                                            cACCGGACCATTGATGCCATCTACCTGTCCACGGCCGAAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCTCCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaaaaa                              <  1:639338/151‑1 (MQ=255)
                                                                                                             aCCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGCCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATaaaaa                              >  2:281583/1‑150 (MQ=255)
                                                                                                               cGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACg                           >  2:386313/1‑151 (MQ=255)
                                                                                                               cGGACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACg                           <  1:867182/151‑1 (MQ=255)
                                                                                                                ggACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACgg                          >  1:418071/1‑151 (MQ=255)
                                                                                                                ggACCATTGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACgg                          <  2:835626/151‑1 (MQ=255)
                                                                                                                       tGATGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTg                    <  1:617032/150‑1 (MQ=255)
                                                                                                                         aTGCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGt                 >  2:816027/1‑151 (MQ=255)
                                                                                                                           gCCATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTa                >  1:116312/1‑150 (MQ=255)
                                                                                                                            ccATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTaa               >  2:28377/1‑150 (MQ=255)
                                                                                                                             cATCTACCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaat             <  2:418026/151‑1 (MQ=255)
                                                                                                                                  aCCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATc         <  2:446590/150‑1 (MQ=255)
                                                                                                                                  aCCTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCg        >  1:642717/1‑151 (MQ=255)
                                                                                                                                   ccTGTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAaata            <  1:890824/146‑1 (MQ=255)
                                                                                                                                      gTCCACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCt     <  1:576563/150‑1 (MQ=255)
                                                                                                                                        ccACGGCCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCtttt  >  2:591075/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CGGGGAGTCTGGTGGGAGTATATCACTGCTTCTTATTTAAGGCGTTTCTTTAATCCCATCAACGAATGCCTGCAACTCCTCCACCGTGCCATTAAATACATTAAAATCCACCGGACCATTGATGCCATCTACCTGTCCACGGTCGGAATGCTGCCAGAAGCGCCAGGCCATTCCGTCATTGTCCGGACGACGTTGATAGTAGTGAGCCACCCACCACGGATATTCATTGAAATAGCCCGCCAGATTGGTGTGATAAAAAACGGCTCCTGAGTAAATAATCGGCTTTT  >  NC_000913/2182797‑2183083

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