Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,839,731 A→G T70A (ACC→GCC)  ascF → beta‑glucoside specific PTS enzyme IIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,839,7310AG98.5% 234.9 / ‑6.1 66T70A (ACC→GCC) ascFbeta‑glucoside specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (34/31);  total (35/31)
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

TCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATG  >  NC_000913/2839597‑2839865
                                                                                                                                      |                                                                                                                                      
tCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGtt                                                                                                                                   >  1:783992/1‑140 (MQ=255)
      aCTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCTGTTTCACaa                                                                                                                             <  2:1079455/140‑1 (MQ=255)
      aCTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACaa                                                                                                                             <  1:518915/140‑1 (MQ=255)
          tATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCt                                                                                                                         >  2:509693/1‑140 (MQ=255)
            tGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCttt                                                                                                                       >  1:229980/1‑140 (MQ=255)
               cgcgGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAg                                                                                                                    >  2:49158/1‑140 (MQ=255)
                    ttGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAgt                                                                                                               >  2:846242/1‑140 (MQ=255)
                                   aaaGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGcc                                                                                                >  1:766185/1‑140 (MQ=255)
                                   aaaGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGcc                                                                                                <  1:764060/140‑1 (MQ=255)
                                       atgatgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGga                                                                                            <  2:229980/140‑1 (MQ=255)
                                       atgatgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGga                                                                                            <  1:49158/140‑1 (MQ=255)
                                        tgatgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGgag                                                                                           >  2:135839/1‑140 (MQ=255)
                                        tgatgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGgag                                                                                           >  1:400671/1‑140 (MQ=255)
                                        tgatgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGgag                                                                                           >  2:424832/1‑140 (MQ=255)
                                           tgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGGAGTCGTCAGCCTGCTGCCGGGAGata                                                                                        >  1:765485/1‑140 (MQ=255)
                                           tgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGcc                                                                                                >  2:688883/1‑132 (MQ=255)
                                           tgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGcc                                                                                                <  1:688883/132‑1 (MQ=255)
                                           tgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGata                                                                                        >  1:1000085/1‑140 (MQ=255)
                                           tgCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGata                                                                                        >  1:188648/1‑140 (MQ=255)
                                            gCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGatat                                                                                       <  2:400671/140‑1 (MQ=255)
                                               cTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCGGGAAGTCGTCAGCCTGCTGcc                                                                                                <  2:252999/128‑1 (MQ=255)
                                               cTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCGGGAAGTCGTCAGCCTGCTGcc                                                                                                >  1:252999/1‑128 (MQ=255)
                                               cTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGcc                                                                                                <  1:20273/128‑1 (MQ=255)
                                               cTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGcc                                                                                                >  2:20273/1‑128 (MQ=255)
                                                 tATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGc                                                                                  >  1:678616/1‑140 (MQ=255)
                                                     gACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGc                                                                              >  2:748306/1‑140 (MQ=255)
                                                      aCAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGca                                                                             <  1:424832/140‑1 (MQ=255)
                                                      aCAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGca                                                                             <  2:1000085/140‑1 (MQ=255)
                                                      aCAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGca                                                                             <  1:1053492/140‑1 (MQ=255)
                                                          cccGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGc                                                                         >  2:1043528/1‑140 (MQ=255)
                                                                   aaaaaCCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTa                                                                <  2:188648/140‑1 (MQ=255)
                                                                         cATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTTACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCCCACCCGGGGGGAAAACCaa                                                         >  1:888463/1‑140 (MQ=255)
                                                                          aTCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCaa                                                         >  2:555601/1‑140 (MQ=255)
                                                                           tCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGt                                                                 <  1:1037114/131‑1 (MQ=255)
                                                                           tCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGt                                                                 >  2:1037114/1‑131 (MQ=255)
                                                                             cccGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACt                                                      <  1:748306/140‑1 (MQ=255)
                                                                                   gTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCt                                                <  1:135839/140‑1 (MQ=255)
                                                                                    tGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTa                                               <  1:509693/140‑1 (MQ=255)
                                                                                    tGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTa                                               <  1:555601/140‑1 (MQ=255)
                                                                                      cTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACg                                             >  2:681849/1‑140 (MQ=255)
                                                                                           cGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCa                                        <  2:678616/140‑1 (MQ=255)
                                                                                           cGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCa                                        <  1:681849/140‑1 (MQ=255)
                                                                                            gTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCAt                                       <  2:766185/140‑1 (MQ=255)
                                                                                            gTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCAt                                       <  1:630299/140‑1 (MQ=255)
                                                                                                      gTGACAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCgggg                             <  2:765485/140‑1 (MQ=255)
                                                                                                       tGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCAGTGGGTAAACCCAAACTCACGCTACGTCGTATTGGTGCggggg                            >  2:155268/1‑139 (MQ=255)
                                                                                                         aCAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATc                          >  2:1033071/1‑140 (MQ=255)
                                                                                                          cAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCAt                                       >  1:1073526/1‑126 (MQ=255)
                                                                                                          cAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCAt                                       <  2:1073526/126‑1 (MQ=255)
                                                                                                          cAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATcc                         >  2:217429/1‑140 (MQ=255)
                                                                                                          cAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGAGcc                         >  2:1020717/1‑140 (MQ=255)
                                                                                                          cAACCAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGAGATcc                         >  2:967638/1‑140 (MQ=255)
                                                                                                             ccAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCg                      <  1:217429/140‑1 (MQ=255)
                                                                                                              cAGTGTCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGa                     >  2:254144/1‑140 (MQ=255)
                                                                                                                  gtCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATgcg                 <  1:150294/140‑1 (MQ=255)
                                                                                                                  gtCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATgcg                 <  2:55318/140‑1 (MQ=255)
                                                                                                                  gtCAGGTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATgcg                 <  2:783992/140‑1 (MQ=255)
                                                                                                                       gTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGAt            >  2:713244/1‑140 (MQ=255)
                                                                                                                       gTGATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGAt            >  2:916948/1‑140 (MQ=255)
                                                                                                                         gATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATgcg                 <  2:799302/133‑1 (MQ=255)
                                                                                                                         gATTATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATgcg                 >  1:799302/1‑133 (MQ=255)
                                                                                                                           ttATCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGc        >  2:782243/1‑140 (MQ=255)
                                                                                                                             aTCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCAc      >  2:361216/1‑140 (MQ=255)
                                                                                                                             aTCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCAc      <  1:782243/140‑1 (MQ=255)
                                                                                                                             aTCGGCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCAc      <  1:1043528/140‑1 (MQ=255)
                                                                                                                                 gCAATGCCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATg  >  2:218699/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                      
TCACGCACTGTATGACGCGGTTGCGCTTTGTTATCAAAGATGATGCACTTATCGACAGCCCGACGTTAAAAACCATCCCCGGCGTGCTCGGCGTGGTACGTAGTGACAACCAGTGTCAGGTGATTATCGGCAATACCGTTTCACAAGCCTTTCAGGAAGTCGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATG  >  NC_000913/2839597‑2839865

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: