Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,701,985 C→A 6.8% L139M (CTG→ATG) ‡ PP_3267 → Clp protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,701,9850CA 6.8% 230.9 / 5.2 73L139M (CTG→ATG) ‡PP_3267Clp protease
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base A (2/3);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

GGCAATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCG  >  NC_002947/3701838‑3702132
                                                                                                                                                      |                                                                                                                                                    
ggCAATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGcc                                                                                                                                                      <  8:49037/148‑1 (MQ=255)
          tctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGATCTGCGCATACAGGCGCAg                                                                                                                                           >  4:296091/1‑149 (MQ=255)
          tctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                              >  2:71419/1‑146 (MQ=255)
           ctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                              >  7:364705/1‑145 (MQ=255)
              ccTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCATACAGGCGCAGGTc                                                                                                                                        <  6:249704/148‑1 (MQ=255)
               cTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATAGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATc                                                                                                                                        <  1:70103/148‑1 (MQ=255)
               cTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGATCTGCGCATACAGGCGCAGGTCg                                                                                                                                       <  4:64517/148‑1 (MQ=255)
               cTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGAGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGt                                                                                                                                      >  6:39265/1‑149 (MQ=255)
                  gCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGATCTGCGCATACAGGCGCAGGTCGtaa                                                                                                                                    <  3:296091/148‑3 (MQ=255)
                    cGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCg                                                                                                                                 >  6:275979/1‑148 (MQ=255)
                             aCGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACg                                                                                                                        >  8:16274/1‑148 (MQ=255)
                               gCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGa                                                                                                                       >  8:45134/1‑147 (MQ=255)
                                cACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAg                                                                                                                    >  1:8923/1‑149 (MQ=255)
                                 aCGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcgc                                                                                                                                             >  2:17785/1‑124 (MQ=255)
                                 aCGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcgc                                                                                                                                             <  1:17785/124‑1 (MQ=255)
                                  cGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGgc                                                                                                                  <  8:364705/149‑1 (MQ=255)
                                          aGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGa                                                                                                          <  4:236064/149‑1 (MQ=255)
                                                 gCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCATGAAAGGCGCGGCCGAGGTTTTg                                                                                                   <  3:84104/149‑1 (MQ=255)
                                                 gCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                                   >  6:167487/1‑149 (MQ=255)
                                                  cAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTg                                                                                                   >  8:57104/1‑148 (MQ=255)
                                                   aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                 <  7:45134/149‑1 (MQ=255)
                                                   aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                 >  6:147270/1‑149 (MQ=255)
                                                   aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                 <  7:145758/149‑1 (MQ=255)
                                                     gTTTCTTTTCCGCCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATg                                                                                               >  4:265875/1‑149 (MQ=255)
                                                     gTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTGACCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                 >  1:112083/1‑149 (MQ=255)
                                                     gTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATg                                                                                               >  7:143113/1‑149 (MQ=255)
                                                              ccACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGt                                                                                       >  8:329601/1‑148 (MQ=255)
                                                                aCCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTAt                                                                                    <  7:57104/149‑1 (MQ=255)
                                                                    gTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTaa                                                                                >  2:116935/1‑149 (MQ=255)
                                                                     tttCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAAt                                                                               <  1:302035/149‑1 (MQ=255)
                                                                     tttCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCCAAt                                                                               <  5:39265/149‑1 (MQ=255)
                                                                       tCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATg                                                                              <  7:16274/148‑1 (MQ=255)
                                                                        cATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGTGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATgg                                                                             <  1:268574/148‑1 (MQ=255)
                                                                        cATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATgg                                                                             <  8:255044/148‑1 (MQ=255)
                                                                         aTTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGcc                                                                           >  7:350907/1‑149 (MQ=255)
                                                                         aTTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                            <  2:112083/148‑1 (MQ=255)
                                                                         aTTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                            <  5:275979/148‑1 (MQ=255)
                                                                         aTTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGc                                                                            <  7:329601/148‑1 (MQ=255)
                                                                             gACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATc                                                                       >  7:16705/1‑149 (MQ=255)
                                                                                cTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATcc                                                                      >  1:146579/1‑147 (MQ=255)
                                                                                cTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCaa                                                                    <  5:259278/149‑1 (MQ=255)
                                                                                     aTGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCTCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATa                                                               <  7:177696/149‑1 (MQ=255)
                                                                                         ttCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTgg                                                           <  1:356883/149‑1 (MQ=255)
                                                                                             gTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGcgcc                                                       >  7:311261/1‑149 (MQ=255)
                                                                                                  ccATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATgg                                                                             <  6:46003/122‑1 (MQ=255)
                                                                                                  ccATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATgg                                                                             >  5:46003/1‑122 (MQ=255)
                                                                                                    aTTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGcgcctcccgcg                                                >  6:355705/1‑142 (MQ=255)
                                                                                                     ttCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAgg                                               <  7:277752/149‑1 (MQ=255)
                                                                                                          ggATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTg                                          >  7:93480/1‑149 (MQ=255)
                                                                                                          ggATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTg                                          <  8:143113/149‑1 (MQ=255)
                                                                                                                   cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATgg                                                                             >  6:105245/1‑107 (MQ=255)
                                                                                                                   cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATgg                                                                             <  5:105245/107‑1 (MQ=255)
                                                                                                                   cGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACg                                   <  3:270137/149‑1 (MQ=255)
                                                                                                                    gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                                  <  1:71419/149‑1 (MQ=255)
                                                                                                                    gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                                  <  2:379220/149‑1 (MQ=255)
                                                                                                                    gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                                  >  4:160629/1‑149 (MQ=255)
                                                                                                                     aGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                                  >  3:248296/1‑148 (MQ=255)
                                                                                                                            gATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATcc                          <  2:8923/149‑1 (MQ=255)
                                                                                                                            gATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATc                           >  8:11582/1‑148 (MQ=255)
                                                                                                                              tGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCAGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCa                         >  1:74346/1‑149 (MQ=255)
                                                                                                                                 aGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATg                       <  7:11582/148‑1 (MQ=255)
                                                                                                                                   ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGaaa                    >  6:157739/1‑149 (MQ=255)
                                                                                                                                    cTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAAt                   >  7:141382/1‑149 (MQ=255)
                                                                                                                                     tGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATc                  >  3:272697/1‑149 (MQ=255)
                                                                                                                                        aTTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac               <  8:141382/149‑1 (MQ=255)
                                                                                                                                         ttACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac               >  1:29540/1‑148 (MQ=255)
                                                                                                                                          tACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcga                >  5:56358/1‑146 (MQ=255)
                                                                                                                                          tACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcga                <  6:56358/146‑1 (MQ=255)
                                                                                                                                          caCGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac               >  6:25048/2‑147 (MQ=255)
                                                                                                                                             gACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGAt          <  8:350907/149‑1 (MQ=255)
                                                                                                                                              aCCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATg         >  6:288482/1‑149 (MQ=255)
                                                                                                                                                cTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGACCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGcc       >  1:81783/1‑149 (MQ=255)
                                                                                                                                                 tGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGga                                >  7:316318/1‑123 (MQ=255)
                                                                                                                                                 tGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGga                                <  8:316318/123‑1 (MQ=255)
                                                                                                                                                     ccTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATcg  >  2:312907/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GGCAATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGT‑‑CCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCACGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCG  >  NC_002947/3701838‑3702132

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: