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% 282.6 / NA 81G152G (GGC→GGT) bglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/42);  total (39/42)

AACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAG  >  NC_000913/3901631‑3901898
                                                                                                                                      |                                                                                                                                     
aaCCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCa                                                                                                                                   <  1:1255187/139‑1 (MQ=255)
 aCCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCag                                                                                                                                  >  2:1349977/1‑139 (MQ=255)
   cAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCagag                                                                                                                                <  1:204067/139‑1 (MQ=255)
          tAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                         <  1:32915/139‑1 (MQ=255)
               tACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGccc                                                                                                                    <  1:1175695/139‑1 (MQ=255)
                       aCCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTa                                                                                                            >  1:1206008/1‑139 (MQ=255)
                       aCCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTa                                                                                                            >  1:794141/1‑139 (MQ=255)
                         cAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAgg                                                                                                          <  2:1205704/139‑1 (MQ=255)
                         cAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAgg                                                                                                          >  2:262718/1‑139 (MQ=255)
                          aaCGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGt                                                                                                         <  2:833978/139‑1 (MQ=255)
                             gTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGaa                                                                                                      <  1:730106/139‑1 (MQ=255)
                              ttACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAAc                                                                                                     >  1:406089/1‑139 (MQ=255)
                                aCCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTTGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCa                                                                                                   >  2:358780/1‑139 (MQ=255)
                                      catcatCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGagag                                                                                             >  1:323281/1‑139 (MQ=255)
                                       atcatcACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCa                                                                                                                                   >  1:1041403/1‑100 (MQ=255)
                                       atcatcACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCa                                                                                                                                   <  2:1041403/100‑1 (MQ=255)
                                         catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGa                                                                                                                <  2:184912/117‑1 (MQ=255)
                                         catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGa                                                                                                                >  1:184912/1‑117 (MQ=255)
                                         catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                          >  1:257183/1‑139 (MQ=255)
                                         catcaCAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGGGAACCATGAGAGGCg                                                                                          <  2:978772/139‑1 (MQ=255)
                                                    cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGaa                                                                                                      >  1:1449059/1‑116 (MQ=255)
                                                    cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGaa                                                                                                      <  2:1449059/116‑1 (MQ=255)
                                                    cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAAcc                                                                               >  1:1102985/1‑139 (MQ=255)
                                                      tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                          >  1:213239/1‑126 (MQ=255)
                                                      tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                          <  2:213239/126‑1 (MQ=255)
                                                      tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAg                                                                             >  1:1048392/1‑139 (MQ=255)
                                                       cAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGa                                                                            <  1:285494/139‑1 (MQ=255)
                                                       cAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGa                                                                            <  1:262718/139‑1 (MQ=255)
                                                        aaCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                           <  1:716452/139‑1 (MQ=255)
                                                        aaCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                           <  2:1340493/139‑1 (MQ=255)
                                                        aaCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                           <  2:1102985/139‑1 (MQ=255)
                                                         aCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCa                                                                                                                                   >  1:471003/1‑82 (MQ=255)
                                                         aCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCa                                                                                                                                   <  2:471003/82‑1 (MQ=255)
                                                         aCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCCGAAc                                                                          >  2:1103281/1‑139 (MQ=255)
                                                              tttGCCGTTTTCGTTGTAGACCCGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTaa                                                                     >  2:1316074/1‑139 (MQ=255)
                                                              tttGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                     <  1:86396/91‑1 (MQ=255)
                                                              tttGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                     >  2:86396/1‑91 (MQ=255)
                                                                tGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATa                                                                   <  2:1373992/139‑1 (MQ=255)
                                                                tGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGA‑TTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAg                                                                  <  1:773551/139‑1 (MQ=255)
                                                                 gCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAg                                                                  <  2:1181285/139‑1 (MQ=255)
                                                                 gCCGTTTTCGTTCTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGCGAACCGATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAg                                                                  <  2:239941/139‑1 (MQ=255)
                                                                     ttttCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAg                                                                                                           <  2:1372106/94‑1 (MQ=255)
                                                                     ttttCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAg                                                                                                           >  1:1372106/1‑94 (MQ=255)
                                                                       ttCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTg                                                            <  2:900233/139‑1 (MQ=255)
                                                                       ttCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTg                                                            <  1:1453723/139‑1 (MQ=255)
                                                                                      ttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAtt                                             >  1:320020/1‑139 (MQ=255)
                                                                                        gttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                          >  1:1066576/1‑92 (MQ=255)
                                                                                        gttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                          <  2:1066576/92‑1 (MQ=255)
                                                                                              aCTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTg                                                            >  2:976358/1‑116 (MQ=255)
                                                                                              aCTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTg                                                            <  1:976358/116‑1 (MQ=255)
                                                                                              aCTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCttt                                     <  1:1349977/139‑1 (MQ=255)
                                                                                               cTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAgg                                                                                                          <  2:513778/69‑1 (MQ=255)
                                                                                               cTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAgg                                                                                                          >  1:513778/1‑69 (MQ=255)
                                                                                                 gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAg                                                   <  1:1440378/122‑1 (MQ=255)
                                                                                                 gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAg                                                   >  2:1440378/1‑122 (MQ=255)
                                                                                                  cAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTcaca                                 >  2:1158142/1‑139 (MQ=255)
                                                                                                     aTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAAt                              <  2:1206008/139‑1 (MQ=255)
                                                                                                            cAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAtt                                                                                                              <  2:507084/52‑1 (MQ=255)
                                                                                                            cAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGAtt                                                                                                              >  1:507084/1‑52 (MQ=255)
                                                                                                            cAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTca                                   >  2:1164996/1‑127 (MQ=255)
                                                                                                            cAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTca                                   <  1:1164996/127‑1 (MQ=255)
                                                                                                               ccGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAAt                    >  1:1063694/1‑139 (MQ=255)
                                                                                                                        tCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCagag           >  1:611663/1‑139 (MQ=255)
                                                                                                                         cGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaa                      >  2:915506/1‑127 (MQ=255)
                                                                                                                         cGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaa                      <  1:915506/127‑1 (MQ=255)
                                                                                                                         cGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGt          <  2:257183/139‑1 (MQ=255)
                                                                                                                         cGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGt          <  2:1048392/139‑1 (MQ=255)
                                                                                                                             aCCAAAGCGACCCAGAGAACCAATCGCCCATGCTTTAGGTGCACCATGAGAGGCGGTCCCCCAACAAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGGCATg      >  2:999365/1‑139 (MQ=255)
                                                                                                                             aCCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATg      >  2:13388/1‑139 (MQ=255)
                                                                                                                               cAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCa              >  1:1121077/1‑129 (MQ=255)
                                                                                                                               cAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCa              <  2:1121077/129‑1 (MQ=255)
                                                                                                                               cAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGga    >  2:1180467/1‑139 (MQ=255)
                                                                                                                                 aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                           <  1:38639/66‑1 (MQ=255)
                                                                                                                                 aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                           >  2:38639/1‑66 (MQ=255)
                                                                                                                                 aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCag             <  2:590459/128‑1 (MQ=255)
                                                                                                                                 aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCag             >  1:590459/1‑128 (MQ=255)
                                                                                                                                 aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaag  >  1:1463457/1‑139 (MQ=255)
                                                                                                                                  aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAgg                                                                                            >  2:795120/1‑48 (MQ=255)
                                                                                                                                  aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAgg                                                                                            <  1:795120/48‑1 (MQ=255)
                                                                                                                                  aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaa                      >  2:789413/1‑118 (MQ=255)
                                                                                                                                  aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaa                      <  1:789413/118‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
AACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAG  >  NC_000913/3901631‑3901898

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: ‑