Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,901,765 G→A G152G (GGC→GGT)  bglH ← carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,901,7650GA100.0% 265.9 / NA 76G152G (GGC→GGT) bglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/40);  total (36/40)

ACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGA  >  NC_000913/3901632‑3901899
                                                                                                                                     |                                                                                                                                      
aCCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCag                                                                                                                                   <  1:512453/139‑1 (MQ=255)
 ccAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCaga                                                                                                                                  <  1:43303/139‑1 (MQ=255)
         tAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                          <  2:287291/139‑1 (MQ=255)
               aCTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCa                                                                                                                    >  1:147913/1‑139 (MQ=255)
                   ctgACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAt                                                                                                                >  1:1123964/1‑139 (MQ=255)
                   ctgACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAt                                                                                                                >  1:1217511/1‑139 (MQ=255)
                    tgACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGCGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAtt                                                                                                               <  1:932958/139‑1 (MQ=255)
                      aCCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTa                                                                                                             >  2:1196823/1‑139 (MQ=255)
                       ccAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAg                                                                                                            <  1:455745/139‑1 (MQ=255)
                        cAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAgg                                                                                                           <  2:1217511/139‑1 (MQ=255)
                         aaCGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGt                                                                                                          <  2:851111/139‑1 (MQ=255)
                                      atcatcACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAgg                                                                                             <  1:1201506/139‑1 (MQ=255)
                                        catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCag                                                                                                                                   >  1:626610/1‑99 (MQ=255)
                                        catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCag                                                                                                                                   <  2:626610/99‑1 (MQ=255)
                                             caACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTccc                                                                                      >  1:70498/1‑139 (MQ=255)
                                              aaCGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTcccc                                                                                     >  1:599712/1‑139 (MQ=255)
                                              aaCGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTcccc                                                                                     >  1:236749/1‑139 (MQ=255)
                                                cGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                          >  2:159149/1‑100 (MQ=255)
                                                cGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                          <  1:159149/100‑1 (MQ=255)
                                                   cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                           <  2:829577/128‑1 (MQ=255)
                                                   cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCCCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                           >  1:829577/1‑128 (MQ=255)
                                                     tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGaa                                                                                                       <  1:400369/114‑1 (MQ=255)
                                                     tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGaa                                                                                                       >  2:400369/1‑114 (MQ=255)
                                                     tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAg                                                                              <  2:147913/139‑1 (MQ=255)
                                                        aCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAAc                                                                           >  1:134261/1‑139 (MQ=255)
                                                            gTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTa                                                                       >  1:996289/1‑139 (MQ=255)
                                                               tGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATa                                                                    <  1:1196823/139‑1 (MQ=255)
                                                                 ccGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGc                                                                  <  2:331024/139‑1 (MQ=255)
                                                                    ttttCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGt                                                               >  2:542159/1‑139 (MQ=255)
                                                                               gACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAg                                                    >  2:1062840/1‑139 (MQ=255)
                                                                                    gttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAtt                                                                                                               <  1:658305/75‑1 (MQ=255)
                                                                                    gttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAtt                                                                                                               >  2:658305/1‑75 (MQ=255)
                                                                                     ttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTccc                                                                                      <  1:582581/99‑1 (MQ=255)
                                                                                     ttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTccc                                                                                      >  2:582581/1‑99 (MQ=255)
                                                                                     ttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAtt                                              >  2:642260/1‑139 (MQ=255)
                                                                                       gttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGt                                                                                                          >  1:38856/1‑77 (MQ=255)
                                                                                       gttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGt                                                                                                          <  2:38856/77‑1 (MQ=255)
                                                                                       gttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                           >  2:559007/1‑92 (MQ=255)
                                                                                       gttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                           <  1:559007/92‑1 (MQ=255)
                                                                                                gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGccc                                                                                                                     <  2:473142/57‑1 (MQ=255)
                                                                                                gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGccc                                                                                                                     >  1:473142/1‑57 (MQ=255)
                                                                                                gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAAc                                                                                 <  1:318612/93‑1 (MQ=255)
                                                                                                gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAAc                                                                                 >  2:318612/1‑93 (MQ=255)
                                                                                                gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACt                                                        >  2:1043988/1‑118 (MQ=255)
                                                                                                gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACt                                                        <  1:1043988/118‑1 (MQ=255)
                                                                                                 cAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTg                                                                                                         <  1:716652/68‑1 (MQ=255)
                                                                                                 cAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTg                                                                                                         >  2:716652/1‑68 (MQ=255)
                                                                                                  aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGt                                                                                                          >  2:481434/1‑66 (MQ=255)
                                                                                                  aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGt                                                                                                          <  1:481434/66‑1 (MQ=255)
                                                                                                  aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTg                                                             >  1:1248666/1‑111 (MQ=255)
                                                                                                  aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTg                                                             <  2:1248666/111‑1 (MQ=255)
                                                                                                  aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACaa                                 >  2:571523/1‑139 (MQ=255)
                                                                                                  aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACaa                                 >  1:1055306/1‑139 (MQ=255)
                                                                                                  aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGAGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCAAATTTCATCTCTCACaa                                 >  1:609474/1‑139 (MQ=255)
                                                                                                     tCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAAtt                              <  2:628397/139‑1 (MQ=255)
                                                                                                             gCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaaa                      <  2:194567/139‑1 (MQ=255)
                                                                                                             gCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaaa                      <  2:1055306/139‑1 (MQ=255)
                                                                                                              ccGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                          <  2:686256/38‑1 (MQ=255)
                                                                                                              ccGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                          >  1:686256/1‑38 (MQ=255)
                                                                                                                   gTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTc                <  2:996289/139‑1 (MQ=255)
                                                                                                                   gTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTc                >  1:248735/1‑139 (MQ=255)
                                                                                                                    tATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAtt                                                                                                               >  2:126288/1‑43 (MQ=255)
                                                                                                                    tATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAtt                                                                                                               <  1:126288/43‑1 (MQ=255)
                                                                                                                    tATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCa               <  2:1198125/139‑1 (MQ=255)
                                                                                                                     aTTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCag              <  1:571523/139‑1 (MQ=255)
                                                                                                                      ttCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAg                                                    <  2:123888/100‑1 (MQ=255)
                                                                                                                      ttCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAg                                                    >  1:123888/1‑100 (MQ=255)
                                                                                                                            aCCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTcccc                                                                                     <  1:924658/61‑1 (MQ=255)
                                                                                                                            aCCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTcccc                                                                                     >  2:924658/1‑61 (MQ=255)
                                                                                                                            aCCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATg       <  1:542159/139‑1 (MQ=255)
                                                                                                                              cAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATg       >  1:84048/1‑137 (MQ=255)
                                                                                                                              cAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATg       <  2:84048/137‑1 (MQ=255)
                                                                                                                              cAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGga     <  1:642260/139‑1 (MQ=255)
                                                                                                                              cAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGga     <  2:134261/139‑1 (MQ=255)
                                                                                                                                aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaag   >  2:332303/1‑139 (MQ=255)
                                                                                                                                 aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaaga  <  2:248735/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
ACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGA  >  NC_000913/3901632‑3901899

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