Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 144,812 T→G 5.8% L79* (TTA→TGA)  yadI → putative PTS Enzyme IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913144,8120TG 5.8% 510.3 / 17.0 137L79* (TTA→TGA) yadIputative PTS Enzyme IIA
Reads supporting (aligned to +/- strand):  ref base T (65/64);  new base G (4/4);  total (69/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 153 total aligned reads displayed.

GGGCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTT  >  NC_000913/144668‑144958
                                                                                                                                                |                                                                                                                                                  
    cGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCcac                                                                                                                                             >  3:312961/1‑148 (MQ=255)
                            ctctAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGAtt                                                                                                                     >  1:340115/1‑148 (MQ=255)
                             tctAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGAttt                                                                                                                    >  3:556404/1‑148 (MQ=255)
                                aaTATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTg                                                                                                                 >  3:921693/1‑148 (MQ=255)
                                 atatGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTgg                                                                                                                <  4:556404/148‑1 (MQ=255)
                                       tCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCAc                                                                                                          >  3:258335/1‑148 (MQ=255)
                                        cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGtt                                                                                                                                >  3:209141/1‑125 (MQ=255)
                                        cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGtt                                                                                                                                <  4:209141/125‑1 (MQ=255)
                                        cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTAATAATGAAGCca                                                                                                                                              >  3:334308/1‑111 (MQ=255)
                                            cGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTcc                                                                                                                                    <  4:658209/117‑1 (MQ=255)
                                              cATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCg                                                                                                   <  4:312961/148‑1 (MQ=255)
                                                tgatgTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTaa                                                                                                                                                  <  1:141982/99‑1 (MQ=255)
                                                tgatgTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTaa                                                                                                                                                  >  2:141982/1‑99 (MQ=255)
                                                      gCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTcc                                                                                                                                    <  2:834601/107‑1 (MQ=255)
                                                      gCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTcc                                                                                                                                    >  1:834601/1‑107 (MQ=255)
                                                            ctctGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAacac                                                                                                                                       >  3:957961/1‑98 (MQ=255)
                                                              ctGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGa                                                                                                                            >  3:262406/1‑107 (MQ=255)
                                                              ctGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGa                                                                                                                            <  4:262406/107‑1 (MQ=255)
                                                                 cATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTc                                                                                                                                     <  3:676659/95‑1 (MQ=255)
                                                                 cATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTc                                                                                                                                     >  4:676659/1‑95 (MQ=255)
                                                                 cATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTa                                                                                                <  3:778894/132‑1 (MQ=255)
                                                                 cATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTAACGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTa                                                                                                >  4:778894/1‑132 (MQ=255)
                                                                   tGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCt                                                                              >  4:1013613/1‑148 (MQ=255)
                                                                       gCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGcc                                                                          <  2:340115/148‑1 (MQ=255)
                                                                        cGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGaa                                                                                                                           >  1:264806/1‑98 (MQ=255)
                                                                        cGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGaa                                                                                                                           <  2:264806/98‑1 (MQ=255)
                                                                        cGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTg                                                                                                                 <  1:499931/108‑1 (MQ=255)
                                                                        cGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTg                                                                                                                 >  2:499931/1‑108 (MQ=255)
                                                                              ctGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGa                                                                                           <  1:70373/124‑1 (MQ=255)
                                                                              ctGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGa                                                                                           >  2:70373/1‑124 (MQ=255)
                                                                                   gAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGa                                                                                                                            >  4:1046053/1‑86 (MQ=255)
                                                                                         gtCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTgg                                                                                                                <  1:328660/92‑1 (MQ=255)
                                                                                         gtCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGtt                                                        >  1:1126565/1‑148 (MQ=255)
                                                                                          tCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTc                                                       <  4:258335/148‑1 (MQ=255)
                                                                                           cATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGa                                                                                                                            >  4:268130/1‑78 (MQ=255)
                                                                                           cATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGa                                                                                                                            <  3:268130/78‑1 (MQ=255)
                                                                                           cATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAg                                                                                                                          <  3:641143/80‑1 (MQ=255)
                                                                                           cATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCa                                                                                                           <  2:564139/95‑1 (MQ=255)
                                                                                           cATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCa                                                                                                           >  1:564139/1‑95 (MQ=255)
                                                                                            aTCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCag                                                     <  1:608757/148‑1 (MQ=255)
                                                                                            aTCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCag                                                     <  3:1013613/148‑1 (MQ=255)
                                                                                             tcttctTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAgatgat                                                                                  <  3:541946/118‑1 (MQ=255)
                                                                                             tcttctTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCtt                                                                             >  3:415851/1‑123 (MQ=255)
                                                                                              cttcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGaa                                                                                                                           <  4:373444/76‑1 (MQ=255)
                                                                                              cttcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATg                                                               >  2:1055018/1‑136 (MQ=255)
                                                                                              cttcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATg                                                               <  1:1055018/136‑1 (MQ=255)
                                                                                                 cttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAAcc                                                                  >  2:885482/1‑130 (MQ=255)
                                                                                                 cttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAAcc                                                                  <  1:885482/130‑1 (MQ=255)
                                                                                                   tAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTcg                                              >  4:547117/1‑148 (MQ=255)
                                                                                                   tAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTAATAATGaa                                                                                                                                                  >  4:852070/1‑48 (MQ=37)
                                                                                                    aaCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTc                                                       >  2:34343/1‑138 (MQ=255)
                                                                                                     aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGa                                                                                                                       >  1:1033872/1‑73 (MQ=255)
                                                                                                     aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGa                                                                                                                       <  2:1033872/73‑1 (MQ=255)
                                                                                                     aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTa                                                                                                <  4:1097230/96‑1 (MQ=255)
                                                                                                     aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTa                                                                                                >  3:1097230/1‑96 (MQ=255)
                                                                                                      ccGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGa                                                                                                                       >  4:238984/1‑72 (MQ=255)
                                                                                                      ccGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGa                                                                                                                       <  3:238984/72‑1 (MQ=255)
                                                                                                       cGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGaa                                                                                          >  1:900383/1‑100 (MQ=255)
                                                                                                       cGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTACTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGaa                                                                                          <  2:900383/100‑1 (MQ=255)
                                                                                                         atatAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTgg                                                                                                                >  2:783301/1‑76 (MQ=255)
                                                                                                         atatAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTgg                                                                                                                <  1:783301/76‑1 (MQ=255)
                                                                                                         atatAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACa                                                                                        >  4:584685/1‑100 (MQ=255)
                                                                                                          tataGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGttt                                                >  3:580629/1‑139 (MQ=255)
                                                                                                           ataGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTg                                                                                                                 <  2:205353/73‑1 (MQ=255)
                                                                                                           ataGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTg                                                                                                                 >  1:205353/1‑73 (MQ=255)
                                                                                                               ccGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAAt                                                                                              >  4:117445/1‑88 (MQ=255)
                                                                                                                cGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTc                                 >  1:383601/1‑148 (MQ=255)
                                                                                                                  gAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCa                                                                 >  1:842849/1‑114 (MQ=255)
                                                                                                                  gAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCa                                                                 <  2:842849/114‑1 (MQ=255)
                                                                                                                      gcCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTg                           <  2:1126565/148‑1 (MQ=255)
                                                                                                                         aCCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGaa                                                                                          >  1:436087/1‑82 (MQ=255)
                                                                                                                         aCCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGaa                                                                                          <  2:436087/82‑1 (MQ=255)
                                                                                                                              aTCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAg                                                                                       >  1:358300/1‑80 (MQ=255)
                                                                                                                              aTCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGcc                   <  2:383601/148‑1 (MQ=255)
                                                                                                                               tCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGttt                                                >  2:699785/1‑118 (MQ=255)
                                                                                                                                  cgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTg                           >  2:192707/1‑136 (MQ=255)
                                                                                                                                  cgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAg               >  4:535975/1‑148 (MQ=255)
                                                                                                                                  cgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAg               <  1:260406/148‑1 (MQ=255)
                                                                                                                                  cgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAg               >  2:260406/1‑148 (MQ=255)
                                                                                                                                  cgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTCATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTg                           <  1:192707/136‑1 (MQ=255)
                                                                                                                                    tGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATg                                                               <  3:259345/98‑1 (MQ=255)
                                                                                                                                    tGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATg                                                               >  4:259345/1‑98 (MQ=255)
                                                                                                                                     ggCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTg            <  3:535975/148‑1 (MQ=255)
                                                                                                                                        ttCATTATTAAGCCTCAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTAtt                                    <  4:935199/121‑1 (MQ=255)
                                                                                                                                         tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCt                                                                                                                  <  3:698093/42‑1 (MQ=255)
                                                                                                                                         tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCt                                                                                                                  >  4:698093/1‑42 (MQ=255)
                                                                                                                                         tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCa                                                                 <  3:878500/91‑1 (MQ=255)
                                                                                                                                         tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCa                                                                 >  4:878500/1‑91 (MQ=255)
                                                                                                                                         tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTc                                               >  3:1033627/1‑109 (MQ=255)
                                                                                                                                         tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGagta        >  2:54554/1‑148 (MQ=255)
                                                                                                                                          cATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTggg                         <  3:617249/130‑1 (MQ=255)
                                                                                                                                          cATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTggg                         >  4:617249/1‑130 (MQ=255)
                                                                                                                                             tattaAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAgg              >  4:118237/1‑138 (MQ=255)
                                                                                                                                             tattaAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAgg              <  3:118237/138‑1 (MQ=255)
                                                                                                                                             tattaAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCt    >  4:412665/1‑148 (MQ=255)
                                                                                                                                              attaAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCtt   <  3:412665/148‑1 (MQ=255)
                                                                                                                                               ttaAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAAt                                                                                              <  3:1140510/56‑1 (MQ=255)
                                                                                                                                               ttaAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTcgcg                                            <  4:528628/106‑1 (MQ=255)
                                                                                                                                                taAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTg                           >  2:729794/1‑122 (MQ=255)
                                                                                                                                                taAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTg                           <  1:729794/122‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GGGCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTTT  >  NC_000913/144668‑144958

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