Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,498,216 C→T 7.7% H39H (CAC→CAT) ‡ PP_2193 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,498,2160CT 7.7% 183.3 / 4.5 65H39H (CAC→CAT) ‡PP_2193ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (2/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GACTGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGAT  >  NC_002947/2498069‑2498357
                                                                                                                                                    |                                                                                                                                             
gaCTGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACg                                                                                                                                              >  6:175711/1‑149 (MQ=255)
 aCTGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACg                                                                                                                                              >  7:107343/1‑148 (MQ=255)
     aTCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCCTGCACGCAgg                                                                                                                                          >  5:157945/1‑149 (MQ=255)
     aTCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAggg                                                                                                                                         <  5:175711/149‑1 (MQ=255)
         gCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑‑CCCCTGCGTGCAgggggccg                                                                                                                                    <  2:211181/149‑5 (MQ=255)
         gCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑‑CCCCTGCGTGCAgggggcc                                                                                                                                     >  6:23742/1‑145 (MQ=255)
          cAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCcgc                                                                                                                                    <  2:202582/149‑1 (MQ=255)
          cAACCCATCCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑‑CCCCTGCGTGCAgggggccga                                                                                                                                   <  5:23742/149‑6 (MQ=255)
            aCCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGa                                                                                                                                  >  7:53351/1‑149 (MQ=255)
               cATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACgg                                                                                                                               <  1:57787/149‑1 (MQ=255)
               cATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACAATGCCTTGATTTCTTGTTGTGCCGTGGCCCTGTTGATCGG‑CCCCATGCACGCAGGGGCCGCGACgg                                                                                                                               >  6:328528/1‑149 (MQ=255)
                 ttCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCt                                                                                                                             <  6:85356/149‑1 (MQ=255)
                      aTTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGTCTGc                                                                                                                           <  8:53351/146‑1 (MQ=255)
                        ttACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGaa                                                                                                                       >  2:66172/1‑148 (MQ=255)
                         tACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAAc                                                                                                                      <  7:538/148‑1 (MQ=255)
                            gACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAc                                                                                                                  >  3:196161/1‑149 (MQ=255)
                                 gTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCa                                                                                                               <  4:162614/147‑1 (MQ=255)
                                  tCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAccc                                                                                                                >  2:111180/1‑145 (MQ=255)
                                  tCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGAt                                                                                                            <  1:45304/149‑1 (MQ=255)
                                  tCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGAt                                                                                                            >  5:145618/1‑149 (MQ=255)
                                  tCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGAGGTTCTTTTCACCATGGCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAccc                                                                                                                <  1:111180/145‑1 (MQ=255)
                                      gAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAg                                                                                                        <  2:225792/149‑1 (MQ=255)
                                          cTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCGTGCAg                                                                                                                                           >  7:262738/1‑110 (MQ=255)
                                          cTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCGTGCAg                                                                                                                                           <  8:262738/110‑1 (MQ=255)
                                          cTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAgg                                                                                                       >  3:139408/1‑146 (MQ=255)
                                          cTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAgg                                                                                                       <  4:139408/146‑1 (MQ=255)
                                             aTTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCCGATGCAGGCCCGc                                                                                                  >  5:185051/1‑148 (MQ=255)
                                                 gTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACa                                                                                                                     >  8:129867/1‑125 (MQ=255)
                                                  tGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTAGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGCACAGAACCAGATCCAGCCCCGCCAGGTc                                                                                            >  8:280869/1‑149 (MQ=255)
                                                      tctGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAg                                                                                         >  7:193839/1‑148 (MQ=255)
                                                           tGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGAc                                                                                   <  4:157662/149‑1 (MQ=255)
                                                                aaGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTg                                                                                                                            <  7:229284/103‑1 (MQ=255)
                                                                aaGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTg                                                                                                                            >  8:229284/1‑103 (MQ=255)
                                                                 aGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCa                                                                             >  1:320580/1‑149 (MQ=255)
                                                                  gCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCa                                                                             <  6:157945/149‑1 (MQ=255)
                                                                          aCCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCg                                                                    >  7:22063/1‑149 (MQ=255)
                                                                            cgccgTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGcc                                                                          <  3:64442/141‑1 (MQ=255)
                                                                            cgccgTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGcc                                                                          >  4:64442/1‑141 (MQ=255)
                                                                              ccgTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAc                                                                                                                  >  2:91214/1‑99 (MQ=255)
                                                                              ccgTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAc                                                                                                                  <  1:91214/99‑1 (MQ=255)
                                                                              ccgTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTgg                                                                <  7:280869/149‑1 (MQ=255)
                                                                                    cgGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAgg                                                          >  3:257755/1‑149 (MQ=255)
                                                                                       ttCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGgcgc                                                       <  1:66172/149‑1 (MQ=255)
                                                                                            ttCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGAc                                                   <  7:64351/148‑1 (MQ=255)
                                                                                                     ccTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTc                                                                                            <  6:304897/98‑1 (MQ=255)
                                                                                                     ccTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTc                                                                                            >  5:304897/1‑98 (MQ=255)
                                                                                                     ccTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGcc                                         <  8:22063/149‑1 (MQ=255)
                                                                                                      cTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGcc                                         >  5:267140/1‑148 (MQ=255)
                                                                                                      cTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCg                                        <  6:145618/149‑1 (MQ=255)
                                                                                                           tGTCTTGTAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAg                                   <  5:328528/149‑1 (MQ=255)
                                                                                                                tGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccgg                              >  5:277927/1‑149 (MQ=255)
                                                                                                                 gCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCTGGGGCCGCGACGGCTGCCGaa                                                                                                                       <  1:316511/59‑1 (MQ=255)
                                                                                                                 gCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggc                             >  8:205882/1‑149 (MQ=255)
                                                                                                                     tGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                          >  7:55677/1‑148 (MQ=255)
                                                                                                                     tGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGc                         <  4:257755/149‑1 (MQ=255)
                                                                                                                             cccTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgctg                 <  6:277927/149‑1 (MQ=255)
                                                                                                                             cccTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgctg                 >  5:215150/1‑149 (MQ=255)
                                                                                                                             cccTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgctg                 <  4:196161/149‑1 (MQ=255)
                                                                                                                                 gTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGGCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTaca             >  3:72157/1‑149 (MQ=255)
                                                                                                                                 gTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCACCCGCAGGCCCAGGCGCTGACCACATTCGCCGCCCAGGCCGGCCTGCACCTGCGGGaca             >  8:110797/1‑149 (MQ=255)
                                                                                                                                 gTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGAc                                                              >  7:240540/1‑100 (MQ=255)
                                                                                                                                 gTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGAc                                                              <  8:240540/100‑1 (MQ=255)
                                                                                                                                 gTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTaca             <  8:239197/149‑1 (MQ=255)
                                                                                                                                      tCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTAcacca          >  8:52650/1‑147 (MQ=255)
                                                                                                                                             cccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGAt  <  2:320580/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GACTGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGG‑CCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGAT  >  NC_002947/2498069‑2498357

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