Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 352,830 A→G 15.2% M301V (ATG→GTG)  prpE → propionate‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779352,8300AG15.2% 183.8 / 3.4 78M301V (ATG→GTG) prpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (13/51);  new base G (11/3);  total (24/54)
Fisher's exact test for biased strand distribution p-value = 6.44e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACG  >  NC_007779/352596‑353060
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
cgGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTggg                                                                                                                                                                                                                                                 >  1:942005/1‑238 (MQ=255)
   gTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGa                                                                                                                                                                                                                                               <  1:724484/238‑1 (MQ=255)
     aCCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACt                                                                                                                                                                                                                                             <  1:1014631/238‑1 (MQ=255)
     aCCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACt                                                                                                                                                                                                                                             <  1:90523/238‑1 (MQ=255)
     aCCGGTGGCGTGGCTCGAATCCACCCACGCCGCCTGCGTTCGCTACGTTTCCGGGACGGACTGCATACCTAAAGGCGTGAGGCGTGACGTCGGCGGAGAGGCGGTGGCGCTGGCGACGTCGATTGGCAACTTTTATAGCGGCAAAGCGGGCAGCGTGTTCTTTGGCGCGGCGGTTATCGGCGGGGGGGCGGGTCATTCGTGTTTCGTTTACGCGCCGCTGCTGGCGGGGATGGCGTCt                                                                                                                                                                                                                                             <  2:882862/238‑1 (MQ=255)
         gtggcgtggcTGTAATCCTACGAAAGCTCCTGCATTCTCTTCACTACCGGCAGCACCGGAAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATCGGCTGTGGAGCTGGCGACCTCTATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCCTCGGATATCGGCTTGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGTGATGGCGACTATCg                                                                                                                                                                                                                                        <  2:1341871/238‑1 (MQ=255)
               tggcTGGAATCCAACGAAACCTCCTGCATTCTCTACCCTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACg                                                                                                                                                                                                                                  <  2:1150383/238‑1 (MQ=255)
                   tGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGa                                                                                                                                                                                                                             <  2:701768/239‑1 (MQ=255)
                          cAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGa                                                                                                                                                                                                                     <  1:203202/239‑1 (MQ=255)
                           aaCGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGAc                                                                                                                                                                                                                   <  1:1606005/239‑1 (MQ=255)
                                aaCCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTgg                                                                                                                                                                                                               <  2:783125/238‑1 (MQ=255)
                                    tcctGCCTTCTCTACGCGTCCGGCACGACCGGCAAACCGAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑a                                                                                                                                                                                                        <  2:1505514/239‑1 (MQ=255)
                                    tcctGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑a                                                                                                                                                                                                        <  1:670048/239‑1 (MQ=255)
                                              tctgcgcTTCCGGCCCGCCCGGACAACCTAAAGGCGTTCAGCGTGACGTCGGCGGACATGCGGTGGCGATGGCGGCCTCGATGGACGCCATTGTTGGCGGCGAGGCGGGTAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCgtg                                                                                                                                                                                              <  2:1318064/233‑1 (MQ=255)
                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCgtg                                                                                                                                                                                              <  1:681842/239‑1 (MQ=255)
                                              tctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCgtg                                                                                                                                                                                              <  1:1211674/239‑1 (MQ=255)
                                               ctACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCgtgt                                                                                                                                                                                             <  1:247651/239‑1 (MQ=255)
                                                   acTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTgtggtg                                                                                                                                                                                         <  1:452571/239‑1 (MQ=255)
                                                    cTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGGGGCGGCTATCGGTTTCGAAGGGTTGCCGACCCGGCCGG‑‑ACCgcggcggggggggg                                                                                                                                                                                        >  1:628707/1‑231 (MQ=255)
                                                     ttCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGa                                                                                                                                                                                       <  1:1769994/239‑1 (MQ=255)
                                                     ttCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGa                                                                                                                                                                                       <  1:339572/239‑1 (MQ=255)
                                                          gCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGTTTTCGAAGGATTGCCGGCCTGGCCGG‑‑GCTGCGGCGTGTGGGGGGCaaa                                                                                                                                                                                   >  1:1248921/1‑237 (MQ=255)
                                                            aCGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAAtcgt                                                                                                                                                                              <  1:637543/239‑1 (MQ=255)
                                                               aCCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAGGGGGGCTATCGGTTACGAAGGATTGCCGACCCGGCCGG‑‑GCTGCGGCGGGTGGGGGGCAAacgtcgg                                                                                                                                                                           >  1:1701352/1‑238 (MQ=255)
                                                                ccGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCgag                                                                                                                                                                          <  1:270229/239‑1 (MQ=255)
                                                                 cGGCAAGCCGTAAGGCGTGCAGCGGGGCGTCGGCGGATATGGGGGGGCGCTGGCGACGTCGATGGACACCATTTTTGGCGGCAAGGCGGGCAGCGTGTTCTTTGGCGCATCGGATATCGGCGGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCgag                                                                                                                                                                          <  2:51617/238‑1 (MQ=255)
                                                                     aaaCCGAAAGGGGTGCGGCGTGGCGTCGGCGGATATGCGGGGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAat                                                                                                                                                                      <  2:942005/238‑1 (MQ=255)
                                                                        ccTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTgcg                                                                                                                                                                                                   <  1:2877/208‑1 (MQ=255)
                                                                          tAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGGATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAgg                                                                                                                                                               <  1:1416202/239‑1 (MQ=255)
                                                                                cGTGCGGCGTGGCGGCGGCGGAGATGCGGTGGCGCTGGCGACCTCGATGGGCACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGGCTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGc                                                                                                                                                          <  2:1047878/238‑1 (MQ=255)
                                                                                cGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGc                                                                                                                                                          <  1:842934/238‑1 (MQ=255)
                                                                                  tGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGTTTTCGAAGGATTGCCGACCCGGCCGG‑‑ACTGCGGCGGGGGGTGGGAAATCGGCGGGAAAAATTAGGTTAGCCgg                                                                                                                                                      >  1:332979/1‑239 (MQ=255)
                                                                                  tGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGGGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGGGTATCGGTTTCCCAAGGATGCCGAGCCTGCCCGG‑‑CTGCGGGGGGTGGGGGGTGATAGTCGAGAAATAACAgggatgccgg                                                                                                                                                       >  2:1746422/1‑231 (MQ=255)
                                                                                   gCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGTTTTCGAAGGATTGCCGACCTGGCCGG‑‑GCTGCGGCGGGTGGGGGGCAATCGGCGGGAAATATCAAGGTAGCCggg                                                                                                                                                     >  1:1889971/1‑238 (MQ=255)
                                                                                     aGCGGGGCGTCGCGGGGTATGCGGTGGCGCTGGCGCCCTCGATGGGCACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTGGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGAt                                                                                                                                                   <  2:1808958/238‑1 (MQ=255)
                                                                                       cGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGTTTACGAAAGATTGCCGACCTGGCCGG‑‑ACTTCGGCGGGTGGGGGACAATCGTCGAGAAATATCAGGTTAGCCCGATGtt                                                                                                                                                >  1:377772/1‑239 (MQ=255)
                                                                                            cGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCa                                                                                                                                            <  1:544751/238‑1 (MQ=255)
                                                                                                                cgcTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGGGTGGGGGGCAATCGTCGGGAAATATCAGGTTAGGCGGATGTTCTCAGCGCCGACCGCCATTCGCGt                                                                                                                        >  1:972237/1‑238 (MQ=255)
                                                                                                                 gtgggCGTCCTGGAGGGACCCCATTTTGGGCGGCAAGGCGGCCGGCTGTTTCTTTGGCGTGTCGGATATGGTCTGGGTGGGGGGGCTTTCTTTTATCGTTTAGGCGCCGCTGCGGGCGGGGGTGGGGTTTATCTTTTACGAGGTATGCGCGCCCGGGCCGG‑‑ACTGGGGTGTGGGGTGGCCAAACTACGAGAAATATCAGGTTAGCCGGATGTTCTCGGCGCCGCCCGCCATTCGCGTGc                                                                                                                      <  2:1236685/236‑1 (MQ=255)
                                                                                                                   tGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCt                                                                                                                     <  1:1110188/238‑1 (MQ=255)
                                                                                                                       gACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGaaa                                                                                                               <  1:1553573/238‑1 (MQ=255)
                                                                                                                             gATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGGGGACAATCGTCGAGAAATATCAGGGTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTc                                                                                                         >  1:1418429/1‑238 (MQ=255)
                                                                                                                               tGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTccc                                                                                                       <  1:1183079/238‑1 (MQ=255)
                                                                                                                                    acCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGGGTGGTGGGCAATCGTCGAGAAATATCAGGGTTGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAAATCCCCACCg                                                                                                 >  1:94466/1‑238 (MQ=255)
                                                                                                                                      cATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGGGGGGGGGATGGTTTTCGAAAGGTTTCCGCCCCTGGCCGG‑‑ATGGGGGGGGGGGTGGGTAACTCGCTCGAAAAAATCCGGTTAGCGGGGGGTGTGCACCCCCGACCCCCCTTCGTGTTTCGAAAAAATtccctccc                                                                                                  >  2:585874/1‑236 (MQ=255)
                                                                                                                                         ttttGGCGGCAAAGCGGGCGGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGaaa                                                                                           <  1:1602479/239‑1 (MQ=255)
                                                                                                                                                  taaggggggCGGGGTGTTCTTTTGGGCGGCGGTTTTCGGGTGGGGGTTGGGGTTTTCTTATATCGTTGACGGGCCGCTGCGGGCGGGGATGGCGACTATCGTTTCCGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCaaa                                                                                   <  2:993420/232‑1 (MQ=255)
                                                                                                                                                              cGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAAtcgtcg                                                                                                                                                                            <  1:1875953/143‑1 (MQ=255)
                                                                                                                                                               gtgtTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGCCCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctctc                                                                        <  1:1705987/236‑1 (MQ=255)
                                                                                                                                                               gtgtTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctctc                                                                        <  1:29810/236‑1 (MQ=255)
                                                                                                                                                                 gtTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCTTTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGc                                                                   <  1:55088/239‑1 (MQ=255)
                                                                                                                                                                    cTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTgg                                                                <  1:1783885/239‑1 (MQ=255)
                                                                                                                                                                     ttttGCGCAGCGGATTTGGGCTGGGTGGGGGGGCTTTCGTATATCGTTTGCGCGCCGCGTCTGGCGGGGAGGGCGGTTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑CCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCGGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGa                                                               <  2:94466/239‑1 (MQ=255)
                                                                                                                                                                      tttGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGaa                                                              <  1:1608644/239‑1 (MQ=255)
                                                                                                                                                                       tgccgcgTCGGTTGTCGGGGGGGGGGTGGGGATTTCGTTAGTCGTTTACGCGCCGCTGCGGGGGGGGGGGGCGACTATGTTTAACGAAGGGTTGCCGGCCTGGCCGG‑‑ACGGCGGCGTGTGGGGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAg                                                             <  2:1318410/236‑1 (MQ=255)
                                                                                                                                                                        tGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTgg                                                                                                                                                                                                               <  2:1599211/102‑1 (MQ=255)
                                                                                                                                                                        tGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTgg                                                                                                                                                                                                               >  1:1599211/1‑102 (MQ=255)
                                                                                                                                                                          ggcATCGGAGAGGGGGGGGGTGGGGGGACTATCGGATATCGTGGACGGGCCGCGGCGGGCGGGGAGGGCGCCTATCGTTTACGAAGGATGCCCGCCTGGGCGGG‑‑CCGGCGGCGGGTGGGGGACAAGCGTGGAGAAAGACCAGGTAAGCCGGATGCTCGCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTg                                                           <  2:1379803/237‑1 (MQ=255)
                                                                                                                                                                               cGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTAt                                                     <  1:1666800/239‑1 (MQ=255)
                                                                                                                                                                                  atgtCGGGTGGGTGGTGGGGCATTCGTAGATCGTGTACCCGCCGCGGCTGGCGGGGATGGCGACTGTCGTTTACGAGGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGGGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATct                                                   <  2:628707/235‑1 (MQ=255)
                                                                                                                                                                                  atatCGGGTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATct                                                   <  1:1350284/238‑1 (MQ=255)
                                                                                                                                                                                  atatCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATct                                                   <  1:1865928/238‑1 (MQ=255)
                                                                                                                                                                                                   gggCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTgg                                  >  1:1863866/1‑238 (MQ=255)
                                                                                                                                                                                                        ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGACTATCGTTTTTGAATGATTGTCCCCCTGGCCGGCCACTGGGCGGGGGGGGGAAATCCGCGTGAAATATCAGGGTTACCCGGGGTTCTCCGCGCCCCCCCCCACCTGCGGGCTCGAAAAAATACCCACCCCCGCGATTCGCCAACACGCTCTCCCTCTCCTGGGATTGCTCTATTCGGGCGGAGAAACGCAGGGCgcg                             >  2:1102907/1‑236 (MQ=255)
                                                                                                                                                                                                        ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGGCGTTTGTTTTCTAAGGCTTGACGGCCCTGCCCG‑‑GGTGCGGCGGGTGGGGGGAAATCGTCGCGAAATATCAGGTTTGCCGGGATGTTTCTCGGCCGACCGCCATTCGCGTGGTGCAAAAAATTCCTTACTCGCGAATTAGCAAAAATCAGATCTCGCCCCTGGAACGGCTCCTTCTGGCTGGAGCACACCCGGACgga                             >  2:43592/1‑236 (MQ=255)
                                                                                                                                                                                                        ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAg                             >  1:1604106/1‑238 (MQ=255)
                                                                                                                                                                                                             atatCGTTTACGCGCCGCTGCTGGCGGGGAAGGGGACCATCGCTTACGAAGGGGTGCCGACCTTGGCGG‑‑ACTGCGGCGTGTGGTGGGCAATCATCGAGAAATATCAGGTTAGGCGGGAGGTCTAAGCCCCGAGCCCCAATCGCGTGCTGGAAAAATTCCCCATCCCTGAAAATCGCAAACACGAACTCCCCGCGGAGGAAGAGCTCTTACTGGCTGGGAAAACGCTGGACGAGCCgcc                        >  2:1724134/1‑236 (MQ=255)
                                                                                                                                                                                                                   tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCCAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCCGGACGAGCCGACCGCCCGt                 >  1:1885557/1‑239 (MQ=255)
                                                                                                                                                                                                                     tACGCGCCGCTGCTGGCGGGGATGGGGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑GCTGCGGGGTGTGGGGGGCAATCGTTGAGAAATATCTGGGTAGCCGGGTGTTTTCCGCGCCGACCGCCCTTCGCGTGCTGAAAAAAATTCCTACCGCTGAAAATCGCAAACACGATCTCTCGTCGCTGGAAATGCTCTTTTTGGGTGGGGAACCCGTGGACGGGCCGACCGCCAGTTg               >  2:38107/1‑239 (MQ=255)
                                                                                                                                                                                                                     tACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑GCTGCGGCGTGTGGGGGAGCATCGTCCAGAAAAATCCGGGTAGCCGGGTGTTCTCCGCGCCGACCGCCATTCGCGTGCTGAAAAAATTTCCCACCGCCGAAGATCCCCAACAACGTCTCTCCTCTCTGAAAGTGCTCTATCTGGGTTGGGAACAGATGGGCGAGCCGACCGCCacgtg               >  2:1657484/1‑235 (MQ=255)
                                                                                                                                                                                                                     tACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGc                                                                                                                                                          >  1:1274448/1‑105 (MQ=255)
                                                                                                                                                                                                                     tACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGc                                                                                                                                                          <  2:1274448/105‑1 (MQ=255)
                                                                                                                                                                                                                      ccgcgCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTgg              <  1:1751994/238‑1 (MQ=255)
                                                                                                                                                                                                                       gcggccggTGGTGGCGGGGGTGGGCGTTTGTTGTTTCTGGAGGTTTGCCGCCTTGGCCGG‑‑CTGGGGGGGGGGGGTGAACAATCGTCGAAAAATATCGGGTTGCCGGTTTTTTCCGGGCCCCGCCCGCCATTCGCGTGCTGAAAAAATTCCCTCCCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGAAAGGGCTCTATCTGGCTGGAGACCCGCTGGACGAGCCGACCGCCAGTTgg              <  2:1701352/236‑1 (MQ=255)
                                                                                                                                                                                                                           ccgcTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGGGCTCTATCTGGCTGGAGAAACGCTGGACGAGCCGACCGCCAGTTGGGTGAg         >  1:1460048/1‑239 (MQ=255)
                                                                                                                                                                                                                            cgcTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGGTAGCCGGATGTTCTCAGCACCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGc        >  1:181336/1‑239 (MQ=255)
                                                                                                                                                                                                                             gctgctGGGGGGGAGGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCCCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCa       <  1:626171/239‑1 (MQ=255)
                                                                                                                                                                                                                                  tgGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGCCCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACg  <  1:1172801/239‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
CGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGG‑‑ACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACG  >  NC_007779/352596‑353060

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: