Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,785,595 G→A A112A (GCG→GCA gpmM → 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,785,5950GA98.2% 185.9 / NA 56A112A (GCG→GCAgpmM2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (29/25);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCG  >  NC_000913/3785353‑3785839
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
aCCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGTAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCCGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATCAAACGAAAAACGCAGGTAAAGCATTACACAt                                                                                                                                                                                                                                              >  1:128666/1‑251 (MQ=255)
      gTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGTAATTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATggg                                                                                                                                                                                                                                        <  1:289656/251‑1 (MQ=255)
       tAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGt                                                                                                                                                                                                                                       >  2:205320/1‑251 (MQ=255)
         aTGGATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTCACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCt                                                                                                                                                                                                                                     >  2:139444/1‑251 (MQ=255)
          tGGATGCACTGTGGGCCACTCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCCTCGTGTATCAGGACCTGACTCGTCTGGACGTTGAACTCAAAGATCGCGCTTTCTTTGCTAATCTAGTGCTGACTGGCGCGGTAGATAAAGCGAACAACGCAGGTAAAGCGGTACACATTATGGGTCTg                                                                                                                                                                                                                                    >  2:58041/1‑251 (MQ=255)
           ggATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTg                                                                                                                                                                                                                                    <  1:132069/250‑1 (MQ=255)
           ggATGCACTGTGGGCCAATCGTCCGCATACCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGc                                                                                                                                                                                                                                   >  2:406285/1‑251 (MQ=255)
                       ggCCAATCGTCCGCATGCCCTGATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAggcg                                                                                                                                                                                                                       <  1:184790/251‑1 (MQ=255)
                                                 aCGCTTCCGGTCTGGAAGTCGGTCTGACTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGACACGAAGATCAcatc                                                                                                                                                                                             >  2:59556/1‑251 (MQ=255)
                                                         ggTCTGGAAGTCGGTCTGCCTGCCCGTCCGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACGCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTGTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGAACACAGCCACGAAGa                                                                                                                                                                                                    <  2:361375/236‑1 (MQ=255)
                                                         ggTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAGGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGa                                                                                                                                                                                                    >  1:361375/1‑236 (MQ=255)
                                                             tGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTACAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGt                                                                                                                                                                                  >  1:21572/1‑250 (MQ=255)
                                                                 aGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAAc                                                                                                                                                                             <  1:306634/251‑1 (MQ=255)
                                                                     ggTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGc                                                                                                                                                                         >  1:278536/1‑251 (MQ=255)
                                                                      gTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGc                                                                                                                                                                         >  1:216167/1‑250 (MQ=255)
                                                                                 cGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAggc                                                                                                                                                                                                                        >  1:138723/1‑192 (MQ=255)
                                                                                 cGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAggc                                                                                                                                                                                                                        <  2:138723/192‑1 (MQ=255)
                                                                                 cGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACgcgg                                                                                                                                                             >  1:122191/1‑251 (MQ=255)
                                                                                  gTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACgcggc                                                                                                                                                            >  1:329258/1‑251 (MQ=255)
                                                                                    cAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGgcgc                                                                                                                                                          >  2:406125/1‑251 (MQ=255)
                                                                                                      gTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTAGACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCa                                                                                                                                        <  2:121398/251‑1 (MQ=255)
                                                                                                              cGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGAGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTCGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGCACGCGGCGCAGAAAAAATCTACCTGCCCGCATTcc                                                                                                                                >  1:262234/1‑251 (MQ=255)
                                                                                                                          ggccggccgCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTcgcg                                                                                                                    <  2:21572/250‑1 (MQ=255)
                                                                                                                            ccggccgCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGac                                                                                                                  >  1:123910/1‑251 (MQ=255)
                                                                                                                                         gtATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCa                                                                                                     <  1:395276/251‑1 (MQ=255)
                                                                                                                                           aTCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGt                                                                                                   >  2:321877/1‑251 (MQ=255)
                                                                                                                                               ggACCTGCCTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCTCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATTGGGCTGCGCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCTATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTTCCCGCACGCATTCCGTGACGGGCTCGACACCCCGCCGCGGAGTGCTg                                                                                               <  2:432369/251‑1 (MQ=255)
                                                                                                                                               ggACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAAGCCAGTGCTGACTGGCGCTGTAGATAAAGCGAAAAACGCAGGTAAATCAGTACACATTATGGTTCTGCTCGCCGCAGGCGGCGTACACAGCCACGAAGAGCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCATAAAAAATCTACCTGCACGCATTCCTTGACGGGCGCGACACCCCGCTGCGCAGTGCTg                                                                                               <  2:432372/251‑1 (MQ=255)
                                                                                                                                                gACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAACAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTg                                                                                               >  1:101811/1‑250 (MQ=255)
                                                                                                                                                    tGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATcc                                                                                          >  1:275578/1‑251 (MQ=255)
                                                                                                                                                      acgcGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTc                                                                                        <  2:156849/248‑1 (MQ=255)
                                                                                                                                                         cGTCTGGACGTTGAAATCAAAGATCGCTCTATCTTTGCTAATCCAGTTCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATTGGTCTGCTCTCCGCAGGCGTCGTACACAGCCACGAAGATCACATCATGTCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGc                                                                                      <  1:47706/250‑1 (MQ=255)
                                                                                                                                                         cGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGc                                                                                      <  1:316068/250‑1 (MQ=255)
                                                                                                                                                          gtgtGGACGTTGAAATAAAAGATCGCTCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATGAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTCCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCAGGGCGATGGTAGAACTTGCAGCTGAACGCGGCTCAGAAAAAATCTACCTGCACTCATTCCTGGACGGTCGCGACACCTCGCCGCGCAGTGCTGAATCCTCGCTg                                                                                    <  2:376139/248‑1 (MQ=255)
                                                                                                                                                            cTGGACGTTGAAATCAAAGCTCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTCCACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTc                                                                                                                                 >  2:32580/1‑204 (MQ=255)
                                                                                                                                                            cTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACAGTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACAGCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCa                                                                                                                                                         <  2:390532/180‑1 (MQ=255)
                                                                                                                                                            cTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACGGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTc                                                                                                                                 <  1:32580/204‑1 (MQ=255)
                                                                                                                                                            cTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGAATGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCa                                                                                                                                                         >  1:390532/1‑180 (MQ=255)
                                                                                                                                                              ggACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGaaa                                                                                 >  1:373766/1‑250 (MQ=255)
                                                                                                                                                                        aTCAAAGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTgaaga                                                                      <  1:406125/251‑1 (MQ=255)
                                                                                                                                                                           aaaGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCCGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGAGACACCCCGCCGCGCAGTGCTGAATACTCGCTGAAAAAATTTgaagaga                                                                    >  1:391970/1‑248 (MQ=255)
                                                                                                                                                                            aaGATCGCGCTTTCTTTGCTAATCCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGagaaa                                                                  <  2:34558/251‑1 (MQ=255)
                                                                                                                                                                                            tGCTAAACCAGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAgggt                                                   <  2:77194/250‑2 (MQ=255)
                                                                                                                                                                                                    cagTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCcg                                           >  1:316434/3‑250 (MQ=255)
                                                                                                                                                                                                    cagTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCcg                                           <  2:316434/248‑1 (MQ=255)
                                                                                                                                                                                                      gTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGAAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGCTGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGACGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAACAATTTGAAGAGAAATTTGCTGCGCAGGGAAAAGGCCGCGt                                        >  2:375823/1‑251 (MQ=255)
                                                                                                                                                                                                          tGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCg                                    >  2:122332/1‑251 (MQ=255)
                                                                                                                                                                                                               ggCGCGGTAGATAAAGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTcc                                 <  2:117355/249‑1 (MQ=255)
                                                                                                                                                                                                                           aaaGCGAAAAACGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGtt                    >  1:408418/1‑250 (MQ=255)
                                                                                                                                                                                                                             aGCGAAAAACGCAGGGACAACAGTACACATTATAGGTCTGCTATCCGCAGGCGACGTACACAGCCCCGAAGATCACATCATGCCGATGGTAGAACTGGCAGATGAACGCGGCGCAGAAAAAATCTACCTGCACGAATTCATTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCAATCATTGGTCGTtact                 >  2:69855/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaCGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGAt          >  2:3917/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aCGCAGGTAAAGCAGAACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATTGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCTCCGCGCAGTTCTGAATCCTCGCTGAAAAAATTTGAAGATAAAGTTGCCGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGGTACTACGCGATg         <  2:262234/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cGCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGATGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATgg        <  1:406285/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gCAGGTAAAGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGACCTGGCAGCTGAACGCGGCGCAGAACAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGAAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATgg        >  1:424434/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ggtaaaGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACcg    <  2:329258/248‑1 (MQ=255)
                                                                                                                                                                                                                                             aaaGCAGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACCCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTTGAAGAGAAATTTGCTGCGCAGGGTAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACcgcg  >  1:349473/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
ACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGCTAATCCGGTGCTGACTGGCGCGGTAGATAAAGCGAAAAACGCAGGCAAAGCGGTACACATTATGGGTCTGCTCTCCGCAGGCGGCGTACACAGCCACGAAGATCACATCATGGCGATGGTAGAACTGGCAGCTGAACGCGGCGCAGAAAAAATCTACCTGCACGCATTCCTTGACGGTCGCGACACTCCGCCGCGCAGTGCTGAATCCTCGCTGAAAAAATTCGAAGAAAAATTTGCCGCGCTGGGCAAAGGCCGCGTAGCGTCCATCATTGGTCGTTACTACGCGATGGACCGCG  >  NC_000913/3785353‑3785839

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