Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,226,127 C→A 5.1% G122C (GGC→TGC)  gltI ← glutamate/aspartate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,226,1270CA 5.1% 252.4 / 3.0 79G122C (GGC→TGC) gltIglutamate/aspartate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (36/39);  new base A (2/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

GATGCGCTCGGAGGTGGTACCGGCGGTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGTTGCTTTTTCAGCGCGT  >  NC_002947/1225990‑1226273
                                                                                                                                          |                                                                                                                                                  
gATGCGTTGGTATGTGGTAGCGGCGGTTGTCACCTCGTTCTTGCCGGGCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGc                                                                                                                                          <  6:5470/148‑1 (MQ=255)
  tGCGCTCGGAGGTGGTACCGGCGGTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTgg                                                                                                                                       >  8:332299/1‑149 (MQ=255)
      cTCGGAGGTGGTACCGGCGGTGGTTACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCAAACCTGTTGCTGGCGCt                                                                                                                                   >  5:275065/1‑149 (MQ=255)
      cTCGGAGGTGGTACCGGCGGTGGTTACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCAAACCTGTTGCTGGCGCt                                                                                                                                   >  5:269116/1‑149 (MQ=255)
         ggAGGTGGTACCGGCGGTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGTGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCg                                                                                                                                 <  2:62309/148‑1 (MQ=255)
           aGGTGGTACCGGCGGTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACg                                                                                                                              <  1:403478/149‑1 (MQ=255)
              tggtACCGGCGGTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACgttg                                                                                                                           >  6:167275/1‑149 (MQ=255)
                gtACCGGCGGTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACgttgtt                                                                                                                         <  3:325976/149‑1 (MQ=255)
                gtACCGGCGGTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTAAGCCGGCTGACCGTCCTTGACCTTGGTCTACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACgttgtt                                                                                                                         >  3:642/1‑149 (MQ=255)
                 tACCGGCGGTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGttgttg                                                                                                                        <  5:53682/149‑1 (MQ=255)
                    cggcggTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTtggtg                                                                                                                     <  3:351985/149‑1 (MQ=255)
                    cggcggTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTtggtg                                                                                                                     <  8:360172/149‑1 (MQ=255)
                     ggcggTGGTCACCACGTTCTTGCAGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTtggtgg                                                                                                                    <  2:329514/149‑1 (MQ=255)
                      gcggTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTtggtggt                                                                                                                   >  4:168793/1‑149 (MQ=255)
                            gtCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACgttgtt                                                                                                                         <  1:8474/137‑1 (MQ=255)
                            gtCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACgttgtt                                                                                                                         >  2:8474/1‑137 (MQ=255)
                             tCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCca                                                                                                            <  3:12039/149‑1 (MQ=255)
                              caccacGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCca                                                                                                            >  8:270731/1‑148 (MQ=255)
                               accacGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCcaca                                                                                                          <  2:114014/149‑1 (MQ=255)
                               accacGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTAGCTGGCGCTCGACGTTGTTGGTGGTGGAGCcaca                                                                                                          >  8:331704/1‑149 (MQ=255)
                                 cacGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACt                                                                                                        >  8:41573/1‑149 (MQ=255)
                                  acGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTc                                                                                                       <  1:286732/149‑1 (MQ=255)
                                     ttgttGCCGGCGAGGTCCGGGAAGTCGTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAg                                                                                                    <  8:193415/146‑1 (MQ=255)
                                      tCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCa                                                                                                     <  5:240236/147‑1 (MQ=255)
                                      tCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCa                                                                                                     >  6:240236/1‑147 (MQ=255)
                                       cTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACGTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAgg                                                                                                   <  2:229413/148‑1 (MQ=255)
                                       cTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAgg                                                                                                   <  4:130520/148‑1 (MQ=255)
                                        ttGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTc                                                                                                 >  5:171089/1‑149 (MQ=255)
                                          gCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAg                                                                                                    >  2:55972/1‑144 (MQ=255)
                                          gCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAg                                                                                                    <  1:55972/144‑1 (MQ=255)
                                               ccAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTtggtg                                                                                                                     <  7:322424/122‑1 (MQ=255)
                                               ccAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTtggtg                                                                                                                     >  8:322424/1‑122 (MQ=255)
                                               ccAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCAAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTg                                                                                          <  6:275065/149‑1 (MQ=255)
                                               ccAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCAAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTg                                                                                          <  6:269116/149‑1 (MQ=255)
                                                 aGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGcc                                                                                        >  7:191985/1‑149 (MQ=255)
                                                 aGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGcc                                                                                        <  5:167275/149‑1 (MQ=255)
                                                        ggAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGc                                                                                 <  6:171089/149‑1 (MQ=255)
                                                          aaGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCAc                                                                               <  7:332299/149‑1 (MQ=255)
                                                              cTTTGTATTCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTg                                                                                  <  8:61900/142‑1 (MQ=255)
                                                              cTTTGTATCCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTg                                                                                  >  7:61900/1‑142 (MQ=255)
                                                               tttGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGa                                                                                                                >  1:186093/1‑111 (MQ=255)
                                                               tttGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGa                                                                                                                <  2:186093/111‑1 (MQ=255)
                                                               tttGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCg                                                                          >  4:349447/1‑149 (MQ=255)
                                                               tttGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCg                                                                          >  2:23217/1‑149 (MQ=255)
                                                                 tGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGc                                                                           <  1:157057/146‑1 (MQ=255)
                                                                 tGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGc                                                                           <  7:270731/146‑1 (MQ=255)
                                                                    aTGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTc                                                                                                 >  7:354737/1‑121 (MQ=255)
                                                                    aTGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTc                                                                                                 <  8:354737/121‑1 (MQ=255)
                                                                       ccGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACgg                                                                  <  3:349447/149‑1 (MQ=255)
                                                                        cGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACggg                                                                 >  8:362432/1‑149 (MQ=255)
                                                                         ggCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGt                                                                <  8:191985/149‑1 (MQ=255)
                                                                                cGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAg                                                         >  2:237889/1‑149 (MQ=255)
                                                                                  tCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGTAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAgg                                                        >  1:95133/1‑149 (MQ=255)
                                                                                  ctgttGACCAAGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGt                                                       >  5:370626/4‑149 (MQ=255)
                                                                                    cTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGa                                                     >  5:379805/1‑149 (MQ=255)
                                                                                     ttGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGAc                                                    >  1:82066/1‑149 (MQ=255)
                                                                                         ccTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAgg                                                >  6:28625/1‑149 (MQ=255)
                                                                                          cTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGGGAGGTGACCAGGt                                               <  4:642/149‑1 (MQ=255)
                                                                                          cTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGt                                               <  7:362432/149‑1 (MQ=255)
                                                                                              gTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGt                                            >  1:8800/1‑148 (MQ=255)
                                                                                               tCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTaa                                          >  7:324940/1‑149 (MQ=255)
                                                                                                   cAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGa                                      <  3:168793/149‑1 (MQ=255)
                                                                                                          gTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACTGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGAt                               <  3:326430/149‑1 (MQ=255)
                                                                                                          gTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGAt                               <  1:237889/149‑1 (MQ=255)
                                                                                                           tACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATg                              >  5:225263/1‑149 (MQ=255)
                                                                                                            aCCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTCTTGGTGGTGGAGCCACACTCCAGGTCGACCGTGCAGTTCTGCACCAGCGGAATACGGGTTAGCGAGGTGACCAGGTTGTATCGGACCTTGATGt                             >  6:71863/1‑149 (MQ=255)
                                                                                                                  ttCGTAGATTCCGACCTAGAATCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGTTTg                                             <  3:38377/127‑1 (MQ=255)
                                                                                                                    cGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGcc                      >  1:102668/1‑148 (MQ=255)
                                                                                                                    cGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGc                       <  7:370059/147‑1 (MQ=255)
                                                                                                                     gAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGcc                      >  2:384963/1‑147 (MQ=255)
                                                                                                                      aaGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAg                   <  7:41573/149‑1 (MQ=255)
                                                                                                                        gATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGtt                 <  5:185465/149‑1 (MQ=255)
                                                                                                                        gATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGtt                 >  4:256401/1‑149 (MQ=255)
                                                                                                                           gCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGTTGCt              >  4:187195/1‑149 (MQ=255)
                                                                                                                           gCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGAGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGTTGCt              >  1:299337/1‑149 (MQ=255)
                                                                                                                             cgaccgaGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGTTGCttt            <  1:23217/149‑1 (MQ=255)
                                                                                                                              gaccgaGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGTTGCtttt           >  8:117001/1‑149 (MQ=255)
                                                                                                                                 cgaGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGTTGCTTTTTCa        <  3:187195/149‑1 (MQ=255)
                                                                                                                                        gCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACAAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGTTGCTTTTTCAGCGCGt  >  3:68724/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
GATGCGCTCGGAGGTGGTACCGGCGGTGGTCACCACGTTCTTGCCGGCCAGGTCCGGGAAGTCTTTGTATGCCGGCTGACCGTCCTTGACCTTGGTCAACAGGCGGGTACCGACTTCGAAGATGCCGACCGAGAAGCCAACCTGTTGCTGGCGCTCGACGTTGTTGGTGGTGGAGCCACACTCCAGGTCGACGGTGCCGTTCTGCACCAGCGGAATACGGGTTTGCGAGGTGACCAGGTTGTAACGGACCTTGATGTCCGTGCCCAGTTGCTTTTTCAGCGCGT  >  NC_002947/1225990‑1226273

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: