Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,261,554 Δ1 bp 5.8% coding (1087/1116 nt) PP_1104 ← AotO‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,261,5540A. 5.8% 279.1 / 9.0 86coding (1087/1116 nt)PP_1104AotO‑related protein
Reads supporting (aligned to +/- strand):  ref base A (41/40);  new base . (3/2);  total (44/42)
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

CAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTT  >  NC_002947/1261407‑1261697
                                                                                                                                                     |                                                                                                                                               
cAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATa                                                                                                                                               <  2:809711/150‑1 (MQ=255)
 aGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATAc                                                                                                                                              <  1:1705435/150‑1 (MQ=255)
   ccGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTATGCCGCAG‑CGGC‑TACCt                                                                                                                                            <  1:137376/150‑1 (MQ=255)
      cgcTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGtt                                                                                                                                         <  1:864948/150‑1 (MQ=255)
       gcTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCGAACAGGTATGCCGCAG‑CGGC‑TACCTGttg                                                                                                                                        <  1:2084923/150‑2 (MQ=255)
       gcTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTc                                                                                                                                        <  2:1512147/150‑1 (MQ=255)
         ttGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGc                                                                                                                                      >  2:2258640/1‑150 (MQ=255)
                 ccAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCcaca                                                                                                                              >  1:2036501/1‑150 (MQ=255)
                  cAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACAt                                                                                                                             >  1:1492979/1‑150 (MQ=255)
                  cAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACAt                                                                                                                             >  1:1402243/1‑150 (MQ=255)
                   aGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATg                                                                                                                            >  2:557066/1‑150 (MQ=255)
                         tGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTCTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGAt                                                                                                                      <  1:2180276/150‑1 (MQ=255)
                          ggCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATg                                                                                                                     <  2:533575/150‑1 (MQ=255)
                           gCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAGGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGcgacga                                                                                                                    >  2:2082489/1‑150 (MQ=255)
                           gCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTATACAGACCACTTGCATGCTTCAATGGTTTTACCGGATCAATTGACCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATgg                                                                                                                    >  2:409943/1‑149 (MQ=255)
                                     tAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGCGCTTCTACCGTCATCCGGCGAACCGGTATTGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGc                                                                                                          <  1:892141/150‑1 (MQ=255)
                                      agcagcAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCg                                                                                                         <  2:2036501/150‑1 (MQ=255)
                                           cagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGCCCATGCCAGCGGTGGc                                                                                                    <  2:1928841/150‑1 (MQ=255)
                                              gCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCAt                                                                                                 >  1:1837186/1‑150 (MQ=255)
                                              gCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCAt                                                                                                 <  2:1205344/150‑1 (MQ=255)
                                                aTGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATgc                                                                                               >  1:2052886/1‑150 (MQ=255)
                                                aTGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CAACATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATgc                                                                                               >  2:447488/1‑150 (MQ=255)
                                                   aGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGcc                                                                                            >  1:307091/1‑150 (MQ=255)
                                                   aGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGcc                                                                                            <  2:1958051/150‑1 (MQ=255)
                                                    gggCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCg                                                                                           <  2:294719/150‑1 (MQ=255)
                                                    gggAAGGCGGCCCTTCGAAGGCTGTGCTGCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCAAACAGCTC‑GCCGATGACGGC‑TCCATCGCCGTCAGCCACACCCGCTAGGACCATGCCAGCGGTGGCCATGCGCCg                                                                                           >  2:1729689/1‑150 (MQ=255)
                                                           cGGTCCTTCGAATGCTGTTCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGCATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGCTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTa                                                                                    <  2:1179778/150‑1 (MQ=255)
                                                           cGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTa                                                                                    <  2:266615/150‑1 (MQ=255)
                                                            ggTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTAc                                                                                   >  2:1916447/1‑150 (MQ=255)
                                                             gTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATgg                                                                                                                           <  2:1447845/109‑1 (MQ=255)
                                                             gTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATgg                                                                                                                           >  1:1447845/1‑109 (MQ=255)
                                                              tCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTAcgc                                                                                 >  1:181040/1‑150 (MQ=255)
                                                                 ttCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCa                                                                              >  2:995775/1‑150 (MQ=255)
                                                                        gCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagca                                                                       <  1:1993636/150‑1 (MQ=255)
                                                                          tgtgCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACagcagc                                                                     >  2:134166/1‑150 (MQ=255)
                                                                               ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTc                                                                                                                                        >  2:1568680/1‑78 (MQ=255)
                                                                               ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTc                                                                                                                                        <  1:1568680/78‑1 (MQ=255)
                                                                               ttCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCgg                                                                >  2:1142610/1‑150 (MQ=255)
                                                                                 ggtaaaaCCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggct                                                              >  2:1738932/7‑150 (MQ=255)
                                                                                  tACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCTCCGGACGCTACGCGCATACAGCAGCCCCggctg                                                             >  2:476589/1‑150 (MQ=255)
                                                                                  tACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctg                                                             >  1:557179/1‑150 (MQ=255)
                                                                                     aGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctggct                                                          >  1:1874788/1‑150 (MQ=255)
                                                                                       aCCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATATCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGt                                                        <  1:1380445/150‑1 (MQ=255)
                                                                                       aCCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGt                                                        <  2:1219042/150‑1 (MQ=255)
                                                                                        ccACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGtt                                                       <  1:2053229/150‑1 (MQ=255)
                                                                                        ccACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGtt                                                       >  2:923727/1‑150 (MQ=255)
                                                                                           cTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGc                                                                                                                         >  2:1399489/1‑81 (MQ=255)
                                                                                           cTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGc                                                                                                                         <  1:1399489/81‑1 (MQ=255)
                                                                                           cTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagc                                                                        >  1:273887/1‑130 (MQ=255)
                                                                                           cTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagc                                                                        <  2:273887/130‑1 (MQ=255)
                                                                                             tGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCa                                                  <  1:923727/150‑1 (MQ=255)
                                                                                                cGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAggg                                               <  2:1054130/150‑1 (MQ=255)
                                                                                                      cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                         <  1:1338106/150‑1 (MQ=255)
                                                                                                       aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGAc                                        <  2:2052886/150‑1 (MQ=255)
                                                                                                       aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGAc                                        >  2:1959342/1‑150 (MQ=255)
                                                                                                                  cGGATCAATTGGCCAACAGGTA‑GCCGCGG‑CGGCATACCTGTCCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGa                             >  1:1904112/1‑150 (MQ=255)
                                                                                                                   ggATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAg                            <  1:2258640/150‑1 (MQ=255)
                                                                                                                    gATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGt                           >  1:830816/1‑150 (MQ=255)
                                                                                                                    gATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGt                           <  1:620068/150‑1 (MQ=255)
                                                                                                                    gATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGt                           >  2:1499300/1‑150 (MQ=255)
                                                                                                                     aTCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGc                                                                                                              <  2:957515/66‑1 (MQ=255)
                                                                                                                     aTCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGc                                                                                                              >  1:957515/1‑66 (MQ=255)
                                                                                                                        aaTTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGt                       >  1:1975788/1‑150 (MQ=255)
                                                                                                                         aTTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTc                      >  2:1887821/1‑150 (MQ=255)
                                                                                                                         aTTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTAGCCATGCCCCCGACGCTACGCGCATACAGCAGCCCCGTCTGGCTGCTGCGCAGGGTGGTGGCGCGTGCACTGAGTGGGTc                      >  2:2307970/1‑150 (MQ=255)
                                                                                                                            ggCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCTCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGAt                   <  1:476589/150‑1 (MQ=255)
                                                                                                                            ggCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGAt                   <  1:557066/150‑1 (MQ=255)
                                                                                                                            ggCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGAt                   <  2:1492979/150‑1 (MQ=255)
                                                                                                                             gCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAgg                                                >  1:1335102/1‑120 (MQ=255)
                                                                                                                             gCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAgg                                                <  2:1335102/120‑1 (MQ=255)
                                                                                                                               cAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGAc                <  2:181040/150‑1 (MQ=255)
                                                                                                                                aaCAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACt               >  2:1980649/1‑150 (MQ=255)
                                                                                                                                   aggtatgCCGCAG‑CGGC‑TACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCAATTACTTCg            >  1:715065/7‑150 (MQ=255)
                                                                                                                                   aggtatgCCGCAG‑CGGC‑TACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATCCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCAGTTGCGCAGGGGGGTGACTCGGTCACTGAGCGGGGCGATGACTTcc            >  2:731526/7‑149 (MQ=255)
                                                                                                                                   aGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCg            >  2:1522917/1‑150 (MQ=255)
                                                                                                                                     gTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGc          <  2:1975788/150‑1 (MQ=255)
                                                                                                                                     gTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGc          <  2:1874788/150‑1 (MQ=255)
                                                                                                                                     gTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGc          >  1:1672669/1‑150 (MQ=255)
                                                                                                                                       a‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGa        <  1:1522917/150‑1 (MQ=255)
                                                                                                                                         gCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGAt       >  1:733642/1‑150 (MQ=255)
                                                                                                                                          ccGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTgg                         >  2:2102786/1‑131 (MQ=255)
                                                                                                                                          ccGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTgg                         <  1:2102786/131‑1 (MQ=255)
                                                                                                                                           cGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGCGTGGGTCGATGACTTCGGCGATCa     >  2:207332/1‑150 (MQ=255)
                                                                                                                                              tG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGtt  <  2:1837186/150‑1 (MQ=255)
                                                                                                                                                   gCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctgg                                                            >  2:534689/1‑88 (MQ=255)
                                                                                                                                                   gCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctgg                                                            <  1:534689/88‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CAGCCGCGCTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTG‑CGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTT  >  NC_002947/1261407‑1261697

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