Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,942,955 T→G 9.2% D313A (GAC→GCC) ‡ PP_3478 ← secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,942,9550TG 9.2% 186.8 / 9.6 65D313A (GAC→GCC) ‡PP_3478secretion protein
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base G (2/4);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

ACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTG  >  NC_002947/3942808‑3943103
                                                                                                                                                   |                                                                                                                                                    
aCGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATACGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTc                                                                                                                                                     <  8:140680/149‑1 (MQ=255)
  gCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTcgc                                                                                                                                                   >  3:3156/1‑149 (MQ=255)
    cGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTACATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTagcg                                                                                                                                                  >  1:222408/1‑148 (MQ=255)
     gAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCa                                                                                                                                                >  4:251053/1‑149 (MQ=255)
       aGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCGACACGTCGGACGCGGTGCGCATg                                                                                                                                              <  5:183880/149‑1 (MQ=255)
          tgtTTGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTc                                                                                                                                           <  2:138985/149‑1 (MQ=255)
            ttGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCGACACGTCGGACGCGGTGCGCATggccg                                                                                                                                          <  6:346860/148‑4 (MQ=38)
              gggCAATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATGACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGTTCAgg                                                                                                                                        <  2:222408/148‑1 (MQ=255)
               ggCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAgg                                                                                                                                        <  6:17673/147‑1 (MQ=255)
               ggCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAgg                                                                                                                                        >  5:17673/1‑147 (MQ=255)
                  cATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGtt                                                                                                                                   <  3:213849/149‑1 (MQ=255)
                    ttGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGa                                                                                                                                 <  3:66340/149‑1 (MQ=255)
                      gCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAg                                                                                                                                <  1:86506/148‑1 (MQ=255)
                            cAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCg                                                                                                                         >  1:277225/1‑149 (MQ=255)
                                gTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTc                                                                                                                                                     >  6:161914/1‑117 (MQ=255)
                                gTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTc                                                                                                                                                     <  5:161914/117‑1 (MQ=255)
                                              aGGTGGCCACCGCCATCACCTCCACTTCCAGCACCACCTCGGGGGTGGACTGGGCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCgg                                                                                                       <  4:202309/149‑1 (MQ=255)
                                               ggTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCggg                                                                                                      >  6:342706/1‑149 (MQ=255)
                                                      aCCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAg                                                                                               <  1:153245/149‑1 (MQ=255)
                                                        cTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCa                                                                                             >  8:4273/1‑149 (MQ=255)
                                                          ccATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAgg                                                                                                                                        >  4:297926/1‑104 (MQ=255)
                                                           cATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATgg                                                                                          >  7:205169/1‑149 (MQ=255)
                                                           cATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATg                                                                                           <  4:3156/148‑1 (MQ=255)
                                                           cATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATAGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATgg                                                                                          >  2:55260/1‑149 (MQ=255)
                                                               aCCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTGCTGCGACTGCAGCAGCTTCTCGGCCATGCTCACCGCGTCCGACGTGTCGAGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCtt                                                                                      <  3:249350/149‑1 (MQ=255)
                                                               aCCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCtt                                                                                      >  6:58195/1‑149 (MQ=255)
                                                                ccTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGaa                                                                                                                 <  4:196496/121‑1 (MQ=255)
                                                                ccTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGaa                                                                                                                 >  3:196496/1‑121 (MQ=255)
                                                                ccTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCtt                                                                                      >  3:186890/1‑148 (MQ=255)
                                                                 cTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAg                                                                                                                                <  8:124530/105‑1 (MQ=255)
                                                                 cTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAg                                                                                                                                >  7:124530/1‑105 (MQ=255)
                                                                 cTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAAc                                                                                                                    <  8:271694/117‑1 (MQ=255)
                                                                 cTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAAc                                                                                                                    >  7:271694/1‑117 (MQ=255)
                                                                    cACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTc                                                                                                                          <  3:111343/108‑1 (MQ=255)
                                                                    cACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTc                                                                                                                          >  4:111343/1‑108 (MQ=255)
                                                                     aCTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCAt                                                                                >  3:225803/1‑149 (MQ=255)
                                                                       ttCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGt                                                                              >  2:236535/1‑149 (MQ=255)
                                                                         ccAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGtt                                                                             >  6:307997/1‑148 (MQ=255)
                                                                              accaccTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAAcgc                                                                       >  7:109892/1‑149 (MQ=255)
                                                                                caccTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCgg                                                                     >  3:31325/1‑149 (MQ=255)
                                                                                    tCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGa                                                                                    >  2:64747/1‑130 (MQ=255)
                                                                                    tCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGa                                                                                    <  1:64747/130‑1 (MQ=255)
                                                                                    tCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGtt                                                                 >  7:240119/1‑149 (MQ=255)
                                                                                      ggggTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTTAGCATGGTCTTGACCATGTTCAACGCGGTGTTgg                                                               >  4:181351/1‑149 (MQ=255)
                                                                                      ggggTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTg                                                                <  7:153644/148‑1 (MQ=255)
                                                                                                        cGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAg                                               <  8:109892/147‑1 (MQ=255)
                                                                                                        cGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGt                                             <  8:205169/149‑1 (MQ=255)
                                                                                                           cTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTc                                                           >  5:294038/1‑132 (MQ=255)
                                                                                                           cTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTc                                                           <  6:294038/132‑1 (MQ=255)
                                                                                                                  agcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTgcg                                   <  4:186890/149‑1 (MQ=255)
                                                                                                                       ctcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCAt                                <  3:251053/147‑1 (MQ=255)
                                                                                                                          ggCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAc                             >  1:200775/1‑147 (MQ=255)
                                                                                                                           gCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATg                                                                               <  5:81984/96‑1 (MQ=255)
                                                                                                                           gCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATg                                                                               >  6:81984/1‑96 (MQ=255)
                                                                                                                            ccATGCGCGACACGTCGGACGCGGTGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCa                                                <  7:290822/126‑1 (MQ=38)
                                                                                                                            ccATGCGCGACACGTCGGACGCGGTGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCa                                                >  8:290822/1‑126 (MQ=38)
                                                                                                                             cATGCGCGACACGTCGGACGCGGTGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGt                                                                              >  6:13880/1‑95 (MQ=25)
                                                                                                                             cATGCGCGACACGTCGGACGCGGTGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGt                                                                              <  5:13880/95‑1 (MQ=25)
                                                                                                                              aTGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTc                       <  5:307997/149‑1 (MQ=255)
                                                                                                                                 cgcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTg                    >  6:323287/1‑149 (MQ=255)
                                                                                                                                 cgcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTg                    <  5:58195/149‑1 (MQ=255)
                                                                                                                                 cgcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTg                    >  5:101512/1‑149 (MQ=255)
                                                                                                                                 cgcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTg                    <  4:31325/149‑1 (MQ=255)
                                                                                                                                  gcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTgg                   <  1:236535/149‑1 (MQ=255)
                                                                                                                                              aCGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGtc       >  2:237102/1‑149 (MQ=255)
                                                                                                                                                   tCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTg  >  1:144699/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTG  >  NC_002947/3942808‑3943103

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