Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 986,671 A→G R104R (CGT→CGC ompF ← outer membrane porin 1a (Ia,b,F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,6710AG100.0% 307.3 / NA 87R104R (CGT→CGCompFouter membrane porin 1a (Ia,b,F)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (47/40);  total (47/40)

AGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATG  >  NC_000913/986535‑986806
                                                                                                                                        |                                                                                                                                       
agTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGcg                                                                                                                                       >  1:325442/1‑139 (MQ=255)
    aTCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGtt                                                                                                                                   >  2:475888/1‑139 (MQ=255)
                     caccacCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGc                                                                                                                  >  2:572788/1‑139 (MQ=255)
                       ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGt                                                                                                                >  2:808593/1‑139 (MQ=255)
                          ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAg                                                                                                             >  1:1462854/1‑139 (MQ=255)
                            aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcg                                                                                                           >  2:1604399/1‑139 (MQ=255)
                            aaaTTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAACCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcg                                                                                                           >  2:757969/1‑139 (MQ=255)
                                   ggCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCag                                                                                                    <  2:441471/139‑1 (MQ=255)
                                    gcagcaTATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCaga                                                                                                   <  1:475888/139‑1 (MQ=255)
                                        caTATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttg                                                                                               >  1:1321573/1‑139 (MQ=255)
                                             cGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAcc                                                                                          <  2:1357544/139‑1 (MQ=255)
                                               gtgtAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCt                                                                                        <  1:572788/139‑1 (MQ=255)
                                                 gtAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                                      <  2:1321573/139‑1 (MQ=255)
                                                 gtAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                                      <  2:810583/139‑1 (MQ=255)
                                                  tAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGa                                                                                     <  1:309694/139‑1 (MQ=255)
                                                     cccAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  <  2:543767/139‑1 (MQ=255)
                                                     cccAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  >  1:543767/1‑139 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGttt                                                                                                                                  >  1:1239710/1‑86 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGttt                                                                                                                                  <  2:1239710/86‑1 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGtt                                                                                                                              >  1:1476308/1‑90 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGtt                                                                                                                              <  2:1476308/90‑1 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCtt                                                                                                       >  2:961767/1‑113 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCtt                                                                                                       <  1:961767/113‑1 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCt                                                                                                        <  1:405884/112‑1 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCt                                                                                                        >  2:405884/1‑112 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCaga                                                                                                   >  1:1352138/1‑117 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCaga                                                                                                   <  2:1352138/117‑1 (MQ=255)
                                                          tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTatat                                                                             >  1:822386/1‑139 (MQ=255)
                                                               catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCGGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  1:808593/139‑1 (MQ=255)
                                                               catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  2:325442/139‑1 (MQ=255)
                                                                  aaaCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTg                                                                     >  1:447936/1‑139 (MQ=255)
                                                                    aCCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGAc                                                                   >  2:1360884/1‑139 (MQ=255)
                                                                        cacCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        >  2:349979/1‑130 (MQ=255)
                                                                        cacCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  1:349979/130‑1 (MQ=255)
                                                                        cacCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        >  2:266197/1‑130 (MQ=255)
                                                                        cacCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  1:266197/130‑1 (MQ=255)
                                                                         acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATaa                                                              >  1:991071/1‑139 (MQ=255)
                                                                            gTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCg                                                           >  1:1319808/1‑139 (MQ=255)
                                                                             tAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                              <  1:1564586/103‑1 (MQ=255)
                                                                             tAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                              >  2:1564586/1‑103 (MQ=255)
                                                                              aGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGt                                                         <  2:447936/139‑1 (MQ=255)
                                                                                   cGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGAt                                                    <  1:1288661/139‑1 (MQ=255)
                                                                                     gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  1:398354/117‑1 (MQ=255)
                                                                                     gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        >  2:398354/1‑117 (MQ=255)
                                                                                       cGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAccc                                                                                         >  1:529104/1‑98 (MQ=255)
                                                                                       cGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAccc                                                                                         <  2:529104/98‑1 (MQ=255)
                                                                                        gTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAcc                                                                                          <  1:14810/96‑1 (MQ=255)
                                                                                        gTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAcc                                                                                          >  2:14810/1‑96 (MQ=255)
                                                                                        gTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGaa                                               >  2:110072/1‑139 (MQ=255)
                                                                                         tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        >  2:507930/1‑113 (MQ=255)
                                                                                         tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  1:507930/113‑1 (MQ=255)
                                                                                          aaTCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCa                                                                 <  2:1361032/119‑1 (MQ=255)
                                                                                          aaTCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCa                                                                 >  1:1361032/1‑119 (MQ=255)
                                                                                                   aaCCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGtt                                                                                                                              >  2:819648/1‑49 (MQ=255)
                                                                                                   aaCCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGtt                                                                                                                              <  1:819648/49‑1 (MQ=255)
                                                                                                   aaCCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGt                                    >  1:827435/1‑139 (MQ=255)
                                                                                                      cAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcgc                                                                                                          >  1:770636/1‑66 (MQ=255)
                                                                                                      cAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcgc                                                                                                          <  2:770636/66‑1 (MQ=255)
                                                                                                      cAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  <  2:986796/90‑1 (MQ=255)
                                                                                                      cAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                                  >  1:986796/1‑90 (MQ=255)
                                                                                                      cAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTc                                 <  2:1319808/139‑1 (MQ=255)
                                                                                                        aCGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        >  2:623316/1‑98 (MQ=255)
                                                                                                        aCGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  1:623316/98‑1 (MQ=255)
                                                                                                         cGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAg                                                                                                             >  2:661876/1‑60 (MQ=255)
                                                                                                         cGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAg                                                                                                             <  1:661876/60‑1 (MQ=255)
                                                                                                             aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGAt                                                    >  1:164031/1‑113 (MQ=255)
                                                                                                             aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGAt                                                    <  2:164031/113‑1 (MQ=255)
                                                                                                                       aGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGAc                >  2:762015/1‑139 (MQ=255)
                                                                                                                         aCCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                              <  2:544486/59‑1 (MQ=255)
                                                                                                                         aCCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                              >  1:544486/1‑59 (MQ=255)
                                                                                                                         aCCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCAt                                                                >  1:1426734/1‑89 (MQ=255)
                                                                                                                         aCCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCAt                                                                <  2:1426734/89‑1 (MQ=255)
                                                                                                                         aCCAGCGAATGCCCGGCGCGTTTTGATACAAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTCCCCTGGACGTTATATTCCCACTGACCATAACCGGCCAGATCGGCATTCAATTGAGTTTCCCCTTTAAAACCAAGACgg              >  1:609823/1‑139 (MQ=255)
                                                                                                                            cgcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTAGAAGTTATATTCCCACTGACCATAACCGTTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCa           >  2:959142/1‑139 (MQ=255)
                                                                                                                             gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        >  2:281330/1‑77 (MQ=255)
                                                                                                                             gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  1:281330/77‑1 (MQ=255)
                                                                                                                             gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                 >  1:972826/1‑132 (MQ=255)
                                                                                                                             gcgAATGCCAGGCGCGTTTTGTTACCAGCTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                 <  2:972826/132‑1 (MQ=255)
                                                                                                                              cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                                      <  2:29786/46‑1 (MQ=255)
                                                                                                                              cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                                      >  1:29786/1‑46 (MQ=255)
                                                                                                                              cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        <  1:102710/76‑1 (MQ=255)
                                                                                                                              cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                        >  2:102710/1‑76 (MQ=255)
                                                                                                                              cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATa         >  1:1132473/1‑139 (MQ=255)
                                                                                                                                aaTGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCa                                                                                                     >  1:1236205/1‑45 (MQ=255)
                                                                                                                                aaTGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCa                                                                                                     <  2:1236205/45‑1 (MQ=255)
                                                                                                                                     cAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg  <  2:1132473/139‑1 (MQ=255)
                                                                                                                                     cAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg  >  2:1569738/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATG  >  NC_000913/986535‑986806

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: