Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,554,034 C→T 5.2% L312L (CTG→CTA)  PP_4041 ← glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,554,0340CT 5.2% 228.2 / 2.0 76L312L (CTG→CTA) PP_4041glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base C (40/32);  new base T (2/2);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

TGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGCCAGCAGGGTCATGGTCGCGTCGCG  >  NC_002947/4553885‑4554180
                                                                                                                                                     |                                                                                                                                                  
tGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCAACATGGCGTTCGCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGc                                                                                                                                                    >  1:737812/1‑150 (MQ=255)
   cGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAgg                                                                                                                                                 >  1:707636/1‑150 (MQ=255)
    ggggCAACCAGCCCTTGATCGCATGAGTCAAAGCATCGGCCAGTTCGACATAGAAGTCAAGTTGCCTCTGCGACGAGGCGGAGTTACCCACCCGCACAGGCTAAGAATGCTCATAACCGGACAGCCACAGCAAGGTGTACTCCTGCAGGt                                                                                                                                                >  2:649865/1‑150 (MQ=255)
          aCCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCt                                                                                                                                          >  2:459702/1‑150 (MQ=255)
            ccGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCa                                                                                                                                                   <  2:912426/139‑1 (MQ=255)
            ccGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCa                                                                                                                                                   >  1:912426/1‑139 (MQ=255)
                        gCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGCCCACCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCa                                                                                                                            <  1:682315/150‑1 (MQ=255)
                        accTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCa                                                                                                                            >  2:977565/2‑150 (MQ=255)
                               tCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGa                                                                                                                     >  2:1414390/1‑150 (MQ=255)
                                    gcatccGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGatctgc                                                                                                                <  2:1166621/150‑1 (MQ=255)
                                      atcgGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGatcttcat                                                                                                              <  2:731482/150‑1 (MQ=255)
                                               tgcGCCATAGCTGTCAGGTTGCTTCTGCGTCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAgg                                                                                                     <  2:142513/148‑1 (MQ=255)
                                                 cGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGt                                                                                                   <  2:737812/150‑1 (MQ=255)
                                                 cGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGt                                                                                                   <  2:1375629/150‑1 (MQ=255)
                                                   ccATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTa                                                                                                 >  2:1067199/1‑150 (MQ=255)
                                                         aTGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTc                                                                                                                                               >  1:484595/1‑98 (MQ=255)
                                                         aTGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTc                                                                                                                                               <  2:484595/98‑1 (MQ=255)
                                                              cAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGccagccag                                                                                                                                <  2:685703/108‑1 (MQ=255)
                                                              cAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGccagccag                                                                                                                                >  1:685703/1‑108 (MQ=255)
                                                                gccgCCTCGGCGCAGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCAcc                                                                                        <  2:1256768/143‑1 (MQ=255)
                                                                gccgCCTCGGCGCAGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCAcc                                                                                        >  1:1256768/4‑146 (MQ=255)
                                                                gTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGGATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAgc                                                                                    >  1:1009115/1‑150 (MQ=255)
                                                                gTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAgc                                                                                    <  1:1153259/150‑1 (MQ=255)
                                                                 ttGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAgcg                                                                                   >  2:273112/1‑150 (MQ=255)
                                                                      tctGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCgcagca                                                                              >  2:363840/1‑150 (MQ=255)
                                                                        tGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCagcagc                                                                            <  1:363840/150‑1 (MQ=255)
                                                                                   ggcgTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGccagccag                                                                                                                                >  1:452160/1‑87 (MQ=255)
                                                                                   ggcgTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGccagccag                                                                                                                                <  2:452160/87‑1 (MQ=255)
                                                                                    gcgTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGa                                                                                      >  1:717841/1‑128 (MQ=255)
                                                                                    gcgTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGa                                                                                      <  2:717841/128‑1 (MQ=255)
                                                                                    gcgTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAg                                                                >  2:301282/1‑150 (MQ=255)
                                                                                    gcgTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAg                                                                <  1:459702/150‑1 (MQ=255)
                                                                                           ccacccGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGatctgcat                                                                                                              >  2:1292759/1‑97 (MQ=255)
                                                                                           ccacccGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGatctgcat                                                                                                              <  1:1292759/97‑1 (MQ=255)
                                                                                           ccacccGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTggg                                                         >  1:1086976/1‑150 (MQ=255)
                                                                                           ccacccGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTggg                                                         >  2:601122/1‑150 (MQ=255)
                                                                                               ccGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCtc                                                     >  1:1334212/1‑150 (MQ=255)
                                                                                               ccGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCtc                                                     >  2:731810/1‑150 (MQ=255)
                                                                                                 gCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCGTACATGATCTGCATCTGCTCAGGGTTACGGGCCACCGAGCGCAGCAGCGAAGAGCGCCGGGACTGGGCCtcgt                                                   >  2:567114/1‑150 (MQ=255)
                                                                                                 gCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGCTCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCtcgt                                                   <  2:697540/150‑1 (MQ=255)
                                                                                                  cacaGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCtcgtc                                                  <  2:28938/150‑1 (MQ=255)
                                                                                                     aGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGc                                                        >  2:1244013/1‑141 (MQ=255)
                                                                                                     aGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGc                                                        <  1:1244013/141‑1 (MQ=255)
                                                                                                     aGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGaa                                               <  2:1334212/150‑1 (MQ=255)
                                                                                                     aGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGaa                                               <  1:1414390/150‑1 (MQ=255)
                                                                                                     aGGCTCAGAATGCTCATAACCGGCCAGTCACGTCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGaa                                               <  1:649865/150‑1 (MQ=255)
                                                                                                             aaTGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCCAGCCTGAGACTAGTCGACGTAGCTaa                                       <  2:19910/150‑1 (MQ=255)
                                                                                                              aTGCTCATAACCGGCCAGCCGCGGCAGGGTGTCCTCCTTCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATGTGCTCAGGGGTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGAAGACAAg                                      >  2:575239/1‑150 (MQ=255)
                                                                                                                 cTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGTAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTcgcg                                                                    <  1:51195/117‑1 (MQ=255)
                                                                                                                 cTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGTAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTcgcg                                                                    <  1:51311/117‑1 (MQ=255)
                                                                                                                 cTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGTAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTcgcg                                                                    >  2:51195/1‑117 (MQ=255)
                                                                                                                 cTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGTAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTcgcg                                                                    >  2:51311/1‑117 (MQ=255)
                                                                                                                  tCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCtcgt                                                   >  2:742104/1‑133 (MQ=255)
                                                                                                                  tCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCtcgt                                                   <  1:742104/133‑1 (MQ=255)
                                                                                                                   cATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCGCGTCGCAGGCGCCAAgggtgt                                 >  2:143183/1‑146 (MQ=255)
                                                                                                                    atgaCCGACCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCAt                                <  1:143183/147‑1 (MQ=255)
                                                                                                                    aTAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCAt                                <  1:301282/150‑1 (MQ=255)
                                                                                                                    aTAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCAt                                >  1:500466/1‑150 (MQ=255)
                                                                                                                      aaCCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCAc                                                                                         <  1:1067241/91‑1 (MQ=255)
                                                                                                                      aaCCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCAc                                                                                         >  2:1067241/1‑91 (MQ=255)
                                                                                                                       aCCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTg                                                           >  2:630704/1‑120 (MQ=255)
                                                                                                                       aCCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTg                                                           <  1:630704/120‑1 (MQ=255)
                                                                                                                        ccGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCGCGGCGAAGTGGCCAAGGTTCATGAAg                            >  2:647218/1‑150 (MQ=255)
                                                                                                                         cGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAgg                           >  1:1399447/1‑150 (MQ=255)
                                                                                                                          ggCCAGCCACGGCAAGGGTTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGc                          <  2:339850/150‑1 (MQ=255)
                                                                                                                                acaCGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGCcagcag                    >  2:978117/2‑150 (MQ=255)
                                                                                                                                  aCGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCGGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGACCCCGATCTGAGTAGCCAAGGTTCATGAAGGCCAGCAggg                  >  2:824540/1‑150 (MQ=255)
                                                                                                                                  aCGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGaa                             >  2:123398/1‑139 (MQ=255)
                                                                                                                                  aCGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGaa                             <  1:123398/139‑1 (MQ=255)
                                                                                                                                  aCGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGCCAGCAggg                  >  1:116203/1‑150 (MQ=255)
                                                                                                                                  aCGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGCCAGCAggg                  >  1:311615/1‑150 (MQ=255)
                                                                                                                                                tCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAg                                      <  1:645157/116‑1 (MQ=255)
                                                                                                                                                tCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAg                                      >  2:645157/1‑116 (MQ=255)
                                                                                                                                                tCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGCCAGCAGGGTCATGgtcgcgtcg    >  2:1025406/1‑150 (MQ=255)
                                                                                                                                                tCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGCCAGCAGGGTCATGgtcgcgtcg    >  2:453646/1‑150 (MQ=255)
                                                                                                                                                 ccTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGCCAGCAGGGTCATGgtcgcgtcgc   <  1:601122/150‑1 (MQ=255)
                                                                                                                                                  cTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGCCAGCAGGGTCATGGtcgcgtcgcg  >  2:457549/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTGTACTCCTGCAGGTCACGCTCGCCAGCCAGGCCATACATGATCTGCATCTGCTCAGGGTTACCGGCCACCGAGCGCAGCAGCCACTCGCGCCAGGCCTGGGCCTCGTCGAAGTAGCCAAGGTTCATGAAGGCCAGCAGGGTCATGGTCGCGTCGCG  >  NC_002947/4553885‑4554180

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑