Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,382,820 T→C 100% Y292H (TAT→CAT)  ompG → outer membrane porin G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,382,8200TC99.1% 353.1 / ‑4.7 106Y292H (TAT→CAT) ompGouter membrane porin G
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (44/61);  total (45/61)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

Only 100 of 108 total aligned reads displayed.

TTTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATA  >  NC_000913/1382724‑1382916
                                                                                                |                                                                                                
ttttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATg                                                                                               >  2:720892/1‑100 (MQ=255)
ttttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCAAGACGAAGGCGACAGTGATAAATTCCATCATg                                                                                               >  1:692595/1‑100 (MQ=255)
  ttACGGTTATGATTTTCATAACGGACTTTCCGTTTCGCTTGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCa                                                                                             <  1:1965894/100‑1 (MQ=255)
  ttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATgta                                                                                             <  1:2062651/100‑3 (MQ=255)
  ttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCa                                                                                             <  1:2177944/100‑1 (MQ=255)
   tACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAg                                                                                            >  2:2414250/1‑100 (MQ=255)
    aCGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAgg                                                                                           <  1:1830726/100‑1 (MQ=255)
     cGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAgg                                                                                           <  2:2443415/99‑1 (MQ=255)
     cGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgt                                                                                          <  1:2604233/100‑1 (MQ=255)
     cGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAAGCGACAGTGATAAATTCCATCATGCAGgt                                                                                          >  1:466467/1‑100 (MQ=255)
      ggTTATGATTTCCAGAACGTACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtg                                                                                         <  1:672165/100‑1 (MQ=255)
      ggTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtg                                                                                         >  2:2548994/1‑100 (MQ=255)
      ggTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtg                                                                                         <  1:1415031/100‑1 (MQ=255)
      ggTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTACATCATGCAGgt                                                                                          >  1:397365/1‑99 (MQ=255)
       gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtgt                                                                                        >  1:44958/1‑100 (MQ=255)
       gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtgt                                                                                        >  2:2445046/1‑100 (MQ=255)
       gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGACGACAGTGATAAATTCCATCATGCAGgtgt                                                                                        >  1:1167869/1‑100 (MQ=255)
        ttATGTTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACTCGTTTGAGTGGCATGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTc                                                                                       <  1:1850757/100‑1 (MQ=255)
        ttATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTc                                                                                       >  1:2519757/1‑100 (MQ=255)
         tATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGACCACAGTGATAAAGTCCATCATGCACGTGTc                                                                                       >  2:2322630/1‑99 (MQ=255)
           tGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTCGAGTGGCAGGATCACGACGAAGGCGACAGTGTTAAATTCCATCATGCAGGTGTCgg                                                                                     <  2:1941961/99‑1 (MQ=255)
              tttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTa                                                                                 >  2:1291014/1‑100 (MQ=255)
                tCCAGAACGGACTTTCCGTTTCGCTGGCATACGCGTTTGAGTGGCAGGATCACGACGAAGACGACAGTGATAAATTCAATCATGCAGGTGTCGGCGTaaa                                                                               >  1:348349/1‑100 (MQ=255)
                 ccAGATCGGATTTTCCGTTTGGCTGTAATACGCGTTTGAGTGGCAGGATCACGACGATTGTGACAGTGATAAATTCCATCATGCATGTGTGGGAGTaaa                                                                               <  2:2095949/99‑1 (MQ=255)
                 ccAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAt                                                                              >  2:127174/1‑100 (MQ=255)
                  cAGAACGGACTTTCCGTTTCGTTGGAATACGCGTTTGAGTGGCAGGATCAAGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAtt                                                                             <  1:262654/100‑1 (MQ=255)
                    gAACGGACTTTCCGTTTCGCTGGAATACGCGAATGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCCGGTGTCGGCGTAAATTAc                                                                           >  2:1682455/1‑100 (MQ=255)
                     aaCGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACt                                                                          <  1:2208104/100‑1 (MQ=255)
                      aCGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTc                                                                         >  1:1239347/1‑100 (MQ=255)
                      aCGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTc                                                                         <  1:916488/100‑1 (MQ=255)
                         gACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGtt                                                                      >  2:2308035/1‑100 (MQ=255)
                           cTTTCCGTTTCGCTGGAATACGCTTTTGAGTGGCAGCATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCg                                                                        <  2:867340/96‑1 (MQ=255)
                           cTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCt                                                                    <  1:720896/100‑1 (MQ=255)
                           cTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCt                                                                    >  1:2179350/1‑100 (MQ=255)
                            tttCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCt                                                                    >  2:1531471/1‑99 (MQ=255)
                              tCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGAt                                                                 >  2:485733/1‑100 (MQ=255)
                                   ttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGAt                                                                 <  2:1254530/95‑1 (MQ=255)
                                     cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                          <  1:1031744/100‑1 (MQ=255)
                                     cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                          >  1:2266233/1‑100 (MQ=255)
                                      gCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGATAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCt                                                         <  2:2130500/100‑1 (MQ=255)
                                        tGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTaa                                                       <  1:3251/100‑1 (MQ=255)
                                         ggAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTaaa                                                      >  1:97516/1‑100 (MQ=255)
                                              aCGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGcc                                                 <  1:394469/100‑1 (MQ=255)
                                              aCGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGcc                                                 >  1:79914/1‑100 (MQ=255)
                                              aCGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGcc                                                 >  1:83028/1‑100 (MQ=255)
                                               cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                <  1:635721/100‑1 (MQ=255)
                                               cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                >  2:2632363/1‑100 (MQ=255)
                                                    ttGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgc                                           <  1:1071343/100‑1 (MQ=255)
                                                    ttGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgc                                           <  1:991473/100‑1 (MQ=255)
                                                     tGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcg                                          <  2:238009/100‑1 (MQ=255)
                                                       aGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggc                                        >  1:202780/1‑100 (MQ=255)
                                                        gTGGCAGGATCACGACGAAGGCGACAGTGCTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCAGATAATGGGCTAAATTGCCGGATgcgg                                         >  1:201791/1‑98 (MQ=255)
                                                        gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggcg                                       >  1:1152542/1‑100 (MQ=255)
                                                        gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggcg                                       <  1:949915/100‑1 (MQ=255)
                                                        gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggcg                                       <  1:1299540/100‑1 (MQ=255)
                                                         tGGCAGTATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGTATGCGgcgc                                      <  1:1419385/100‑1 (MQ=255)
                                                         tGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGgcgc                                      <  2:44958/100‑1 (MQ=255)
                                                          ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                     >  2:1440812/1‑100 (MQ=255)
                                                          ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                     <  2:364911/100‑1 (MQ=255)
                                                           gCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                     <  1:538408/99‑1 (MQ=255)
                                                            cAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                   <  1:2524417/100‑1 (MQ=255)
                                                            cAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                   <  1:512406/100‑1 (MQ=255)
                                                            cAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                   <  2:615352/100‑1 (MQ=255)
                                                             aGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGt                                  >  1:528712/1‑100 (MQ=255)
                                                             aGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGt                                  <  1:1896558/100‑1 (MQ=255)
                                                              ggATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTa                                 >  1:165198/1‑100 (MQ=255)
                                                                aTCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACt                               <  1:408726/100‑1 (MQ=255)
                                                                 tCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                              <  2:675923/100‑1 (MQ=255)
                                                                  cACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACttt                             >  1:1231349/1‑100 (MQ=255)
                                                                   acgacgAAGGCGACATTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTa                            <  1:1024285/100‑1 (MQ=255)
                                                                     gacgaAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTAt                           >  2:1481038/1‑99 (MQ=255)
                                                                        gaAGGCGACAGTCATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCTCTAGTACTTTATCCGa                       <  2:2010737/100‑1 (MQ=255)
                                                                          aGGCTACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                     <  1:2881989/100‑1 (MQ=255)
                                                                          aGGCGACAGTGATAAATTCCATCATGCAGGTGTTGGCGTAAATTATTTGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                     <  1:851880/100‑1 (MQ=255)
                                                                          aGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAACTTGCCGGATGCGGCGCGAGTACTTTATCCGATc                     <  2:1815982/100‑1 (MQ=255)
                                                                          aGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                     <  2:2449505/100‑1 (MQ=255)
                                                                           ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                    >  2:2714448/1‑100 (MQ=255)
                                                                           ggCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                    >  2:2779457/1‑100 (MQ=255)
                                                                            gCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCta                   <  2:1548034/100‑1 (MQ=255)
                                                                              gACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCtata                 >  2:1981213/1‑100 (MQ=255)
                                                                                cAGTGATAAATTCCATCATGTAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTTTCCGATCTATaaa               <  2:2610132/100‑1 (MQ=255)
                                                                                cAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATaaa               <  1:464220/100‑1 (MQ=255)
                                                                                cAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATaaa               <  1:2430054/100‑1 (MQ=255)
                                                                                 aGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAGATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAAt              <  1:1081498/100‑1 (MQ=255)
                                                                                  gTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg             >  2:1460325/1‑100 (MQ=255)
                                                                                  gTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg             <  1:2646755/100‑1 (MQ=255)
                                                                                    gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa           <  1:289313/100‑1 (MQ=255)
                                                                                     aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCTGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg          <  2:1431835/100‑1 (MQ=255)
                                                                                     aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg          >  2:86280/1‑100 (MQ=255)
                                                                                      tAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAgg         >  1:382021/1‑100 (MQ=255)
                                                                                        aaTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGcc       <  1:1230865/100‑1 (MQ=255)
                                                                                         aTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCg      <  2:1152531/100‑1 (MQ=255)
                                                                                         aTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCg      >  2:627476/1‑100 (MQ=255)
                                                                                          ttCCATCATGCAGTTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCgg     <  1:2179056/100‑1 (MQ=255)
                                                                                          ttCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCgg     <  2:2495752/100‑1 (MQ=255)
                                                                                          ttCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCgg     <  1:447398/100‑1 (MQ=255)
                                                                                            ccATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt   <  1:2399353/100‑1 (MQ=255)
                                                                                            ccATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt   >  2:566264/1‑100 (MQ=255)
                                                                                            ccATCATCCAGGTCTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt   >  2:500375/1‑100 (MQ=255)
                                                                                             catcatGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATa  >  1:999299/1‑100 (MQ=255)
                                                                                                |                                                                                                
TTTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATA  >  NC_000913/1382724‑1382916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: