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.

AAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑GGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTG  >  NC_000913/561591‑561781
                                                                                                    |                                                                                            
aaGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCt                                                                                               >  1:2571166/1‑100 (MQ=255)
   gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑ggg                                                                                           >  2:1215812/1‑100 (MQ=255)
   gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑ggg                                                                                           >  2:1527016/1‑100 (MQ=255)
    ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑gggg                                                                                          <  1:1161996/100‑1 (MQ=255)
    ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑gggg                                                                                          >  2:1207126/1‑100 (MQ=255)
    ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑gggg                                                                                          <  2:1678892/100‑1 (MQ=255)
      acaACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑ggggg‑g                                                                                       >  2:2551449/1‑100 (MQ=255)
       caACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑ggg                                                                                           >  2:3059139/1‑96 (MQ=255)
          ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑ggggg‑ggcgg                                                                                   <  2:861557/100‑5 (MQ=255)
            tgcgTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCAACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑ggggg‑ggcggg                                                                                  <  2:2905021/96‑6 (MQ=255)
            tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGt                                                                                  <  2:3149285/100‑1 (MQ=255)
             gAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGta                                                                                 <  1:2152039/100‑1 (MQ=255)
             gAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGta                                                                                 <  1:1768642/100‑1 (MQ=255)
              aGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGtat                                                                                <  2:1649358/100‑1 (MQ=255)
               gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGtata                                                                               >  1:2183476/1‑100 (MQ=255)
               gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGtata                                                                               <  1:3387845/100‑1 (MQ=255)
               gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGtata                                                                               <  1:3225371/100‑1 (MQ=255)
               gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGtata                                                                               <  2:384812/100‑1 (MQ=255)
                 tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATaa                                                                              >  1:2820723/1‑99 (MQ=255)
                 tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAAc                                                                             <  1:1904825/100‑1 (MQ=255)
                  aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACt                                                                            >  1:1860741/1‑100 (MQ=255)
                    cACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTAc                                                                          >  1:1672466/1‑100 (MQ=255)
                      cggaaGCCAGGGTGAGACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACg                                                                         <  1:875813/96‑1 (MQ=255)
                      cGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGa                                                                        >  1:1023246/1‑100 (MQ=255)
                       gTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGAc                                                                       >  2:1260461/1‑100 (MQ=255)
                          aGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACcgc                                                                    >  2:3302139/1‑100 (MQ=255)
                           gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACcgcg                                                                   <  1:252164/100‑1 (MQ=255)
                           gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACcgcg                                                                   <  2:3341156/100‑1 (MQ=255)
                            ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGa                                                                  >  2:316125/1‑100 (MQ=255)
                            ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGa                                                                  <  1:1057415/100‑1 (MQ=255)
                             cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGAt                                                                 <  2:961448/100‑1 (MQ=255)
                             cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGAt                                                                 <  2:3252717/100‑1 (MQ=255)
                             cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGAt                                                                 <  2:917861/100‑1 (MQ=255)
                              aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATc                                                                <  2:1333595/100‑1 (MQ=255)
                              aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATc                                                                <  1:945300/100‑1 (MQ=255)
                              aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATc                                                                <  1:691603/100‑1 (MQ=255)
                               gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCa                                                               <  1:2671738/100‑1 (MQ=255)
                               gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCa                                                               <  2:1679046/100‑1 (MQ=255)
                               gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCa                                                               >  2:2401913/1‑100 (MQ=255)
                                    aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATg                                                          >  2:1783751/1‑100 (MQ=255)
                                     tACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGa                                                         >  2:804545/1‑100 (MQ=255)
                                     tACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGa                                                         >  1:3008578/1‑100 (MQ=255)
                                     tACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGa                                                         >  2:2847848/1‑100 (MQ=255)
                                      aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGAc                                                        >  2:1215249/1‑100 (MQ=255)
                                      aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGAc                                                        >  1:2991814/1‑100 (MQ=255)
                                      aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGAc                                                        <  1:1350574/100‑1 (MQ=255)
                                        cTCGAATAATGGGTACACGGGCGGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGt                                                      >  2:1974725/1‑100 (MQ=255)
                                         tCGAATAATGGGTACACGGGCAGCGCCACGGCAAACGGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGtt                                                     <  2:612537/100‑1 (MQ=255)
                                             ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                 <  1:916771/100‑1 (MQ=255)
                                                 tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                             >  2:2428361/1‑100 (MQ=255)
                                                 tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                             <  2:3008566/100‑1 (MQ=255)
                                                 tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                             >  1:3141337/1‑100 (MQ=255)
                                                  gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                            <  1:2425916/100‑1 (MQ=255)
                                                  gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                            <  1:2425897/100‑1 (MQ=255)
                                                   ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                           <  2:359784/100‑1 (MQ=255)
                                                   ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                           >  1:1426976/1‑100 (MQ=255)
                                                   ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                           <  1:54247/100‑1 (MQ=255)
                                                     tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                         <  1:316126/100‑1 (MQ=255)
                                                     tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                         <  1:2847858/100‑1 (MQ=255)
                                                     tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                         >  2:104770/1‑100 (MQ=255)
                                                         cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                     >  2:3103345/1‑100 (MQ=255)
                                                         cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                     >  1:811070/1‑100 (MQ=255)
                                                         cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                     >  2:3301759/1‑100 (MQ=255)
                                                           ggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGc                                   <  1:2814044/100‑1 (MQ=255)
                                                             naGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                 >  2:1442206/2‑100 (MQ=255)
                                                             cAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                 <  1:1783755/100‑1 (MQ=255)
                                                              aGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgt                                <  2:2303661/100‑1 (MQ=255)
                                                               gcgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtg                               <  1:1207129/100‑1 (MQ=255)
                                                                 gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             <  1:3301772/100‑1 (MQ=255)
                                                                 gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             >  1:3167364/1‑100 (MQ=255)
                                                                 gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             >  2:3038449/1‑100 (MQ=255)
                                                                 gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             >  2:3035529/1‑100 (MQ=255)
                                                                   ttccgatctcTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGt                           <  1:471688/91‑1 (MQ=255)
                                                                    aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                          >  1:2898805/1‑100 (MQ=255)
                                                                    aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                          >  1:320564/1‑100 (MQ=255)
                                                                    aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCCCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTATGGTc                          >  2:1106514/1‑100 (MQ=255)
                                                                     cGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCg                         >  1:1823634/1‑100 (MQ=255)
                                                                      ggCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCgg                        <  1:2428369/100‑1 (MQ=255)
                                                                       gCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                       >  1:2882492/1‑100 (MQ=255)
                                                                       ctaaaCTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                       <  2:1513014/98‑1 (MQ=255)
                                                                        cAAACTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGc                      >  2:947564/1‑100 (MQ=255)
                                                                         aaaCTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCgggag                     >  1:1513017/1‑98 (MQ=255)
                                                                          aaCTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGCTCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGGGGTCgcgc                      >  2:1175123/1‑98 (MQ=255)
                                                                            cTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGtagatcggaa                  >  2:471687/1‑91 (MQ=255)
                                                                            cTGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGa                  <  1:2401920/100‑1 (MQ=255)
                                                                             tGGCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaa                 >  2:939058/1‑100 (MQ=255)
                                                                               gCAGGCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaaa               <  1:2506425/100‑1 (MQ=255)
                                                                               gCAGGCCGCTCACGGTACGCT‑GGGGGCGCTGTTATAACTACGACCGCGATCAACATGACGTTAACTAGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaaa               <  2:2850398/100‑1 (MQ=255)
                                                                                  ggCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg            <  1:3287994/100‑1 (MQ=255)
                                                                                  ggCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg            >  1:428066/1‑100 (MQ=255)
                                                                                  ggCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg            >  2:766301/1‑100 (MQ=255)
                                                                                  ctccGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg            <  1:1354945/98‑1 (MQ=255)
                                                                                   gCCGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGc           >  2:2939669/1‑100 (MQ=255)
                                                                                    ccGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATggag          >  2:1354942/1‑98 (MQ=255)
                                                                                    ccGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa          <  1:3059151/100‑1 (MQ=255)
                                                                                    ccGCTTACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa          <  2:2714744/100‑1 (MQ=255)
                                                                                        ttACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATaa       >  1:2876844/1‑99 (MQ=255)
                                                                                        ttACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAAc      <  2:2876833/100‑1 (MQ=255)
                                                                                         tACGGTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACg     <  1:1215815/100‑1 (MQ=255)
                                                                                            ggTACGCTGGGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTg  >  2:3386947/1‑100 (MQ=255)
                                                                                                    |                                                                                            
AAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑GGGGG‑GCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTG  >  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: ‑