Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,964,697 G→A L118L (CTC→CTT)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,964,6970GA100.0% 253.8 / NA 79L118L (CTC→CTT) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/40);  total (39/40)

ACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAAC  >  NC_000913/1964565‑1964832
                                                                                                                                    |                                                                                                                                       
aCCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAg                                                                                                                                   >  1:1462220/1‑139 (MQ=255)
               gagaACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCg                                                                                                                    <  2:495538/139‑1 (MQ=255)
                agaACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGa                                                                                                                   >  2:795305/1‑139 (MQ=255)
                agaACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGa                                                                                                                   >  1:1382515/1‑139 (MQ=255)
                    cAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGaccac                                                                                                               >  2:281192/1‑139 (MQ=255)
                       aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCtt                                                                                                            >  2:599390/1‑139 (MQ=255)
                       aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCtt                                                                                                            >  1:758921/1‑139 (MQ=255)
                       aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCtt                                                                                                            >  1:307592/1‑139 (MQ=255)
                           cgcACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCt                                                                                                        >  2:1462791/1‑139 (MQ=255)
                                  aCTTCAGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTgccgccg                                                                                                 <  2:609597/139‑1 (MQ=255)
                                     tCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGccgccgccg                                                                                              <  2:758921/139‑1 (MQ=255)
                                      cTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCcgccgccgc                                                                                             <  1:599390/139‑1 (MQ=255)
                                        gCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcgc                                                                                           <  2:1481507/139‑1 (MQ=255)
                                            tACGCCCGGCACCTTTGGTAATGACCATAAAGTCGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGt                                                                                       >  1:1559412/1‑139 (MQ=255)
                                             aCGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTa                                                                                      <  2:1462220/139‑1 (MQ=255)
                                              cGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTAt                                                                                     >  2:1351958/1‑139 (MQ=255)
                                                 ccGGCACCTTTGGTAATGACCATAAAGTCGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGcc                                                                                  <  2:1559412/139‑1 (MQ=255)
                                                  cGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                                 <  2:642860/139‑1 (MQ=255)
                                                  cGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                                 <  1:281192/139‑1 (MQ=255)
                                                      aCCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTcc                                                                             >  2:1433283/1‑139 (MQ=255)
                                                               aaTGACCATAAAGTTGATGATCACGAGAATGACAAACATCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcg                                                                                            <  1:1311888/115‑1 (MQ=255)
                                                               aaTGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAgg                                                                                                                                      >  1:662467/1‑73 (MQ=255)
                                                               aaTGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAgg                                                                                                                                      <  2:662467/73‑1 (MQ=255)
                                                               aaTGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcg                                                                                            >  2:1311888/1‑115 (MQ=255)
                                                                 tGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcg                                                                                            <  1:1421859/113‑1 (MQ=255)
                                                                 tGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcg                                                                                            >  2:1421859/1‑113 (MQ=255)
                                                                 tGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGAt                                                                  >  1:1500619/1‑139 (MQ=255)
                                                                  gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTgccgccgcc                                                                                               <  1:170410/109‑1 (MQ=255)
                                                                  gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTgccgccgcc                                                                                               >  2:170410/1‑109 (MQ=255)
                                                                  gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATa                                                                 >  1:1035488/1‑139 (MQ=255)
                                                                  gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATa                                                                 <  2:728403/139‑1 (MQ=255)
                                                                  gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATa                                                                 >  1:728403/1‑139 (MQ=255)
                                                                  gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATa                                                                 <  1:795305/139‑1 (MQ=255)
                                                                   aCCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGaa                                                                                                                                    >  1:639665/1‑71 (MQ=255)
                                                                   aCCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGaa                                                                                                                                    <  2:639665/71‑1 (MQ=255)
                                                                   aCCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATa                                                                 >  2:151484/1‑138 (MQ=255)
                                                                   aCCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATa                                                                 <  1:151484/138‑1 (MQ=255)
                                                                         aaGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTg                                                          <  1:1351958/139‑1 (MQ=255)
                                                                             tgatgaTCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATAc                                                                <  1:194003/129‑1 (MQ=255)
                                                                             tgatgaTCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATAc                                                                >  2:194003/1‑129 (MQ=255)
                                                                              gatgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGacc                                                                                                                 >  2:707785/1‑79 (MQ=255)
                                                                              gatgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGacc                                                                                                                 <  1:707785/79‑1 (MQ=255)
                                                                               atgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAAt                                                                     >  2:825420/1‑122 (MQ=255)
                                                                               atgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAAt                                                                     <  1:825420/122‑1 (MQ=255)
                                                                               atgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCa                                                    >  2:1061736/1‑139 (MQ=255)
                                                                               atgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCa                                                    >  1:1003892/1‑139 (MQ=255)
                                                                                   tCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCTCCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGtt                                                <  1:1061736/139‑1 (MQ=255)
                                                                                      cGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGa                                                                                                                            <  2:451170/60‑1 (MQ=255)
                                                                                      cGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGa                                                                                                                            >  1:451170/1‑60 (MQ=255)
                                                                                      cGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                                 <  2:1471317/103‑1 (MQ=255)
                                                                                      cGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                                 <  1:1533426/103‑1 (MQ=255)
                                                                                      cGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                                 >  1:1471317/1‑103 (MQ=255)
                                                                                      cGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                                 >  2:1533426/1‑103 (MQ=255)
                                                                                      cGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATAc                                                                <  1:1534185/120‑1 (MQ=255)
                                                                                      cGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATAc                                                                >  2:1534185/1‑120 (MQ=255)
                                                                                        agaATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTg                                           <  2:1035488/139‑1 (MQ=255)
                                                                                          aaTGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGcc                                         <  2:1382515/139‑1 (MQ=255)
                                                                                            tGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTg                                           <  2:1157211/135‑1 (MQ=255)
                                                                                            tGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTCCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTg                                           >  1:1157211/1‑135 (MQ=255)
                                                                                                 aaCACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCgg                                                                                        <  2:621421/85‑1 (MQ=255)
                                                                                                 aaCACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCgg                                                                                        >  1:621421/1‑85 (MQ=255)
                                                                                                       acGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCTGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCg                            <  1:1462791/139‑1 (MQ=255)
                                                                                                         gATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTg                          <  2:4653/139‑1 (MQ=255)
                                                                                                           tACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAAc                                                                                                                          <  2:1013354/41‑1 (MQ=255)
                                                                                                           tACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAAc                                                                                                                          >  1:1013354/1‑41 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGa                                                                                                                   <  1:783934/43‑1 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGa                                                                                                                   >  2:783934/1‑43 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTgccgccgc                                                                                                <  2:355519/62‑1 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTgccgccgc                                                                                                >  1:355519/1‑62 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTa                                                                                      >  2:1229851/1‑72 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTa                                                                                      <  1:1229851/72‑1 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATg                                                                                    >  2:1408007/1‑74 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATg                                                                                    <  1:1408007/74‑1 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                                 >  1:765204/1‑77 (MQ=255)
                                                                                                                aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                                 <  2:765204/77‑1 (MQ=255)
                                                                                                                  aGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGac                                                                                                                  >  2:1222991/1‑42 (MQ=255)
                                                                                                                  aGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGac                                                                                                                  <  1:1222991/42‑1 (MQ=255)
                                                                                                                    cGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAg               >  1:868819/1‑139 (MQ=255)
                                                                                                                                 aaCAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAAc  <  1:1433283/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
ACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAATGGTCGGAAAC  >  NC_000913/1964565‑1964832

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