Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,938,848 G→A 9.3% P70S (CCG→TCG) ‡ PP_1738 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,938,8480GA 9.3% 208.4 / 9.1 76P70S (CCG→TCG) ‡PP_1738hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/30);  new base A (3/4);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

TCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGATATCCCGCGCT  >  NC_002947/1938701‑1938993
                                                                                                                                                   |                                                                                                                                                 
tCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCgg                                                                                                                                                   <  5:90421/148‑1 (MQ=255)
   ggTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGgcg                                                                                                                                                >  5:321022/1‑148 (MQ=255)
       ttCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGt                                                                                                                                           >  6:95786/1‑149 (MQ=255)
         cATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGtt                                                                                                                                          >  6:169364/1‑148 (MQ=255)
          aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGtt                                                                                                                                          >  3:25269/1‑147 (MQ=255)
          aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCa                                                                                                                                        <  1:313704/149‑1 (MQ=255)
          aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCa                                                                                                                                        <  3:174264/149‑1 (MQ=255)
          aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATAAGCGGCGCTTCCTGATCCTCGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCATGAAGTTGTAGGGCAACGAGAACTCGCGATACGGGCGCAAGTACa                                                                                                                                        >  6:279855/1‑149 (MQ=255)
          aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGAGGATGAAGAGCAGCGGCGCTTCGTGATCGTGGCGGGCCAGGCGGTACTTGAGCTGCTTCACTGGCCAGAAGTGGTCGGGCGACGTGGGGTCGGGCTCGGGgcggcgggg                                                                                                                                          >  6:202405/1‑141 (MQ=255)
             gATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATTATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAgg                                                                                                                                      <  7:364503/148‑1 (MQ=255)
               tgctgctgcTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTTGCGCCCGGAGCGCGAGTTCacgttg                                                                                                                                   <  6:174005/149‑6 (MQ=255)
               tgctgctgcTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCGTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCGTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTg                                                                                                                                   <  8:155912/149‑1 (MQ=255)
                gctgctgctGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTg                                                                                                                                   >  6:253257/1‑148 (MQ=255)
                gctgctgctGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTg                                                                                                                                   >  8:364749/1‑148 (MQ=255)
                         gTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACt                                                                                                                         >  5:64101/1‑149 (MQ=255)
                          tAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGAc                                                                                                                          <  1:370884/147‑1 (MQ=255)
                          tAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACt                                                                                                                         <  6:184728/148‑1 (MQ=255)
                           catcatcgcCGGCACCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTACACTGCCCAGAAGTTTTCGGGCAACGTGAACTCGCGCTCCCGGCGCAAGATCAGGCTGTAGTCGGACt                                                                                                                         >  4:270200/7‑147 (MQ=255)
                            gggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAg                                                                                                                      <  5:317977/149‑1 (MQ=255)
                            gggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAg                                                                                                                      <  5:326995/149‑1 (MQ=255)
                                      gTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtc                                                                                                            >  5:15207/1‑149 (MQ=255)
                                          cGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtcgtcg                                                                                                        <  5:95786/149‑1 (MQ=255)
                                             cGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTg                                                                                                     <  7:364749/149‑1 (MQ=255)
                                                  attaaGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCa                                                                                                 <  5:305927/145‑1 (MQ=255)
                                                        aTCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAgg                                                                                          <  8:365017/149‑1 (MQ=255)
                                                          cAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAggg                                                                                         <  3:34443/148‑1 (MQ=255)
                                                           aGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGgc                                                                                        <  6:378142/148‑1 (MQ=255)
                                                             cggcgCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAATTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCt                                                                                     >  8:175252/1‑149 (MQ=255)
                                                                    tCGTGATCCTGGCGTGCCAGGCGGTACTTGAGCTTCTTCTCTTCGCAGAAGTTGTCGGTCAACTTGATCTCGCGCTCTGTTGGCAAGTTCAGGCTGTAGTGTTTCTGAGTGTTTTCATCGTCGGTGGGCTGCTCAGGGGGCtttttcg                                                                               <  8:8051/148‑6 (MQ=255)
                                                                      gTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCTGGCGCTAGTTCAGGCTGTAGTCGGACTGAGTCATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggc                                                                             <  5:279855/148‑1 (MQ=255)
                                                                      gTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcgtcg                                                                            >  5:16420/1‑149 (MQ=255)
                                                                      gTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcg                                                                            >  1:230412/1‑149 (MQ=255)
                                                                      gTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcg                                                                            <  6:64101/149‑1 (MQ=255)
                                                                       tGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcgg                                                                           >  2:92687/1‑149 (MQ=255)
                                                                        gATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCggcggc                                                                          >  6:164914/1‑149 (MQ=255)
                                                                          tCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgt                                                                        >  6:304016/1‑149 (MQ=255)
                                                                            cTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtgg                                                                      <  2:230412/149‑1 (MQ=255)
                                                                            cTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtg                                                                       <  5:169364/148‑1 (MQ=255)
                                                                             tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGggtggt                                                                     >  8:67494/1‑149 (MQ=255)
                                                                              ggCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggc                                                                             <  1:60068/140‑1 (MQ=255)
                                                                              ggCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggc                                                                             >  2:60068/1‑140 (MQ=255)
                                                                                       ggCGGTACTTGAGCTTCTTCACTGCCCAGAAGTGGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAtggtg                                                           <  5:304016/149‑1 (MQ=255)
                                                                                                 gAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACg                                                 >  7:364623/1‑149 (MQ=255)
                                                                                                 gAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACg                                                 >  6:292246/1‑149 (MQ=255)
                                                                                                        ttcACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTgg                                           >  2:238831/1‑148 (MQ=255)
                                                                                                         tcACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGt                                          >  7:216696/1‑148 (MQ=255)
                                                                                                          cACTGCCCAGAAGTTGTCGGGCAACGTGAACTCCCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGt                                          <  4:25269/147‑1 (MQ=255)
                                                                                                           aCTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGt                                          >  1:45993/1‑146 (MQ=255)
                                                                                                            cTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTaa                                        >  4:376498/1‑147 (MQ=255)
                                                                                                            cTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                                      <  6:15207/149‑1 (MQ=255)
                                                                                                            cTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                                      <  8:216696/149‑1 (MQ=255)
                                                                                                             tGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGaa                                                   <  5:335222/135‑1 (MQ=255)
                                                                                                             tGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGaa                                                   >  6:335222/1‑135 (MQ=255)
                                                                                                                ccAGAAGTTGTCGGGCAACGTGAACTAGCGCTCCGGGCGCAAGTTCAGGCTCTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGCCGCTTGTCCGGCGGCGTTGATGAGATGGAGTCCTCGAACGGATTGCTAAGAggt                                   >  6:271324/1‑147 (MQ=255)
                                                                                                                  aGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTTGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAg                                <  3:270200/149‑1 (MQ=255)
                                                                                                                  aGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAg                                >  7:46875/1‑149 (MQ=255)
                                                                                                                       ttGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAACGTCGCTGGGTAGCTCAGGTCGCTTTTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCAt                            <  5:271324/148‑1 (MQ=255)
                                                                                                                          tCGGGCAACGTGAATTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCt                        <  7:175252/149‑1 (MQ=255)
                                                                                                                          tCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGCTACCCATAGCt                        >  7:187053/1‑149 (MQ=255)
                                                                                                                            gggCCACGTGATCTAGCGCTCCGGGCGCAAGTTCAGGCTGTAGGCGGACTGAGTGTTTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTgg                      <  8:187053/149‑1 (MQ=255)
                                                                                                                             ggCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTgg                      >  8:238409/1‑148 (MQ=255)
                                                                                                                              gCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTgg                      >  3:154261/1‑147 (MQ=255)
                                                                                                                               cAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAg                   >  5:381906/1‑149 (MQ=255)
                                                                                                                               cAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAg                   >  5:387072/1‑149 (MQ=255)
                                                                                                                                aaCGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGc                  <  5:292246/149‑1 (MQ=255)
                                                                                                                                  cctGAACTTGCGCCCGGAGCGCGAGTTCACGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcg                                                                            <  4:238923/87‑1 (MQ=25)
                                                                                                                                  cctGAACTTGCGCCCGGAGCGCGAGTTCACGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcggcg                                                                            >  3:238923/3‑89 (MQ=25)
                                                                                                                                   gTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTc               >  2:270179/1‑149 (MQ=255)
                                                                                                                                   gTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTc               >  4:43714/1‑149 (MQ=255)
                                                                                                                                    tGAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGACAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTg                                            <  7:5429/119‑1 (MQ=255)
                                                                                                                                    tGAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGACAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTg                                            >  8:5429/1‑119 (MQ=255)
                                                                                                                                     gAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGt                                              >  3:178728/1‑116 (MQ=255)
                                                                                                                                     gAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGt                                              <  4:178728/116‑1 (MQ=255)
                                                                                                                                        cTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGCCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGatct          >  7:161803/1‑147 (MQ=255)
                                                                                                                                         tCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGATATc         >  5:351976/1‑149 (MQ=255)
                                                                                                                                                 cGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGATATCCCgcgc   <  5:253257/147‑1 (MQ=255)
                                                                                                                                                 cGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGATATCCCGCGCt  >  7:187801/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTCGATATCCCGCGCT  >  NC_002947/1938701‑1938993

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: ‑