Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 144,806 C→A 5.6% S77* (TCA→TAA)  yadI → putative PTS Enzyme IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913144,8060CA 5.6% 495.1 / 12.3 139S77* (TCA→TAA) yadIputative PTS Enzyme IIA
Reads supporting (aligned to +/- strand):  ref base C (66/65);  new base A (4/4);  total (70/69)
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.99e-01

Only 100 of 143 total aligned reads displayed.

TCAGTGCCGGGCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTA  >  NC_000913/144660‑144952
                                                                                                                                                  |                                                                                                                                                  
tcAGTGCCGGGCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTc                                                                                                                                                    >  1:102384/1‑147 (MQ=255)
tcAGTGCCGGGCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTc                                                                                                                                                    <  2:102384/147‑1 (MQ=255)
        gggCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAg                                                                                                                                           <  4:785276/148‑1 (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)
                                                cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACt                                                                                                                                >  4:180293/1‑119 (MQ=255)
                                                cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACt                                                                                                                                <  3:180293/119‑1 (MQ=255)
                                                cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGtt                                                                                                                          <  4:209141/125‑1 (MQ=255)
                                                cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGtt                                                                                                                          >  3:209141/1‑125 (MQ=255)
                                                cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACg                                                                                                   >  4:245572/1‑148 (MQ=255)
                                                cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTAATAATGAAGCca                                                                                                                                        >  3:334308/1‑111 (MQ=255)
                                                cAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTAATAATGAAGCca                                                                                                                                        <  4:334308/111‑1 (MQ=255)
                                                 aGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGt                                                                                                  <  3:245572/148‑1 (MQ=255)
                                                    cGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTcc                                                                                                                              <  4:658209/117‑1 (MQ=255)
                                                    cGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTcc                                                                                                                              >  3:658209/1‑117 (MQ=255)
                                                      cATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCg                                                                                             <  4:312961/148‑1 (MQ=255)
                                                        tgatgTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTaa                                                                                                                                            >  2:141982/1‑99 (MQ=255)
                                                        tgatgTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTaa                                                                                                                                            <  1:141982/99‑1 (MQ=255)
                                                              gCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGa                                                                                     <  4:921693/148‑1 (MQ=255)
                                                                    ctctGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAacac                                                                                                                                 >  3:957961/1‑98 (MQ=255)
                                                                      ctGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGa                                                                                                                      >  3:262406/1‑107 (MQ=255)
                                                                      ctGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGa                                                                                                                      <  4:262406/107‑1 (MQ=255)
                                                                        gCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCa                                                                                                                                                   <  3:1061152/76‑1 (MQ=255)
                                                                        gCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCa                                                                                                                                                   >  4:1061152/1‑76 (MQ=255)
                                                                        gCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCa                                                                                                                                                   >  2:78111/1‑76 (MQ=255)
                                                                           tGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCt                                                                        >  4:1013613/1‑148 (MQ=255)
                                                                                cGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGaa                                                                                                                     <  2:264806/98‑1 (MQ=255)
                                                                                cGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGaa                                                                                                                     >  1:264806/1‑98 (MQ=255)
                                                                                cGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTg                                                                                                           >  2:499931/1‑108 (MQ=255)
                                                                                 ggACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAga                                                                                >  1:500155/1‑134 (MQ=255)
                                                                                      ctGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGa                                                                                     <  1:70373/124‑1 (MQ=255)
                                                                                      ctGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGa                                                                                     >  2:70373/1‑124 (MQ=255)
                                                                                           gAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGa                                                                                                                      <  3:1046053/86‑1 (MQ=255)
                                                                                           gAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGa                                                                                                                      >  4:1046053/1‑86 (MQ=255)
                                                                                                 gtCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTgg                                                                                                          >  2:328660/1‑92 (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)
                                                                                                   cATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAg                                                                                                                    >  4:641143/1‑80 (MQ=255)
                                                                                                   cATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCa                                                                                                     <  2:564139/95‑1 (MQ=255)
                                                                                                   cATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCa                                                                                                     >  1:564139/1‑95 (MQ=255)
                                                                                                   cATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAgatgat                                                                            <  2:234221/120‑1 (MQ=255)
                                                                                                    aTCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCag                                               <  3:1013613/148‑1 (MQ=255)
                                                                                                    aTCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCag                                               <  1:608757/148‑1 (MQ=255)
                                                                                                     tcttctTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAgatgat                                                                            >  4:541946/1‑118 (MQ=255)
                                                                                                     tcttctTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAgatgat                                                                            <  3:541946/118‑1 (MQ=255)
                                                                                                      cttcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGaa                                                                                                                     <  4:373444/76‑1 (MQ=255)
                                                                                                      cttcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGaa                                                                                                                     >  3:373444/1‑76 (MQ=255)
                                                                                                      cttcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATg                                                         <  1:1055018/136‑1 (MQ=255)
                                                                                                        tcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGaaa                                                              >  1:431289/1‑129 (MQ=255)
                                                                                                        tcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGaaa                                                              <  2:431289/129‑1 (MQ=255)
                                                                                                         cttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAAcc                                                            <  1:885482/130‑1 (MQ=255)
                                                                                                           tAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTAATAATGaa                                                                                                                                            <  3:852070/48‑1 (MQ=37)
                                                                                                           tAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTAATAATGaa                                                                                                                                            >  4:852070/1‑48 (MQ=37)
                                                                                                            aaCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTc                                                 <  1:34343/138‑1 (MQ=255)
                                                                                                            aaCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTc                                                 >  2:34343/1‑138 (MQ=255)
                                                                                                             aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGa                                                                                                                 >  1:1033872/1‑73 (MQ=255)
                                                                                                             aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGa                                                                                                                 <  2:1033872/73‑1 (MQ=255)
                                                                                                             aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTa                                                                                          >  3:1097230/1‑96 (MQ=255)
                                                                                                              ccGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAAt                                                                                        >  1:728414/1‑97 (MQ=255)
                                                                                                               cGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTACTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGaa                                                                                    <  2:900383/100‑1 (MQ=255)
                                                                                                                 atatAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTgg                                                                                                          <  1:783301/76‑1 (MQ=255)
                                                                                                                 atatAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTgg                                                                                                          >  2:783301/1‑76 (MQ=255)
                                                                                                                 atatAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACa                                                                                  <  3:584685/100‑1 (MQ=255)
                                                                                                                  tataGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGttt                                          >  3:580629/1‑139 (MQ=255)
                                                                                                                  tataGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGttt                                          <  4:580629/139‑1 (MQ=255)
                                                                                                                   ataGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTg                                                                                                           >  1:205353/1‑73 (MQ=255)
                                                                                                                   ataGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTg                                                                                                           <  2:205353/73‑1 (MQ=255)
                                                                                                                       ccGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAAt                                                                                        >  4:117445/1‑88 (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)
                                                                                                                                 aCCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGt                  >  1:818998/1‑148 (MQ=255)
                                                                                                                                      aTCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAg                                                                                 <  2:358300/80‑1 (MQ=255)
                                                                                                                                      aTCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAg                                                                                 >  1:358300/1‑80 (MQ=255)
                                                                                                                                       tCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGttt                                          >  2:699785/1‑118 (MQ=255)
                                                                                                                                          cgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTg                     >  2:192707/1‑136 (MQ=255)
                                                                                                                                          cgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAg         <  1:260406/148‑1 (MQ=255)
                                                                                                                                          cgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAg         >  4:535975/1‑148 (MQ=255)
                                                                                                                                          cgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTCATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTg                     <  1:192707/136‑1 (MQ=255)
                                                                                                                                            tGGCTTAATAATGAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAAc                                                                                   <  4:691114/72‑1 (MQ=39)
                                                                                                                                            tGGCTTAATAATGAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAAc                                                                                   >  3:691114/1‑72 (MQ=39)
                                                                                                                                            tGGCTTAATAATGAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCgg           >  3:272822/1‑144 (MQ=255)
                                                                                                                                             ggCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTg      <  3:535975/148‑1 (MQ=255)
                                                                                                                                               cTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGag    <  2:818998/148‑1 (MQ=255)
                                                                                                                                                ttCATTATTAAGCCTCAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTAtt                              <  4:935199/121‑1 (MQ=255)
                                                                                                                                                 tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCt                                                                                                            >  4:698093/1‑42 (MQ=255)
                                                                                                                                                 tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCt                                                                                                            <  3:698093/42‑1 (MQ=255)
                                                                                                                                                 tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCa                                                           <  1:843994/91‑1 (MQ=255)
                                                                                                                                                 tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCa                                                           <  3:878500/91‑1 (MQ=255)
                                                                                                                                                 tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCa                                                           >  2:843994/1‑91 (MQ=255)
                                                                                                                                                 tCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGagta  >  2:54554/1‑148 (MQ=255)
                                                                                                                                                  cATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTggg                   >  4:617249/1‑130 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCAGTGCCGGGCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTA  >  NC_000913/144660‑144952

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