Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,077,916 C→T G323R (GGG→AGG)  putA ← fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,077,9160CT98.7% 242.4 / NA 75G323R (GGG→AGG) putAfused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (35/39);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTC  >  NC_000913/1077670‑1078162
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
ccTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCa                                                                                                                                                                                                                                                     <  2:442680/251‑1 (MQ=255)
 cTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCAt                                                                                                                                                                                                                                                    >  1:540372/1‑251 (MQ=255)
    aaTAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGACGCCAGTTTGATTGAAATGCCCGGCCTTTCAt                                                                                                                                                                                                                                                    >  1:295750/1‑248 (MQ=255)
       aaCAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATg                                                                                                                                                                                                                                               <  1:176993/250‑1 (MQ=255)
       aaCAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGACCGGCCTTTCATAGATGc                                                                                                                                                                                                                                              >  2:558502/1‑251 (MQ=255)
       aaCAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAAAGCGGATGCCGCGCCGACAGTTTTAATGTAATGCCCGGCCTTTCAAAGATGc                                                                                                                                                                                                                                              >  2:474673/1‑251 (MQ=255)
          aaaaCCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCAc                                                                                                                                                                                                                                           <  2:97939/251‑1 (MQ=255)
          aaaaCCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTGTGATTGAAATGCCCGGCCTTTCATAGATGccgc                                                                                                                                                                                                                                           <  1:268771/251‑3 (MQ=255)
               cGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCg                                                                                                                                                                                                                                      >  1:402111/1‑251 (MQ=255)
               cGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGAGCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATGCCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTGCCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCg                                                                                                                                                                                                                                      <  1:291593/251‑1 (MQ=255)
                            gCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTAc                                                                                                                                                                                                                         <  2:397775/251‑1 (MQ=255)
                                gCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGAGCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCg                                                                                                                                                                                                                       <  1:37023/249‑1 (MQ=255)
                                 ccAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGATATCTCCAGGCGATCTGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATTCCACGACCGTTAGACGCTTTACCGATg                                                                                                                                                                                                                    <  2:306544/251‑1 (MQ=255)
                                    gTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTTTTCCATTACCCGGTCATACTGGGCGCGGCTATAAAGCGGATGCAGCGCCGACCGTTTGATTGAAATGGCCGGGCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCg                                                                                                                                                                                                                 >  2:210362/1‑251 (MQ=255)
                                       ccGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGa                                                                                                                                                                                                              <  2:31874/251‑1 (MQ=255)
                                          gCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATc                                                                                                                                                                                                           <  1:349977/251‑1 (MQ=255)
                                           cTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTTGGCGCGGCTATAACGCGGATGCAGCGCCGACAGGTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGCCCGTTAGACGCTTTACCGCTGGCGTGAATcc                                                                                                                                                                                                          >  2:315/1‑250 (MQ=255)
                                                    gagTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGAt                                                                                                                                                                                                 <  2:167933/251‑1 (MQ=255)
                                                       tttttCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCACTTCCATTACACGGACATACTGGGAGCGGGTATAACGCGGATGCAGCGACGACAGTTTGATTGAAATGCTCGGCCTTTCATAGATGCCCCGACCGTTAGACGCTTTAACGATGGCGTGAAAGGCctgctggaagg                                                                                                                                                                                              >  1:378690/1‑246 (MQ=255)
                                                          ttCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCAGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGaa                                                                                                                                                                                            <  1:12097/250‑1 (MQ=255)
                                                              agcagATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTTAAATGCCCGGCCTTTCATAGATGCCCCGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTTCTGATAGGAAACCat                                                                                                                                                                                       >  2:561852/1‑251 (MQ=255)
                                                               gcagATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCata                                                                                                                                                                                      <  1:558502/251‑1 (MQ=255)
                                                                  gATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTCTTAAGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGCAGCTTTACCGAATGCGTGAATCGCCTGCTGATAGGACACCATATAc                                                                                                                                                                                   >  2:126831/1‑251 (MQ=255)
                                                                    tCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATTGCCTGCTGATAGGAAACCATATACGc                                                                                                                                                                                 >  2:243102/1‑251 (MQ=255)
                                                                     cGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGcc                                                                                                                                                                                <  1:561852/251‑1 (MQ=255)
                                                                      gAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCt                                                                                                                                                                               >  1:115612/1‑251 (MQ=255)
                                                                      gAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGAAGCTTTACCGATGGCGTGAATCTCCTGCTGAAAGGAAACCATATACGCCt                                                                                                                                                                               >  2:356793/1‑251 (MQ=255)
                                                                       aGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGGAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAAAGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATTCCACGACCGTTAGACGCTTTACCGATGGCGTTAATCGCCTGTTGATATGAAACCATATACCCCtg                                                                                                                                                                              >  2:536076/1‑251 (MQ=255)
                                                                        gggAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCtgt                                                                                                                                                                             >  1:463147/1‑251 (MQ=255)
                                                                         ggAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCtgtg                                                                                                                                                                            <  1:109144/251‑1 (MQ=255)
                                                                                 ccAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATc                                                                                                                                                                        <  1:170533/247‑1 (MQ=255)
                                                                                 ccAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCt                                                                                                                                                                       <  2:540372/248‑1 (MQ=255)
                                                                                 ccAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTg                                                                                                                                                                      <  1:447112/249‑1 (MQ=255)
                                                                                                  tCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCtt                                                                                                                                                    <  1:474673/250‑1 (MQ=255)
                                                                                                    ggCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGc                                                                                                                                                 <  1:321742/251‑1 (MQ=255)
                                                                                                    ggCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTAACCGGTCATACTGGGCGCGGATATAAAGCGGATGCCGCGCCGCAAGTTTGATTGAAATGCCCGGCCCTTCATAGAAGCCCCGGCCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGAAAAGAAACCATATACCCCTATGCATCTGCGGCGGGCAGCGCGGCtttgc                                                                                                                                                 >  1:345012/1‑248 (MQ=255)
                                                                                                          aaTGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCa                                                                                                                                           <  1:248477/251‑1 (MQ=255)
                                                                                                            tGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCCCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCCCGACCGTTAGACGCTTTTCCGCTGGCTTGAATCCCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGGCAACGCAGCtttcccccgccac                                                                                                                                         >  2:57400/1‑240 (MQ=255)
                                                                                                                   aCCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTgcgg                                                                                                                                                                   <  1:24722/218‑1 (MQ=255)
                                                                                                                   aCCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCACTTTCAATAACCGGGCATACCGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCCCGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGACTGTGAATCTgcgg                                                                                                                                                                   >  2:24722/1‑218 (MQ=255)
                                                                                                                     cAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACTACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCCGCGCGGCGTTGCCCCGCATATCGTAaggg                                                                                                                                >  1:541699/1‑249 (MQ=255)
                                                                                                                              aCTGACTCGCCAGCAGGGTGAGTGATGTCAGACGCGTGTAAAGCTCGTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGACACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAGCGGaaa                                                                                                                       <  2:339099/251‑1 (MQ=255)
                                                                                                                                   cgcgCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGGCATCCTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGACACCATATACGCCTGTGCATCTGCGGCGGGCCACGCGGCTTCGCCCCACACATCGTAAGAGGTACGGAAAACttt                                                                                                                 >  2:620814/1‑251 (MQ=255)
                                                                                                                                          gcagGGTGAGTGATTGCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTcc                                                                                                          <  1:420537/251‑1 (MQ=255)
                                                                                                                                           cagGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCa                                                                                                         <  1:315/251‑1 (MQ=255)
                                                                                                                                           cagGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGACTATCTTAATAGGAAAGGAAACCTTTCTCTTccc                                                                                                         >  2:525705/1‑250 (MQ=255)
                                                                                                                                               gtgagtgaTTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCtt                                                                                                     <  1:37122/251‑1 (MQ=255)
                                                                                                                                                       tttCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCCGCGCGGCTTCCCCCCGCATAACGTAAGAGTAACGGGAACCCTTCCCCTCCCACTTGCgggga                                                                                              >  2:322094/1‑248 (MQ=255)
                                                                                                                                                        ttCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGCCAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGACTGCTGATAGGAAACCCTATACGCCTGTGCATCTGCGGCGGTCGGCGCGGGTTCGCCCAGCATTGCGTAAAAGAAACGGAAACACTTTTCTTCCAACtttgtggnct                                                                                             >  2:601541/1‑242 (MQ=255)
                                                                                                                                                          cAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTg                                                                                           <  1:353177/250‑1 (MQ=255)
                                                                                                                                                             aCGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTa                                                                                       <  1:356793/251‑1 (MQ=255)
                                                                                                                                                              cgcgGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTa                                                                                       >  2:101228/1‑250 (MQ=255)
                                                                                                                                                                     aaaGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAAAGCTTcc                                                                               >  2:158157/1‑251 (MQ=255)
                                                                                                                                                                       aGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGGCAGCGCCGCTTCGCCCAGCATAACGCAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGGTTCcgg                                                                             >  1:58541/1‑250 (MQ=255)
                                                                                                                                                                                    cccGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCa                                                                <  2:115612/251‑1 (MQ=255)
                                                                                                                                                                                             aCTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGaa                                                        <  1:311337/250‑1 (MQ=255)
                                                                                                                                                                                                 ggCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCCGTGACGAACTGCt                                                   >  1:491731/1‑251 (MQ=255)
                                                                                                                                                                                                      ggCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCAccc                                              <  1:483550/251‑1 (MQ=255)
                                                                                                                                                                                                      ggCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGCCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCAccc                                              >  1:185974/1‑251 (MQ=255)
                                                                                                                                                                                                          ataACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGCAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATcc                                          >  1:373998/1‑250 (MQ=255)
                                                                                                                                                                                                          ataACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCa                                          <  1:11546/251‑1 (MQ=255)
                                                                                                                                                                                                               gcgGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGgcgc                                     <  1:204906/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGCTAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCAt                                   >  1:54973/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCAt                                            <  1:190097/242‑1 (MQ=255)
                                                                                                                                                                                                                 ggATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCAt                                            >  2:190097/1‑242 (MQ=255)
                                                                                                                                                                                                                    tGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCGTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGc                                >  1:661248/1‑251 (MQ=255)
                                                                                                                                                                                                                     gcagcgatgaatggttgattgaAATGCCCGGCCTTTCATAGATGCCCCGACCGATAGACGCTTTACCGATGGCGTCAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCCTCGCCCACCCTATCGTAAGAACAACGGAAACCTTTCTCTTCCAGCTTGCCGGCATTGGCCAACGCCTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACACATCATGCGCCTCGcc                               >  2:78865/15‑251 (MQ=255)
                                                                                                                                                                                                                     gCAGCGCCGACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGcc                               <  2:463147/251‑1 (MQ=255)
                                                                                                                                                                                                                             gACAGTTTGATTGAAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCcac                       <  1:545911/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gattgaAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCTTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGTATCGCCATATCCACACCTTTGc               <  2:295750/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gattgaAAGGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTCCCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGc               <  1:635116/251‑1 (MQ=255)
                                                                                                                                                                                                                                      attgaAATGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGGAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTTGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCg              >  2:33190/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aaaTGCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATc          <  2:402111/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gCCCGGCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCg      >  1:587737/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ggCCTTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCCTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTc  >  1:418460/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
CCTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGAGTTTTTCCAGCAGATCGAGGGAGATCTCCAGGCGATCGGACTCTTCGGCGTCAATGTTGATACCAATATCGTACTGACGCGCCAGCAGGGTGAGTGATTTCAGACGCGGGTAAAGCTCTTCCATTACCCGGTCATACTGGGCGCGGCTATAACGCGGATGCAGCGCCGACAGTTTGATTGAAATGCCCGGCCCTTCATAGATGCCACGACCGTTAGACGCTTTACCGATGGCGTGAATCGCCTGCTGATAGGAAACCATATACGCCTGTGCATCTGCGGCGGTCAGCGCGGCTTCGCCCAGCATATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTC  >  NC_000913/1077670‑1078162

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