Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,757,734 T→G 5.5% W264G (TGG→GGG)  PP_1563 → terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,757,7340TG 5.5% 279.3 / 4.4 87W264G (TGG→GGG) PP_1563terminase large subunit
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base G (3/2);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAGCAGG  >  NC_002947/1757595‑1757881
                                                                                                                                            |                                                                                                                                                   
ggAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGGCCCTGCCACCAGTTGATCCGCGATTCCGCGCGATTTCTGGAGGGTGTGATCGACC‑GGCCGGACCTATggccggcc                                                                                                                                             <  5:211074/148‑1 (MQ=255)
ggAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATggccggcc                                                                                                                                             <  3:35577/148‑1 (MQ=255)
  agaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGCCATAGgtccggccggt                                                                                                                                          <  7:286353/149‑4 (MQ=255)
  agaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCctc                                                                                                                                          <  6:164268/149‑1 (MQ=255)
   gaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCctct                                                                                                                                         <  5:262154/149‑1 (MQ=255)
   gaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGAGAATGCTGGAGGGTGTGAGCTACC‑GGCCGGACCTATGGCCGGcccga                                                                                                                                          >  5:283404/1‑146 (MQ=255)
               ggCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGAc                                                                                                                              >  2:124883/1‑148 (MQ=255)
                  ccGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAgg                                                                                                                          <  6:267148/149‑1 (MQ=255)
                  ccGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAgg                                                                                                                          >  4:131797/1‑149 (MQ=255)
                  ccGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAgg                                                                                                                          >  8:33775/1‑149 (MQ=255)
                        gCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCAGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGcgac                                                                                                                     >  7:35656/1‑148 (MQ=255)
                          aaCCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGAt                                                                                                                  >  1:157700/1‑149 (MQ=255)
                          aaCCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGAt                                                                                                                  <  7:192159/149‑1 (MQ=255)
                           aCCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGAtt                                                                                                                 >  4:99938/1‑149 (MQ=255)
                           aCCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGAtt                                                                                                                 >  8:243596/1‑149 (MQ=255)
                            cccATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGCCATAGgtccggccgg                                                                                                                                           >  5:6501/1‑120 (MQ=255)
                            cccATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGCCATAGgtccggccgg                                                                                                                                           <  6:6501/122‑3 (MQ=255)
                             ccagcttGCTGGTTCTCACCTCGTCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGGGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCTACCAGGGCGACGATTgg                                                                                                               <  6:283404/143‑1 (MQ=255)
                              catcatGCTGTTTATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACAAGTTGATCCGCGATTCCGAGCGAATGCTTTAGGGTGTGATCGACC‑GGCCGTACCTATTGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGa                                                                                                              <  7:207925/149‑1 (MQ=255)
                               atcgtGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGCCCAGGGCGACGATTGGAc                                                                                                             >  2:364976/1‑149 (MQ=255)
                               atcatGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGAc                                                                                                             >  2:266831/1‑149 (MQ=255)
                               atcatGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGAc                                                                                                             >  5:334527/1‑149 (MQ=255)
                                tcatGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGAcc                                                                                                            >  8:201014/1‑149 (MQ=255)
                                   tGCTGGTGATCTGCAGGGCGGTCTCAAGTATTGGCGGGCCATGCCACCTGGTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGAGC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGc                                                                                                         <  8:160689/149‑1 (MQ=255)
                                      tGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGa                                                                                                       >  7:116806/1‑148 (MQ=255)
                                      tGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACACGCGa                                                                                                       >  8:140057/1‑148 (MQ=255)
                                        gTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCa                                                                                                           <  1:174360/142‑1 (MQ=255)
                                        gTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCa                                                                                                           >  2:174360/1‑142 (MQ=255)
                                        gTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGaa                                                                                                      <  5:280850/147‑1 (MQ=255)
                                         tGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATc                                                                                                   <  3:131797/149‑1 (MQ=255)
                                          gATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCg                                                                                                  >  2:88786/1‑149 (MQ=255)
                                            tCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCa                                                                                                <  8:116806/149‑1 (MQ=255)
                                              accacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGc                                                                                                         >  1:9801/1‑138 (MQ=255)
                                              accacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGc                                                                                                         <  2:9801/138‑1 (MQ=255)
                                               ccacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGccc                                                                                                                                            <  8:37826/102‑1 (MQ=255)
                                               ccacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGccc                                                                                                                                            >  7:37826/1‑102 (MQ=255)
                                               ccacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTg                                                                                             >  3:216239/1‑149 (MQ=255)
                                                   ggcggGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCa                                                                                         <  7:33775/149‑1 (MQ=255)
                                                     cggGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAg                                                                                       <  1:124883/149‑1 (MQ=255)
                                                        gCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCa                                                                                                           <  5:190958/126‑1 (MQ=255)
                                                        gCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCa                                                                                                           >  6:190958/1‑126 (MQ=255)
                                                            aaGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGaa                                                                                  <  4:362453/147‑1 (MQ=255)
                                                               tATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGAGCAGCGAAATCGCACTGCGCAAGGCGAATCCa                                                                              >  5:285714/1‑148 (MQ=255)
                                                               tATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCa                                                                              >  5:331036/1‑148 (MQ=255)
                                                                  tGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAAc                                                                           <  2:157700/148‑1 (MQ=255)
                                                                    gcggcCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACttt                                                                        >  3:166375/1‑149 (MQ=255)
                                                                    gcggcCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACttt                                                                        >  6:138360/1‑149 (MQ=255)
                                                                     cggcCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                                       <  7:201014/149‑1 (MQ=255)
                                                                      ggcCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTgg                                                                      <  7:243596/149‑1 (MQ=255)
                                                                         cccTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGa                                                                    >  8:178836/1‑148 (MQ=255)
                                                                         cccTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATAGCCC‑GGCCGGACCTAGGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGAGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGa                                                                    >  5:217887/1‑148 (MQ=255)
                                                                           cTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGata                                                                  >  1:318272/1‑148 (MQ=255)
                                                                            tGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATc                                                                >  7:46044/1‑149 (MQ=255)
                                                                              ccaccaGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATcc                                                               <  1:88786/148‑1 (MQ=255)
                                                                                accaGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTc                                                            <  6:334527/149‑1 (MQ=255)
                                                                                     ttGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGa                                                       <  1:266831/149‑1 (MQ=255)
                                                                                        aTCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAgg                                                     <  4:216239/148‑1 (MQ=255)
                                                                                         tCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACCTTGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAgg                                                     <  1:125856/148‑1 (MQ=255)
                                                                                         tCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGAc                                                   <  7:178836/149‑1 (MQ=255)
                                                                                          ccGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACt                                                  >  5:32814/1‑149 (MQ=255)
                                                                                          ccGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACt                                                  >  4:326702/1‑149 (MQ=255)
                                                                                             cgATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTcc                                               <  6:285714/149‑1 (MQ=255)
                                                                                              gATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCt                                              <  4:166375/149‑1 (MQ=255)
                                                                                                ttCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCTAACTTTGGGATATCCGTCGCCGAGGACTTCCTGc                                            <  1:364976/149‑1 (MQ=255)
                                                                                                 tCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATTGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCTCACTGCGCAAGTCGAATCCAAACTTTTGGATATCCGTCGCCGAGGACTTCCTGc                                            <  8:35656/148‑1 (MQ=255)
                                                                                                   cGAGCGTATGCTGGAGGGTCTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGGCGAGGACTTCCTGCtt                                          <  4:322638/148‑1 (MQ=255)
                                                                                                      gCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGTCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTc                                                <  3:27984/139‑1 (MQ=255)
                                                                                                      gCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGACCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTc                                                >  4:27984/1‑139 (MQ=255)
                                                                                                             cTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGc                                <  5:138360/148‑1 (MQ=255)
                                                                                                                gAGGGTGTGATCGAAC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTg                             <  7:140057/148‑1 (MQ=255)
                                                                                                                       tgATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgatgcgatg                     <  8:46044/149‑1 (MQ=255)
                                                                                                                       tgATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgatgcgat                      <  3:99938/148‑1 (MQ=255)
                                                                                                                         aTCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCGAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgat                           <  1:365570/141‑1 (MQ=255)
                                                                                                                         aTCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgat                           >  2:365570/1‑141 (MQ=255)
                                                                                                                          tCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAg                  >  7:84709/1‑149 (MQ=255)
                                                                                                                          tCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAg                  >  2:156313/1‑149 (MQ=255)
                                                                                                                            gACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTc                >  7:354072/1‑149 (MQ=255)
                                                                                                                             aCC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCg               <  3:21905/149‑1 (MQ=255)
                                                                                                                              cc‑ggccggACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCgg              >  5:313956/1‑149 (MQ=255)
                                                                                                                               c‑ggccggACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCTCGGCAGCGTGATGCGATGCAGTCGGc             <  1:353882/149‑1 (MQ=255)
                                                                                                                                gggccggccATAGGTCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCCGc             >  8:231485/2‑149 (MQ=255)
                                                                                                                                  gccggACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCa                                                                                                <  3:305420/64‑1 (MQ=255)
                                                                                                                                  gccggACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCa                                                                                                >  4:305420/1‑64 (MQ=255)
                                                                                                                                      gACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCaa       >  2:217223/1‑149 (MQ=255)
                                                                                                                                         cTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAg      >  2:53885/1‑147 (MQ=255)
                                                                                                                                         cTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAGCa    >  4:230871/1‑149 (MQ=255)
                                                                                                                                           aTGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAGCAgg  >  2:34505/1‑149 (MQ=255)
                                                                                                                                            tGGCGGGCGCTCTACACTATCGACCAGGGCGACGATTGGACGAGCGAAATCGCACTGCGCAAGGCGAAACCAAAGTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAGCAg   <  6:217887/147‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCGAATGCTGGAGGGTGTGATCGACC‑GGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGACCAAGCAGG  >  NC_002947/1757595‑1757881

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