Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 561,696 (G)6→7 coding (2000/2604 nt) sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,6901.G99.0% 378.6 / ‑3.0 102coding (1994/2604 nt)sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (49/52);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

AAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTG  >  NC_000913/561591‑561781
                                                                                                    |                                                                                            
aaGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCt                                                                                               >  1:2571166/1‑100 (MQ=255)
   gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggg                                                                                           >  2:1215812/1‑100 (MQ=255)
   gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggg                                                                                           >  2:1527016/1‑100 (MQ=255)
    ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg                                                                                          <  1:1161996/100‑1 (MQ=255)
    ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg                                                                                          >  2:1207126/1‑100 (MQ=255)
    ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg                                                                                          <  2:1678892/100‑1 (MQ=255)
      acaACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggg                                                                                       >  2:2551449/1‑100 (MQ=255)
       caACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggg                                                                                           >  2:3059139/1‑96 (MQ=255)
          ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcgg                                                                                   <  2:861557/100‑5 (MQ=255)
            tgcgTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCAACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcggg                                                                                  <  2:2905021/96‑6 (MQ=255)
            tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGt                                                                                  <  2:3149285/100‑1 (MQ=255)
             gAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGta                                                                                 <  1:2152039/100‑1 (MQ=255)
             gAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGta                                                                                 <  1:1768642/100‑1 (MQ=255)
              aGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtat                                                                                <  2:1649358/100‑1 (MQ=255)
               gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtata                                                                               >  1:2183476/1‑100 (MQ=255)
               gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtata                                                                               <  1:3387845/100‑1 (MQ=255)
               gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtata                                                                               <  1:3225371/100‑1 (MQ=255)
               gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtata                                                                               <  2:384812/100‑1 (MQ=255)
                 tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATaa                                                                              >  1:2820723/1‑99 (MQ=255)
                 tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAAc                                                                             <  1:1904825/100‑1 (MQ=255)
                  aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACt                                                                            >  1:1860741/1‑100 (MQ=255)
                    cACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTAc                                                                          >  1:1672466/1‑100 (MQ=255)
                      cggaaGCCAGGGTGAGACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACg                                                                         <  1:875813/96‑1 (MQ=255)
                      cGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGa                                                                        >  1:1023246/1‑100 (MQ=255)
                       gTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGAc                                                                       >  2:1260461/1‑100 (MQ=255)
                          aGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcgc                                                                    >  2:3302139/1‑100 (MQ=255)
                           gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcgcg                                                                   <  1:252164/100‑1 (MQ=255)
                           gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcgcg                                                                   <  2:3341156/100‑1 (MQ=255)
                            ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGa                                                                  >  2:316125/1‑100 (MQ=255)
                            ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGa                                                                  <  1:1057415/100‑1 (MQ=255)
                             cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGAt                                                                 <  2:961448/100‑1 (MQ=255)
                             cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGAt                                                                 <  2:3252717/100‑1 (MQ=255)
                             cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGAt                                                                 <  2:917861/100‑1 (MQ=255)
                              aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATc                                                                <  2:1333595/100‑1 (MQ=255)
                              aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATc                                                                <  1:945300/100‑1 (MQ=255)
                              aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATc                                                                <  1:691603/100‑1 (MQ=255)
                               gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCa                                                               <  1:2671738/100‑1 (MQ=255)
                               gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCa                                                               <  2:1679046/100‑1 (MQ=255)
                               gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCa                                                               >  2:2401913/1‑100 (MQ=255)
                                    aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATg                                                          >  2:1783751/1‑100 (MQ=255)
                                     tACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGa                                                         >  2:804545/1‑100 (MQ=255)
                                     tACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGa                                                         >  1:3008578/1‑100 (MQ=255)
                                     tACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGa                                                         >  2:2847848/1‑100 (MQ=255)
                                      aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                        >  2:1215249/1‑100 (MQ=255)
                                      aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                        >  1:2991814/1‑100 (MQ=255)
                                      aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                        <  1:1350574/100‑1 (MQ=255)
                                        cTCGAATAATGGGTACACGGGCGGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                      >  2:1974725/1‑100 (MQ=255)
                                         tCGAATAATGGGTACACGGGCAGCGCCACGGCAAACGGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGtt                                                     <  2:612537/100‑1 (MQ=255)
                                             ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                 <  1:916771/100‑1 (MQ=255)
                                                 tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                             >  2:2428361/1‑100 (MQ=255)
                                                 tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                             <  2:3008566/100‑1 (MQ=255)
                                                 tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                             >  1:3141337/1‑100 (MQ=255)
                                                  gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                            <  1:2425916/100‑1 (MQ=255)
                                                  gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                            <  1:2425897/100‑1 (MQ=255)
                                                   ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                           <  2:359784/100‑1 (MQ=255)
                                                   ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                           >  1:1426976/1‑100 (MQ=255)
                                                   ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                           <  1:54247/100‑1 (MQ=255)
                                                     tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                         <  1:316126/100‑1 (MQ=255)
                                                     tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                         <  1:2847858/100‑1 (MQ=255)
                                                     tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                         >  2:104770/1‑100 (MQ=255)
                                                         cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                     >  2:3103345/1‑100 (MQ=255)
                                                         cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                     >  1:811070/1‑100 (MQ=255)
                                                         cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                     >  2:3301759/1‑100 (MQ=255)
                                                           ggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGc                                   <  1:2814044/100‑1 (MQ=255)
                                                             naGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                 >  2:1442206/2‑100 (MQ=255)
                                                             cAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                 <  1:1783755/100‑1 (MQ=255)
                                                              aGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgt                                <  2:2303661/100‑1 (MQ=255)
                                                               gcgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtg                               <  1:1207129/100‑1 (MQ=255)
                                                                 gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             <  1:3301772/100‑1 (MQ=255)
                                                                 gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             >  1:3167364/1‑100 (MQ=255)
                                                                 gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             >  2:3038449/1‑100 (MQ=255)
                                                                 gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             >  2:3035529/1‑100 (MQ=255)
                                                                   ttccgatctcTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGt                           <  1:471688/91‑1 (MQ=255)
                                                                    aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                          >  1:2898805/1‑100 (MQ=255)
                                                                    aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                          >  1:320564/1‑100 (MQ=255)
                                                                    aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCCCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTATGGTc                          >  2:1106514/1‑100 (MQ=255)
                                                                     cGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCg                         >  1:1823634/1‑100 (MQ=255)
                                                                      ggCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCgg                        <  1:2428369/100‑1 (MQ=255)
                                                                       gCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                       >  1:2882492/1‑100 (MQ=255)
                                                                       ctaaaCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                       <  2:1513014/98‑1 (MQ=255)
                                                                        cAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGc                      >  2:947564/1‑100 (MQ=255)
                                                                         aaaCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCgggag                     >  1:1513017/1‑98 (MQ=255)
                                                                          aaCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGCTCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGGGGTCgcgc                      >  2:1175123/1‑98 (MQ=255)
                                                                            cTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGtagatcggaa                  >  2:471687/1‑91 (MQ=255)
                                                                            cTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGa                  <  1:2401920/100‑1 (MQ=255)
                                                                             tGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaa                 >  2:939058/1‑100 (MQ=255)
                                                                               gCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaaa               <  1:2506425/100‑1 (MQ=255)
                                                                               gCAGGCCGCTCACGGTACGCTGGGGGCGCTGTTATAACTACGACCGCGATCAACATGACGTTAACTAGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaaa               <  2:2850398/100‑1 (MQ=255)
                                                                                  ggCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg            <  1:3287994/100‑1 (MQ=255)
                                                                                  ggCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg            >  1:428066/1‑100 (MQ=255)
                                                                                  ggCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg            >  2:766301/1‑100 (MQ=255)
                                                                                  ctccGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg            <  1:1354945/98‑1 (MQ=255)
                                                                                   gCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGc           >  2:2939669/1‑100 (MQ=255)
                                                                                    ccGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATggag          >  2:1354942/1‑98 (MQ=255)
                                                                                    ccGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa          <  1:3059151/100‑1 (MQ=255)
                                                                                    ccGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa          <  2:2714744/100‑1 (MQ=255)
                                                                                        ttACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATaa       >  1:2876844/1‑99 (MQ=255)
                                                                                        ttACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAAc      <  2:2876833/100‑1 (MQ=255)
                                                                                         tACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACg     <  1:1215815/100‑1 (MQ=255)
                                                                                            ggTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTg  >  2:3386947/1‑100 (MQ=255)
                                                                                                    |                                                                                            
AAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTG  >  NC_000913/561591‑561781

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