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% 313.9 / NA 87R104R (CGT→CGC) ompFouter membrane porin 1a (Ia,b,F)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (46/41);  total (46/41)

CGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATA  >  NC_000913/986541‑986799
                                                                                                                                  |                                                                                                                                
cGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGtttt                                                                                                                          >  1:332593/1‑139 (MQ=255)
       tGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCa                                                                                                                   >  2:985247/1‑139 (MQ=255)
                 ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGt                                                                                                         >  2:983869/1‑139 (MQ=255)
                 ccaccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGt                                                                                                         >  2:760060/1‑139 (MQ=255)
                    ccaAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAg                                                                                                      <  1:985247/139‑1 (MQ=255)
                     caAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                                     <  1:760060/139‑1 (MQ=255)
                     caAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                                     <  2:1210583/139‑1 (MQ=255)
                     caAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                                     >  2:195267/1‑139 (MQ=255)
                               cagcaTATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCagag                                                                                           >  1:1514331/1‑139 (MQ=255)
                                          tgtAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcgc                                                                                                   >  1:209035/1‑120 (MQ=255)
                                          tgtAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcgc                                                                                                   <  2:209035/120‑1 (MQ=255)
                                            tAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcgc                                                                                                   <  2:1280998/118‑1 (MQ=255)
                                            tAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgcgc                                                                                                   >  1:1280998/1‑118 (MQ=255)
                                            tAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGa                                                                              >  1:1083775/1‑139 (MQ=255)
                                              aCCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAc                                                                                                                     >  2:433845/1‑98 (MQ=255)
                                              aCCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAc                                                                                                                     <  1:433845/98‑1 (MQ=255)
                                               cccAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                           >  1:579108/1‑139 (MQ=255)
                                                    tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAc                                                                                                                     <  2:693151/92‑1 (MQ=255)
                                                    tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAc                                                                                                                     >  1:693151/1‑92 (MQ=255)
                                                    tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttg                                                                                        >  2:399641/1‑121 (MQ=255)
                                                    tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttg                                                                                        <  1:399641/121‑1 (MQ=255)
                                                    tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                               >  2:1162832/1‑130 (MQ=255)
                                                    tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                               <  2:555168/130‑1 (MQ=255)
                                                    tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                               <  1:1162832/130‑1 (MQ=255)
                                                    tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                               >  1:555168/1‑130 (MQ=255)
                                                    tGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTatat                                                                      <  1:1060019/139‑1 (MQ=255)
                                                        tcatAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACTCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                           <  2:1036078/130‑1 (MQ=255)
                                                        tcatAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                           >  1:1036078/1‑130 (MQ=255)
                                                         catAAACCACATCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 <  2:1083775/139‑1 (MQ=255)
                                                         catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 >  2:1244198/1‑139 (MQ=255)
                                                         catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 <  2:579108/139‑1 (MQ=255)
                                                         catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 <  1:983869/139‑1 (MQ=255)
                                                         catAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 <  2:332593/139‑1 (MQ=255)
                                                                   acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATaa                                                       >  2:1327982/1‑139 (MQ=255)
                                                                   acCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATaa                                                       <  1:195267/139‑1 (MQ=255)
                                                                          ttACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                               <  2:127691/108‑1 (MQ=255)
                                                                          ttACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTgg                                                                               >  1:127691/1‑108 (MQ=255)
                                                                               gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                                     >  2:1565189/1‑81 (MQ=255)
                                                                               gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                                     <  1:1565189/81‑1 (MQ=255)
                                                                               gCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCg                                           >  2:413671/1‑139 (MQ=255)
                                                                                 cGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAcc                                                                                   <  1:951201/97‑1 (MQ=255)
                                                                                 cGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAcc                                                                                   >  2:951201/1‑97 (MQ=255)
                                                                                   tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                               <  1:1165878/83‑1 (MQ=255)
                                                                                   tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                               >  2:1165878/1‑83 (MQ=255)
                                                                                   tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 >  2:625214/1‑113 (MQ=255)
                                                                                   tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 <  1:625214/113‑1 (MQ=255)
                                                                                          aaGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 >  2:1448817/1‑106 (MQ=255)
                                                                                          aaGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 <  1:1448817/106‑1 (MQ=255)
                                                                                             aaCCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                               <  1:1426709/73‑1 (MQ=255)
                                                                                             aaCCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                               >  2:1426709/1‑73 (MQ=255)
                                                                                                  aCGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAcc                                                                                                                    <  1:776312/47‑1 (MQ=255)
                                                                                                  aCGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTAcc                                                                                                                    >  2:776312/1‑47 (MQ=255)
                                                                                                   cGTCAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTcccc                       >  2:157249/1‑139 (MQ=255)
                                                                                                      cAGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCttt                    <  1:1327982/139‑1 (MQ=255)
                                                                                                       aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                               >  1:590558/1‑63 (MQ=255)
                                                                                                       aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                               <  2:590558/63‑1 (MQ=255)
                                                                                                       aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGaa                                        <  2:336532/118‑1 (MQ=255)
                                                                                                       aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGaa                                        >  1:336532/1‑118 (MQ=255)
                                                                                                       aGCGTATTTAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTa                   >  1:108549/1‑139 (MQ=255)
                                                                                                             tttAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 <  2:888800/87‑1 (MQ=255)
                                                                                                             tttAAGACCCGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCa                                                                 >  1:888800/1‑87 (MQ=255)
                                                                                                                     ccGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGTGTCAGCGCCTTCAGAgttgt                                                                                       >  1:301600/1‑57 (MQ=255)
                                                                                                                     ccGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGTGTCAGCGCCTTCAGAgttgt                                                                                       <  2:301600/57‑1 (MQ=255)
                                                                                                                     ccGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                                     <  2:897198/43‑1 (MQ=255)
                                                                                                                     ccGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAgc                                                                                                     >  1:897198/1‑43 (MQ=255)
                                                                                                                     ccGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCa                                                                                              <  2:1114326/50‑1 (MQ=255)
                                                                                                                     ccGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGAt                                   <  2:1389650/109‑1 (MQ=255)
                                                                                                                     ccGCGAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGAt                                   >  1:1389650/1‑109 (MQ=255)
                                                                                                                     ccGCGAATGCCAGGCACGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCa                                                                                              >  1:1114326/1‑50 (MQ=255)
                                                                                                                      cgcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                       <  2:68981/56‑1 (MQ=255)
                                                                                                                      cgcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAgttgt                                                                                       >  1:68981/1‑56 (MQ=255)
                                                                                                                       gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                           <  1:422051/67‑1 (MQ=255)
                                                                                                                       gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                           >  2:422051/1‑67 (MQ=255)
                                                                                                                       gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTcc                         >  2:690104/1‑117 (MQ=255)
                                                                                                                       gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTcc                         <  1:690104/117‑1 (MQ=255)
                                                                                                                       gcgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCAt   >  2:1638890/1‑139 (MQ=255)
                                                                                                                        cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCg                                                                                                          >  2:1573620/1‑35 (MQ=255)
                                                                                                                        cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCg                                                                                                          <  1:1573620/35‑1 (MQ=255)
                                                                                                                        cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                               <  1:862546/46‑1 (MQ=255)
                                                                                                                        cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTc                                                                                               >  2:862546/1‑46 (MQ=255)
                                                                                                                        cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAcc                                                                                   >  1:1301336/1‑58 (MQ=255)
                                                                                                                        cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTAcc                                                                                   <  2:1301336/58‑1 (MQ=255)
                                                                                                                        cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                           >  2:482828/1‑66 (MQ=255)
                                                                                                                        cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGt                                                                           <  1:482828/66‑1 (MQ=255)
                                                                                                                        cgAATGCCAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATa  >  1:635197/1‑139 (MQ=255)
                                                                                                                              ccAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTc                          <  2:88375/109‑1 (MQ=255)
                                                                                                                              ccAGGCGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTc                          >  1:88375/1‑109 (MQ=255)
                                                                                                                                  |                                                                                                                                
CGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAACCACACCGTAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATA  >  NC_000913/986541‑986799

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