Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,311,986 A→C L255R (CTG→CGG)  ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,311,9860AC98.8% 303.0 / NA 82L255R (CTG→CGG) ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (42/39);  minor base G (0/1);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGA  >  NC_000913/2311752‑2312218
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
gTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTa                                                                                                                                                                                                                                      >  1:777283/1‑239 (MQ=255)
     cATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGAtgt                                                                                                                                                                                                                                 >  2:45625/1‑239 (MQ=255)
           tttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAg                                                                                                                                                                                                                           <  1:43227/239‑1 (MQ=255)
             tgttgAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGcg                                                                                                                                                                                                                         <  2:554679/239‑1 (MQ=255)
             tgttgAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGcg                                                                                                                                                                                                                         <  1:45625/239‑1 (MQ=255)
                tgAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGcgtcg                                                                                                                                                                                                                      >  1:217385/1‑239 (MQ=255)
                            tagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGacc                                                                                                                                                                                                           <  2:405308/238‑1 (MQ=255)
                            tagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGacc                                                                                                                                                                                                           <  1:1148968/238‑1 (MQ=255)
                            tagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGacc                                                                                                                                                                                                           >  1:768400/1‑238 (MQ=255)
                             agCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGaccac                                                                                                                                                                                                         <  1:657132/239‑1 (MQ=255)
                             agCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGaccac                                                                                                                                                                                                         <  1:707384/239‑1 (MQ=255)
                               cACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGaccac                                                                                                                                                                                                         <  2:782010/237‑1 (MQ=255)
                               cACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGaccac                                                                                                                                                                                                         >  1:782010/1‑237 (MQ=255)
                               cACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTAATGAGCAGCACGGTAGAAGTTGTTATCGTCGTATTTTAGAccacca                                                                                                                                                                                                       >  1:622143/1‑239 (MQ=255)
                                    acatcaacatATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTa                                                                                                                                                                                                  >  1:942966/1‑239 (MQ=255)
                                                        atCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGt                                                                                                                                                                              >  1:808606/1‑239 (MQ=255)
                                                        atCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGt                                                                                                                                                                              >  1:88755/1‑239 (MQ=255)
                                                          cTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTa                                                                                                                                                                            <  2:54280/239‑1 (MQ=255)
                                                                              cccAGGTTTTTACCTTTAGCCTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGGCGAACTGGTACTGAGCAACAGCTTCGAAGTTCGGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGtt                                                                                                                                                        <  2:646625/239‑1 (MQ=255)
                                                                              cccAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGtt                                                                                                                                                        <  2:942966/239‑1 (MQ=255)
                                                                                 aGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCAGGTTCGCCGTTACCGATGTAAACAGCGgggt                                                                                                                                                         >  2:615929/1‑235 (MQ=255)
                                                                                        tACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAg                                                                                                                                                                                                                           >  1:992201/1‑162 (MQ=255)
                                                                                        tACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAg                                                                                                                                                                                                                           <  2:992201/162‑1 (MQ=255)
                                                                                        tACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCg                                                                                                                                                             <  2:914894/224‑1 (MQ=255)
                                                                                        tACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCg                                                                                                                                                             >  1:914894/1‑224 (MQ=255)
                                                                                             ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAagcag                                                                                                                                                               >  2:80877/1‑217 (MQ=255)
                                                                                             ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAagcag                                                                                                                                                               <  1:80877/217‑1 (MQ=255)
                                                                                             ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGGACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGGGTAGGTTTCAGCACGGTCGCCGTTACCGATGTCAGCAGCGGTGTTCTGAGCATCAGTACg                                                                                                                                         >  1:975015/1‑239 (MQ=255)
                                                                                                     cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGa                                                                                                                                   >  1:962274/1‑237 (MQ=255)
                                                                                                     cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGa                                                                                                                                   >  1:381523/1‑237 (MQ=255)
                                                                                                          aaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGga                                                                                                                             <  2:768400/238‑1 (MQ=255)
                                                                                                                  gacggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGc                                                                                                                       <  2:962274/236‑1 (MQ=255)
                                                                                                                  gacggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAAATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTg                                                                                                                                     >  1:595883/1‑222 (MQ=255)
                                                                                                                  gacggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAAATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTg                                                                                                                                     <  2:595883/222‑1 (MQ=255)
                                                                                                                      gacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCg                                                                                                                >  1:165047/1‑239 (MQ=255)
                                                                                                                       acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCg                                                                                                                <  2:808606/238‑1 (MQ=255)
                                                                                                                        cgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGcac                                                                                                                     <  1:814282/232‑1 (MQ=255)
                                                                                                                           aGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTAt                                                                                                                                                                                                                   >  2:487789/1‑135 (MQ=255)
                                                                                                                           aGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTAAGCAGCCCGGTAGATGTTGTTAGCGTCGTAt                                                                                                                                                                                                                   <  1:487789/135‑1 (MQ=255)
                                                                                                                             aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTa                                                                                                                                                                            <  2:396200/172‑1 (MQ=255)
                                                                                                                             aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCCCCGTTa                                                                                                                                                                            >  1:396200/1‑172 (MQ=255)
                                                                                                                                  aGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGa                                                                                                         <  1:615929/234‑1 (MQ=255)
                                                                                                                                       aaCTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTa                                                                                               >  1:921940/1‑239 (MQ=255)
                                                                                                                                              aCTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATa                                                                                         <  2:921940/238‑1 (MQ=255)
                                                                                                                                                       cAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAg                                                                                                                                                                                                                                     >  1:738627/1‑89 (MQ=255)
                                                                                                                                                       cAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAg                                                                                                                                                                                                                                     <  2:738627/89‑1 (MQ=255)
                                                                                                                                                       cAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGt                                                                                                                                            <  2:485575/178‑1 (MQ=255)
                                                                                                                                                       cAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGt                                                                                                                                            >  1:485575/1‑178 (MQ=255)
                                                                                                                                                          cTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCg                                                                                                                                                             <  2:939626/158‑1 (MQ=255)
                                                                                                                                                          cTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCg                                                                                                                                                             <  2:181127/158‑1 (MQ=255)
                                                                                                                                                          cTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCg                                                                                                                                                             >  1:181127/1‑158 (MQ=255)
                                                                                                                                                          cTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCg                                                                                                                                                             >  1:939626/1‑158 (MQ=255)
                                                                                                                                                            tCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAAccgcc                                                                          >  1:157082/1‑239 (MQ=255)
                                                                                                                                                                       tgCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGtt                                                                                                                                       >  1:135526/1‑167 (MQ=255)
                                                                                                                                                                       tgCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGtt                                                                                                                                       <  2:135526/167‑1 (MQ=255)
                                                                                                                                                                            tGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATa                                                                                         <  1:93241/208‑1 (MQ=255)
                                                                                                                                                                            tGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATa                                                                                         >  2:93241/1‑208 (MQ=255)
                                                                                                                                                                                   ccAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTc                                                                                                  >  1:337382/1‑192 (MQ=255)
                                                                                                                                                                                   ccAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTc                                                                                                  <  2:337382/192‑1 (MQ=255)
                                                                                                                                                                                            ggAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGaccac                                                                                                                                                                                                         <  1:586088/80‑1 (MQ=255)
                                                                                                                                                                                            ggAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGaccac                                                                                                                                                                                                         >  2:586088/1‑80 (MQ=255)
                                                                                                                                                                                            ggAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTaa                                                                                              >  1:1090068/1‑187 (MQ=255)
                                                                                                                                                                                            ggAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTaa                                                                                              <  2:1090068/187‑1 (MQ=255)
                                                                                                                                                                                             gAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGGTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGCGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACg                                         >  1:643822/1‑239 (MQ=255)
                                                                                                                                                                                                 cTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAg                                                                                                                                                                                                                                     <  2:193830/47‑1 (MQ=255)
                                                                                                                                                                                                 cTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAg                                                                                                                                                                                                                                     >  1:193830/1‑47 (MQ=255)
                                                                                                                                                                                                          ttGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTg                                                                                                                                                     >  2:928786/1‑118 (MQ=255)
                                                                                                                                                                                                          ttGCGTTGTAGGTCTGGGTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTg                                                                                                                                                     <  1:928786/118‑1 (MQ=255)
                                                                                                                                                                                                                  tAGGTCTGGGTGTACTGAGCAGACCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTCGTTACGCCACt                    >  1:789428/1‑239 (MQ=255)
                                                                                                                                                                                                                          ggTGTACTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGCGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAg             >  2:156070/1‑238 (MQ=255)
                                                                                                                                                                                                                            tgtcaggAGCGGCCGGGTGGATGTTTTTGGCGTCTTTTTTAAGACCACCAGTGTGGGTTCCACCACGGTCGCGGTTACCGAGGTAAGCAGGGGTGTTCTGAGCATCAGTAGGTTTGGAGCTGGAGACCCCCCCACCGATCCCAAACCCTTCGTAATAATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGcc           <  2:157082/232‑1 (MQ=255)
                                                                                                                                                                                                                               aCTGAGCAGCCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGGGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTACGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTc        >  1:122529/1‑238 (MQ=255)
                                                                                                                                                                                                                                   agcagcCCGGTAGATGTTGTTAGCGTCGTATTTCa                                                                                                                                                                                                               >  1:104805/1‑35 (MQ=255)
                                                                                                                                                                                                                                   agcagcCCGGTAGATGTTGTTAGCGTCGTATTTCa                                                                                                                                                                                                               <  2:104805/35‑1 (MQ=255)
                                                                                                                                                                                                                                    gcagcCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTa                                                                                                                                                                                                  <  2:152843/47‑1 (MQ=255)
                                                                                                                                                                                                                                    gcagcCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTa                                                                                                                                                                                                  >  1:152843/1‑47 (MQ=255)
                                                                                                                                                                                                                                    gcagcCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATa                                                                                         <  1:29732/152‑1 (MQ=255)
                                                                                                                                                                                                                                    gcagcCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATa                                                                                         >  2:29732/1‑152 (MQ=255)
                                                                                                                                                                                                                                    gcagcCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGa  >  1:649346/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cagcCCGGTAGATGTTGTTAGCGTCGTATTTCAGacc                                                                                                                                                                                                           <  1:699270/37‑1 (MQ=255)
                                                                                                                                                                                                                                     cagcCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGga                                                                                                                             <  1:496651/115‑1 (MQ=255)
                                                                                                                                                                                                                                     cagcCCGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGga                                                                                                                             >  2:496651/1‑115 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
GTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAGGTAGATGTTGTTAGCGTCGTATTTCAGACCACCAGTGTAGGTTTCAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGA  >  NC_000913/2311752‑2312218

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