Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,082,957 G→A 5.7% R45C (CGT→TGT)  purK ← 5‑(carboxyamino)imidazole ribonucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,082,9570GA 5.7% 210.2 / 2.8 70R45C (CGT→TGT) purK5‑(carboxyamino)imidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base G (34/32);  new base A (2/2);  total (36/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.41e-01

GCGCAGCGATTCGGCGCTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCC  >  NC_002947/6082811‑6083099
                                                                                                                                                  |                                                                                                                                               
gcgcAGCGATTCGGCGCTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAg                                                                                                                                              <  1:891805/150‑1 (MQ=255)
  gcAGCGATTCGGCGCTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGAGGGCCTGACCGCCGTAGTCGGCACGCAGGt                                                                                                                                            >  2:206261/1‑150 (MQ=255)
   cAGCGATTCGGCGCTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTg                                                                                                                                           >  1:562252/1‑150 (MQ=255)
     gCGATTCGGCGCTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCt                                                                                                                                         <  1:1133191/150‑1 (MQ=255)
      cGATTCGGCGCTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTc                                                                                                                                        <  2:975676/150‑1 (MQ=255)
         ttCGGCGCTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGcc                                                                                                                                     >  1:1017285/1‑150 (MQ=255)
              cgcTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAggg                                                                                                                                <  2:799459/150‑1 (MQ=255)
               gcTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAgggg                                                                                                                               >  2:144369/1‑150 (MQ=255)
               gcTGGGGTACACGGGCACGAACGGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAgggg                                                                                                                               <  2:802727/150‑1 (MQ=255)
                 tGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAg                                                                                                                                  <  1:789011/145‑1 (MQ=255)
                 tGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAg                                                                                                                                  >  2:789011/1‑145 (MQ=255)
                        cacGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGgcgcagcgca                                                                                                                      <  1:987138/150‑1 (MQ=255)
                                            ggAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGa                                                                                                 <  2:169309/150‑1 (MQ=255)
                                             gAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAggg                                                                                                                                <  2:1201668/119‑1 (MQ=255)
                                             gAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAggg                                                                                                                                >  1:1201668/1‑119 (MQ=255)
                                             gAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAg                                                                                                >  2:65186/1‑150 (MQ=255)
                                                      ggTCTCGGCGGGGACGCTTGCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGaa                                                                                       <  1:206261/150‑1 (MQ=255)
                                                       gTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGT‑GTCGAGGAAGGCGAAg                                                                                      <  2:1223107/150‑1 (MQ=255)
                                                       gTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAg                                                                                      <  2:1017285/150‑1 (MQ=255)
                                                       gTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAg                                                                                      <  2:562252/150‑1 (MQ=255)
                                                          tcGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTc                                                                                   >  2:636382/1‑150 (MQ=255)
                                                                 ggACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCa                                                                            <  1:144369/150‑1 (MQ=255)
                                                                  gACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAg                                                                           >  1:1372055/1‑150 (MQ=255)
                                                                       tttCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCGCAGGGGCGCAGCGCAGGCGTCCGCTCCCGG‑GTCGAGGAAAGCGAACTTCATGCCCAGCgggg                                                                      >  2:386009/1‑150 (MQ=255)
                                                                       tttCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCg                                                                                                              <  2:792150/111‑1 (MQ=255)
                                                                       tttCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCg                                                                                                              >  1:792150/1‑111 (MQ=255)
                                                                       tttCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCgggg                                                                      >  1:1314852/1‑150 (MQ=255)
                                                                       tttCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCgggg                                                                      >  1:455077/1‑150 (MQ=255)
                                                                       tttCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCgggg                                                                      >  1:392361/1‑150 (MQ=255)
                                                                               tCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGaa                                                                                       >  1:508861/1‑125 (MQ=255)
                                                                               tCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGaa                                                                                       <  2:508861/125‑1 (MQ=255)
                                                                                   aGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGAAGCACAGCCG‑CCGCTTCGGGCGACGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAgg                                                          >  1:1600715/1‑150 (MQ=255)
                                                                                      tGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCa                                                       <  2:392361/150‑1 (MQ=255)
                                                                                        aCCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGc                                                     <  2:455077/150‑1 (MQ=255)
                                                                                          cAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAg                                                                                            <  1:1011196/109‑1 (MQ=255)
                                                                                          cAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAg                                                                                            >  2:1011196/1‑109 (MQ=255)
                                                                                                  cTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGaa                                                                                       <  1:335008/106‑1 (MQ=255)
                                                                                                  cTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGaa                                                                                       >  2:335008/1‑106 (MQ=255)
                                                                                                            cAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGaa                                                                                       <  1:181622/96‑1 (MQ=255)
                                                                                                            cAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGaa                                                                                       >  2:181622/1‑96 (MQ=255)
                                                                                                            cAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGgccgcc                                 >  2:557363/1‑150 (MQ=255)
                                                                                                              gCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAg                                          <  2:1448119/139‑1 (MQ=255)
                                                                                                              gCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAg                                          >  1:1448119/1‑139 (MQ=255)
                                                                                                                    gCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAg                                                                                      <  2:1522090/89‑1 (MQ=255)
                                                                                                                    gCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAg                                                                                      >  1:1522090/1‑89 (MQ=255)
                                                                                                                    gCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTa                         >  1:1124737/1‑150 (MQ=255)
                                                                                                                    gCAGGTGGTCCTGATCGCCGTAGTCGGCACACAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTa                         >  2:367364/1‑150 (MQ=255)
                                                                                                                    gCAGGTGGTCCTGATCGCCGTAGTCGGCACACAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTa                         >  2:367645/1‑150 (MQ=255)
                                                                                                                      aGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTaca                       <  2:1314852/150‑1 (MQ=255)
                                                                                                                       ggtggtCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCTGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTacac                      <  2:1124737/150‑1 (MQ=255)
                                                                                                                             ccTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCt                >  2:1429884/1‑150 (MQ=255)
                                                                                                                              cTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCtt               >  2:799927/1‑150 (MQ=255)
                                                                                                                               tGATCGCCGTAGTCGGCACGCAGGTGCNCGCCCAGGGGCGNAGCGCANGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTAATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCCGCTGGCCGCCACCGATTACACCGATCTTc              >  1:1517251/1‑150 (MQ=255)
                                                                                                                                gATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCa             >  2:833684/1‑150 (MQ=255)
                                                                                                                                 aTCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTc                                                                                                   <  2:914878/63‑1 (MQ=255)
                                                                                                                                 aTCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTc                                                                                                   >  1:914878/1‑63 (MQ=255)
                                                                                                                                 aTCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCAt            >  1:839774/1‑150 (MQ=255)
                                                                                                                                 aTCGCCGTAGTCGGCACACAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCAt            <  1:367645/150‑1 (MQ=255)
                                                                                                                                 aTCGCCGTAGTCGGCACACAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCAt            <  1:367364/150‑1 (MQ=255)
                                                                                                                                      cGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATgg          <  1:569300/147‑1 (MQ=255)
                                                                                                                                      cGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATgg          >  2:569300/1‑147 (MQ=255)
                                                                                                                                        tAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGa     <  1:1214991/150‑1 (MQ=255)
                                                                                                                                          gTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGc   <  1:557363/150‑1 (MQ=255)
                                                                                                                                           tCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCaa                                                            <  1:1471501/92‑1 (MQ=255)
                                                                                                                                           tCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCaa                                                            >  2:1471501/1‑92 (MQ=255)
                                                                                                                                           tCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGcc  >  1:137836/1‑150 (MQ=255)
                                                                                                                                             ggCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGaa                                                                                       >  1:1259272/1‑63 (MQ=255)
                                                                                                                                             ggCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGaa                                                                                       <  2:1259272/63‑1 (MQ=255)
                                                                                                                                              gCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTgg                                      >  2:1211557/1‑111 (MQ=255)
                                                                                                                                              gCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTgg                                      <  1:1211557/111‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCGCAGCGATTCGGCGCTGGGGTACACGGGCACGAACTGCGAGAGGAAGGCCACGGTCTCGGCCGGGACGCTTTCGAACTCGAAGGTGACCAGGTCGACTTCGTCGGCCAGCTGACGCAGGTGGTCCTGATCGCCGTAGTCGGCACGCAGGTGCTCGCCCAGGGGCGCAGCGCAGGCGTCCGGTGCCGG‑GTCGAGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCC  >  NC_002947/6082811‑6083099

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