Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 257,753 C→T 5.9% I409I (ATC→ATT)  frsA → fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913257,7530CT 5.9% 210.4 / 2.4 68I409I (ATC→ATT) frsAfermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base T (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

CAGGCTACTGGAAGAACGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCG  >  NC_000913/257605‑257899
                                                                                                                                                    |                                                                                                                                                  
cAGGCTACTGGAAGAACGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATc                                                                                                                                                    >  3:206868/1‑149 (MQ=255)
          gaagaaCGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCt                                                                                                                                          >  5:150886/1‑149 (MQ=255)
            agaaCGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTtg                                                                                                                                        <  3:192541/149‑1 (MQ=255)
             gttcGATCCGTTCAGCCCGGAAGAGGACTCACGCTTATTCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCTGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAATCGCTtgt                                                                                                                                       <  8:36104/146‑1 (MQ=255)
               aCGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGaaaaa                                                                                                                                              <  2:34006/140‑1 (MQ=255)
               aCGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTAAGGAAATCACAGACTGGATCGAAAAAcagatagg                                                                                                                                      >  1:34006/1‑141 (MQ=255)
                    ccGTTCAGCCCTGAAGAGGACTGACGCTTAATCACCTCATGATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCTGAATTTTGAGAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTaaa                                                                                                                                 <  5:211532/148‑1 (MQ=255)
                             ccGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTa                                                                                                                        >  8:44070/1‑148 (MQ=255)
                                  agagGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAAtttt                                                                                                                  >  7:117917/1‑149 (MQ=255)
                                  agagGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAAtttt                                                                                                                  <  3:132002/149‑1 (MQ=255)
                                   gatgACTCACGCATAATCACCTCATCATCTGCTGACGGTAAATTAATAGAGATCCCATTTAACCCGGTTTATCGGAATTATCACAAAGGTCTTCAGGAAAGCACCGACTGGATCGAAAAACGCTTCTGTAAAAAATTTGCTAAAtttag                                                                                                                 >  6:93908/4‑147 (MQ=255)
                                   gagGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTg                                                                                                                 <  1:107760/149‑1 (MQ=255)
                                     ggACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGc                                                                                                                >  1:113161/1‑148 (MQ=255)
                                           aCGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGc                                                                                                                                           >  7:43908/1‑115 (MQ=255)
                                           aCGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGc                                                                                                                                           <  8:43908/115‑1 (MQ=255)
                                           aCGCTTAATCAACTCATAATCTGCTGACGGTAAATTATTATAGATACCATTTAACACTGTGTATCGGAATTTTTACAAAGTTCTTCCTGAAATCACCGACTTGATAGAAAAACTATTGTGTTAAAAATTTGCTAAATTTTGCCAAtt                                                                                                           <  3:53429/147‑1 (MQ=255)
                                               ttAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGAAAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTa                                                                                                      >  8:69777/1‑148 (MQ=255)
                                                   tCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACa                                                                                                 <  4:213342/149‑1 (MQ=255)
                                                   tCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACa                                                                                                 >  4:173708/1‑149 (MQ=255)
                                                     aCCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGTTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTTTGTTAAAAATTTGATAAATTTTGCAAATTTGTTAAAACAGt                                                                                               >  1:181227/1‑149 (MQ=255)
                                                      ccTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGtt                                                                                              <  3:172895/149‑1 (MQ=255)
                                                          atcatcTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGc                                                                                            <  3:98286/147‑1 (MQ=255)
                                                          atcatcTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCAt                                                                                          <  4:195349/149‑1 (MQ=255)
                                                             atcTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGTTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATcac                                                                                       <  2:113161/149‑1 (MQ=255)
                                                             atcTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATca                                                                                        >  2:81482/1‑148 (MQ=255)
                                                                  ctgACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATcaca                                                                                      <  4:206868/145‑1 (MQ=255)
                                                                     aCGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGga                                                                                <  5:226634/148‑1 (MQ=255)
                                                                       ggTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGgaga                                                                              >  5:134754/1‑148 (MQ=255)
                                                                          aaaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTTAAAAAGGTCTTCAGTAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACACGAGATAg                                                                           >  8:131629/1‑148 (MQ=255)
                                                                          aaaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATTGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATc                                                                                         <  7:226970/134‑1 (MQ=255)
                                                                          aaaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATTGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATc                                                                                         >  8:226973/1‑134 (MQ=255)
                                                                          aaaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATTGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATc                                                                                         <  7:226973/134‑1 (MQ=255)
                                                                          aaaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATTGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATc                                                                                         >  8:226970/1‑134 (MQ=255)
                                                                           aaTTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCa                                                                         <  7:197156/149‑1 (MQ=255)
                                                                               attaGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCacaaca                                                                                   >  5:130126/1‑135 (MQ=255)
                                                                               attaGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCacaaca                                                                                   <  6:130126/135‑1 (MQ=255)
                                                                                ttaGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGa                                                                     <  3:173708/148‑1 (MQ=255)
                                                                                ttaGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc                                                                    >  2:90669/1‑149 (MQ=255)
                                                                                     gaTCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc                                                               <  3:170404/149‑1 (MQ=255)
                                                                                       tCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc                                                              >  1:151902/1‑148 (MQ=255)
                                                                                         ccATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGCCGTTACCGAg                                                           >  6:98348/1‑149 (MQ=255)
                                                                                           aTTTAACCCTGTGTATCGGAATTTTGACAAAGGGCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGGGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTg                                                         <  5:145789/149‑1 (MQ=255)
                                                                                           aTTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAACTTTTGCCAATATGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTg                                                         <  5:165083/149‑1 (MQ=255)
                                                                                            tttAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTgg                                                        <  8:117917/149‑1 (MQ=255)
                                                                                                 cccGGTGTATCGGATTTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACAcc                                                   <  7:89281/149‑1 (MQ=255)
                                                                                                     gtgtATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGaa                                               <  1:147165/149‑1 (MQ=255)
                                                                                                        tATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGga                                                                                >  1:220840/1‑113 (MQ=255)
                                                                                                        tATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGga                                                                                <  2:220840/113‑1 (MQ=255)
                                                                                                        tATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAaga                                             <  6:30326/148‑1 (MQ=255)
                                                                                                            ggAATTTAGACAAAGGTCTTCAGGAAATCACAGACAGGATCGAAAAACGCATGTGATAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAagag                                            >  4:160396/1‑145 (MQ=255)
                                                                                                                ttttGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGa                                    >  6:34423/1‑149 (MQ=255)
                                                                                                                  ttGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGAt                                   <  1:81482/148‑1 (MQ=255)
                                                                                                                    gACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCa                                 >  7:60945/1‑148 (MQ=255)
                                                                                                                       aaaGGTCTTCAGGAAATCACCGCCTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCaaaa                              >  2:193029/1‑148 (MQ=255)
                                                                                                                       aaaGGTATTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCAAAATTTTGCCAATTTGGTAAAACAGTAGCATCCCAACTGGAGATAGCAAAGAGGTTACCGAGTGGCCCCCCGAAGAGCAGCTTGATCaaaa                              >  4:184214/1‑148 (MQ=255)
                                                                                                                              ttCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTAcc                      >  7:207757/1‑149 (MQ=255)
                                                                                                                              ttCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTAcc                      >  8:26819/1‑149 (MQ=255)
                                                                                                                                cAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTAc                       <  6:150886/146‑1 (MQ=255)
                                                                                                                                 aGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGc                    <  7:44070/148‑1 (MQ=255)
                                                                                                                                   gAAATCACCGCCTGGATCGAAAAAGGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTATAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGCAGAGCAGATTGATCAAAAAATTTACCGCaa                  >  5:238153/1‑147 (MQ=255)
                                                                                                                                   gAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCacaa                                                                                     >  6:19275/1‑81 (MQ=255)
                                                                                                                                          ccGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCg          <  3:13702/149‑1 (MQ=255)
                                                                                                                                           cGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCg          >  2:193872/1‑148 (MQ=255)
                                                                                                                                              cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGTAGAGCAGATTGAt                                   <  5:87116/120‑1 (MQ=255)
                                                                                                                                              cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGtata      <  8:207757/149‑1 (MQ=255)
                                                                                                                                              cTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGtata      >  6:5723/1‑149 (MQ=255)
                                                                                                                                              cTGGATCGAAAAACGCATGTGTTAAAAATTTGCTAAATTTTGCCAATTTGCTAAAACAGTTGCATCACAACAGGAGATAGCAAAGACGATACCGAGTGGACACCCGAAGAGCAGATTGAt                                   >  6:87116/1‑120 (MQ=255)
                                                                                                                                                 gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTc   <  7:131629/149‑1 (MQ=255)
                                                                                                                                                  aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCg  >  6:110031/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAGGCTACTGGAAGAACGATCCGTTCAGCCCGGAAGAGGACTCACGCTTAATCACCTCATCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCG  >  NC_000913/257605‑257899

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑