Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,587,425 T→C T104A (ACC→GCC)  ydeQ ← putative fimbrial adhesin protein YdeQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,587,4250TC98.5% 228.4 / ‑5.9 66T104A (ACC→GCC) ydeQputative fimbrial adhesin protein YdeQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (30/35);  minor base A (0/1);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCT  >  NC_000913/1587291‑1587550
                                                                                                                                      |                                                                                                                             
gACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACat                                                                                                                          <  2:115506/140‑1 (MQ=255)
gACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACat                                                                                                                          <  1:525977/140‑1 (MQ=255)
  cTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACatta                                                                                                                        <  1:679307/140‑1 (MQ=255)
   ttCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACattat                                                                                                                       <  2:406040/140‑1 (MQ=255)
               tATGGCAACTCCACCGGCCGCCCCAACTGGTGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAAAATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGACGCATTATTCCAGTACAGCg                                                                                                           <  2:31729/140‑1 (MQ=255)
               tATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCg                                                                                                           <  2:102608/140‑1 (MQ=255)
                aTGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGa                                                                                                          <  2:644188/140‑1 (MQ=255)
                  gACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATc                                                                                                        <  1:388766/140‑1 (MQ=255)
                      aCTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACattatt                                                                                                                      >  2:471332/1‑122 (MQ=255)
                      aCTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACattatt                                                                                                                      <  1:471332/122‑1 (MQ=255)
                      aCTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATc                                                                                                        <  1:503813/136‑1 (MQ=255)
                      aCTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATc                                                                                                        >  2:503813/1‑136 (MQ=255)
                         ccACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTtata                                                                                                 <  2:727879/140‑1 (MQ=255)
                          cACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTatat                                                                                                >  2:467889/1‑140 (MQ=255)
                            ccagccgccCCAACTGGGGTAATATAGAGTCTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGa                                                                                              <  2:492696/137‑1 (MQ=255)
                             cGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGAc                                                                                             <  2:286890/140‑1 (MQ=255)
                             cGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATCATTCCAGTACAGCGATCCTTTATATGAc                                                                                             <  1:544653/140‑1 (MQ=255)
                                    cccAACTGGGGTAATATAGAGTTTTAACGGTAATGACATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAggg                                                                                      >  2:21681/1‑140 (MQ=255)
                                              gTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGaa                                                                            >  1:403387/1‑140 (MQ=255)
                                              gTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGaaa                                                                           >  2:127537/1‑140 (MQ=255)
                                                      gagTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGc                                                                    >  1:332916/1‑140 (MQ=255)
                                                      gagTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGc                                                                    >  2:688020/1‑140 (MQ=255)
                                                       agTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGcc                                                                   >  2:213004/1‑140 (MQ=255)
                                                       agTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGcc                                                                   >  2:778607/1‑140 (MQ=255)
                                                         ttttAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCtt                                                                 >  1:545623/1‑140 (MQ=255)
                                                            tAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTa                                                              >  1:131982/1‑140 (MQ=255)
                                                            tAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTa                                                              >  2:698962/1‑140 (MQ=255)
                                                            tAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTa                                                              <  1:213004/140‑1 (MQ=255)
                                                            tAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTa                                                              >  2:366924/1‑140 (MQ=255)
                                                            tAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTa                                                              <  2:219730/140‑1 (MQ=255)
                                                               cGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCa                                                           <  1:688020/140‑1 (MQ=255)
                                                                 gTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTAcc                                                            >  2:372537/1‑137 (MQ=255)
                                                                 gTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTAcc                                                            <  1:372537/137‑1 (MQ=255)
                                                                 gTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAgg                                                         <  2:403387/140‑1 (MQ=255)
                                                                  tAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTAc                                                             <  1:488831/135‑1 (MQ=255)
                                                                  tAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTAc                                                             >  2:488831/1‑135 (MQ=255)
                                                                   aaTGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGtt                                                       >  2:172110/1‑140 (MQ=255)
                                                                          tCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGAt                                               <  1:127537/140‑1 (MQ=255)
                                                                           cGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGAt                                               <  2:539479/140‑1 (MQ=255)
                                                                           cGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGAt                                               <  2:74663/140‑1 (MQ=255)
                                                                           cGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGAt                                               <  1:661294/140‑1 (MQ=255)
                                                                           cGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGAt                                               <  1:172110/140‑1 (MQ=255)
                                                                            ggTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATc                                              <  1:21681/140‑1 (MQ=255)
                                                                             gtgtCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCa                                             >  1:219041/1‑140 (MQ=255)
                                                                                          aaTATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAgc                                >  1:265803/1‑140 (MQ=255)
                                                                                           atatCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAgcc                               >  2:221976/1‑140 (MQ=255)
                                                                                            tatCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAgccg                              >  2:240813/1‑140 (MQ=255)
                                                                                            tatCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAgccg                              >  2:591343/1‑140 (MQ=255)
                                                                                                 aaTACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTa                         <  2:265803/140‑1 (MQ=255)
                                                                                                   tACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAAt                       >  2:36884/1‑140 (MQ=255)
                                                                                                        tGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTc                  >  2:526080/1‑140 (MQ=255)
                                                                                                        tGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTc                  >  2:205472/1‑140 (MQ=255)
                                                                                                        tGGTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGCACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTc                  >  2:737439/1‑140 (MQ=255)
                                                                                                         ggTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTcc                 <  1:526080/140‑1 (MQ=255)
                                                                                                          gTATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCa                <  1:591343/140‑1 (MQ=255)
                                                                                                           tATTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAg               >  1:608183/1‑140 (MQ=255)
                                                                                                            aTTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGc              <  1:778607/140‑1 (MQ=255)
                                                                                                            aTTTGTCGTCAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGc              <  1:36884/140‑1 (MQ=255)
                                                                                                                gtcgtcAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACt          >  2:691933/1‑140 (MQ=255)
                                                                                                                   gtcAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTtat       <  1:221976/140‑1 (MQ=255)
                                                                                                                    tcAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTtata      >  2:664256/1‑140 (MQ=255)
                                                                                                                    tcAACGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTtata      >  1:246958/1‑140 (MQ=255)
                                                                                                                       aCGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGc   <  1:561803/140‑1 (MQ=255)
                                                                                                                        cGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCt  <  2:608183/140‑1 (MQ=255)
                                                                                                                        cGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCt  <  2:246958/140‑1 (MQ=255)
                                                                                                                        cGGGAAGGGGTAGGCCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCt  <  2:219041/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                             
GACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCT  >  NC_000913/1587291‑1587550

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