Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,632,711 G→A R454C (CGT→TGT)  guaB ← inosine 5'‑monophosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,7110GA97.5% 107.6 / NA 40R454C (CGT→TGT) guaBinosine 5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (26/13);  total (26/14)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAG  >  NC_000913/2632493‑2632945
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
aGTCCAGAATGCGGCTGCGATGCTTCTGAATGTTTTCCGTCATTGACGCTTATTCCGAGGAAAGTGAAACAGATAATATAAATCGACCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACAGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACATACTCCGCTTTAGTACACAGTTCGTCGATaagaacggaagagaggcggg                                                                                                                                                                                                             >  2:43928/1‑232 (MQ=255)
aGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGGCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATa                                                                                                                                                                                                                                <  1:43928/232‑1 (MQ=255)
 gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGAcc                                                                                                                                                                                                            <  1:118108/251‑1 (MQ=255)
                 cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAgg                                                                                                                                                                                            >  1:249504/1‑251 (MQ=255)
                   atgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGAACCACAGCCGGTCAGACCCATACAGGAGCGCAGgcc                                                                                                                                                                                          >  2:76498/1‑251 (MQ=255)
                         aTGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTCCAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCAt                                                                                                                                                                                    >  1:350096/1‑251 (MQ=255)
                              tGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATAAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATctgct                                                                                                                                                                               >  1:347167/1‑251 (MQ=255)
                              tGTTTTCCGTCATAGACGCTTATTCAGAGGCAAGGGCCAAAGATACTATAAATCGACCGACATGAAGTCGGGCGAAGAGAATAAGGAGCCACGACAGTAGTTCGGGGACTCTTTAGTAATGGACACGTCGTGACAGTGGCTTTCCTGAATACCCGCACCGCTGATACGTACACACTCCGCTTTAGTACACAGTTCGTCGATCGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGTCCACCCATatgct                                                                                                                                                                               >  2:49200/1‑251 (MQ=255)
                               gTTTTCCGTCATTGACGCTTAGTCCGAGGCAAGTGAAACAGATAATATAAACCGCCCTACATGAAGTCGGGGGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGGCGATAGTACCACAGCCGGGCAGACCCATACAGGAGCGCAGGCCGCCCATctgct                                                                                                                                                                               <  2:249504/250‑1 (MQ=255)
                                   tCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTg                                                                                                                                                                           >  1:173995/1‑250 (MQ=255)
                                     cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAAt                                                                                                                                                                        >  1:330213/1‑251 (MQ=255)
                                        cagtGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGa                                                                                                                                                                      <  1:8947/247‑1 (MQ=255)
                                            gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGAGCAGGCCGCCCATCTGCTGGTGAATGATctct                                                                                                                                                                 >  1:232955/1‑251 (MQ=255)
                                                          gggaaGTGAAACAGATAATATAACTCGCCCGCATGGAAGTCGGGCGAAGAGAACCAGGAGCCCAGGCGTTAGGTCGGGGATTCTTTAGTAATGGTCACGTCTTGAACGGGGCTTTCCTGAATGCCCGCTCCGCTGATACGTACAAACTCCGCTTTAGTCCACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCttt                                                                                                                                                   <  2:367204/248‑1 (MQ=255)
                                                                  aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCt                                                                                                                                            <  1:191603/250‑1 (MQ=255)
                                                                  aaaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCAGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCt                                                                                                                                            >  2:191603/1‑250 (MQ=255)
                                                                      aGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGgccgcc                                                                                                                                                                                       >  1:376893/1‑203 (MQ=255)
                                                                      aGATAATATAAAACGCCCGACATGAAGTCGGGCGAAGAGAATCAGTAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGTCCCATACAGGAGCGCAGgccgcc                                                                                                                                                                                       <  2:376893/203‑1 (MQ=255)
                                                                       gATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGACCAGACGGTAGTTCGGGGACTCTTTAGCAATGGTAACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGATGCTACATAAAAACTCCGCTTTAGTACACAGTTCGTCGCTAGTACCACAACCGGTCAGACCCATACCGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTAcgcg                                                                                                                                       >  1:167858/1‑250 (MQ=255)
                                                                             ataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCg                                                                                                                                >  1:179590/1‑251 (MQ=255)
                                                                                        tccATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTAGGGGACTCTTTAGTACTGGTCAAGTCGTGAAAGTGGCTTTCCTGAATGCCCTCACCGCTGAGACGTAAAAACTCCGCTTTCGTACACAGTTCGTCGATAGTACCACAGCCAGTCAGACCCATACAGGAGCGCAGACCGCACCTCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTCTAGGCTACGCGCCATTCGATACCTTCCgg                                                                                                                     >  2:237269/3‑251 (MQ=255)
                                                                                         aCATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGc                                                                                                                    >  1:298898/1‑251 (MQ=255)
                                                                                                    gggCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGt                                                                                                          >  1:208964/1‑251 (MQ=255)
                                                                                                            agaATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCg                                                                                                  >  1:193345/1‑251 (MQ=255)
                                                                                                              aaTCAGGAGCCCGGACGTTTGTGTGGGGACTCTTTAGTAGTGGTCACGTCGAGAACGTGGCTTTCCTGAATGCCCGCCCCGCTGATTCGTCCAAACTCCGCTTTAGTACACCTTTCGCTGTTAGTACCACAGCCGGTCAGACCCATTCAGGAGCGCATGCCGCTCATCTGCGGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTTTCGTCAGCGtt                                                                                                <  2:86224/251‑1 (MQ=255)
                                                                                                                          agccgggaGTTCGGGTACTCTTTATTAATGGTCACTTCGTGAACGTGGCTTTCCTTAATGCCCGCTCCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTGCCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAAGGATCTCTTTCAGGCGGCCTTTATAGGCTACGCGACCTTCGATTCCTTCCGGCACCAGTTTGTCGGCAGCGGTATCGGTCTGGaa                                                                                    <  2:193345/244‑1 (MQ=255)
                                                                                                                                                   tAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGAc                                                           >  1:220163/1‑251 (MQ=255)
                                                                                                                                                    aaTGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGAc                                                           >  1:135128/1‑250 (MQ=255)
                                                                                                                                                                gtGAACGTGGCTTTCCTGAATGCCCGCCCCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTGCCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCGTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATcg                                                       <  1:358656/242‑1 (MQ=255)
                                                                                                                                                                  gAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGa                                             >  1:115000/1‑250 (MQ=255)
                                                                                                                                                                    ccGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCTATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCATTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAAcc                                          <  2:298898/250‑1 (MQ=255)
                                                                                                                                                                            ttCCTGAATTCCCTCCCCACTGATACGTACAAACCCCGCTTTAGGACACAGTTCGCCGATAGTCCCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATAcc                                   <  2:135128/249‑1 (MQ=255)
                                                                                                                                                                                     gCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCAATCTGCTGGGGAATGATCTCTTTCAGGCGACCTTTATAGGCTACACGACCTTCTATACCTTCCGGCACCATTTTGTCGGCAGCGTTATCGCTCTGGCAATAACGGTCAGAGGCACCTTTGGACTACGCTCCCCGGGAACCCATCCCACGGAAAgctt                        >  2:286388/1‑248 (MQ=255)
                                                                                                                                                                                                tGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACgg             >  2:195091/1‑251 (MQ=255)
                                                                                                                                                                                                    acttACAAACTCCTCTTCAGTACGCAGTTCGTCTATATAACCACAGCCGGTCAGACCCAGACAGTAGCGCAGACCGCCCATCTGCTGTTGAAGGATCTCTTTCAGTCGACCTTTATAGGCTACTCGACCTTCTATACCTGCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTTGAAATAACGGTCAGAGGAACCTTTGTACCTCGCTCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGcca          <  2:350096/247‑2 (MQ=255)
                                                                                                                                                                                                      gTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTgg       >  1:126144/1‑251 (MQ=255)
                                                                                                                                                                                                      gTACAAACTCCGCTTTAATACACAGTTAGTCGATAGTACTACAGCCGGTCCGACCCATACAGGAACGCAGGACGCCCATCTGCTGGTGCATGATCTCTTTCAGGCGCCCTTTATAGGCTCCACGCCCTTCGATAACTTCAGGCACCAGTTGATCGGCAGCGTTATCGCGCTGGCAGTAACGGTCAGGGAAACCTTTGTGCATCGCGACCAGTGAACCCATACCACGGTAAGATTTGTAAGAACGACCCTgg       >  2:190083/1‑251 (MQ=255)
                                                                                                                                                                                                        acaaacgtcGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCTCGCTAAGATTTGTAAGAACGGCCCTGGt      <  1:117063/242‑1 (MQ=255)
                                                                                                                                                                                                           aaCTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGTTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCATGGtcg    >  1:465/1‑247 (MQ=255)
                                                                                                                                                                                                           aaCTCCGCTTTAGCCCAAAGTTCGTCGATAGTACCACAGCCGGTCAGAACCATACAGGAGCGCCGTCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTAAGCGACCTTCGATACCTGCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGACATAACTGTCGGAGGAACCTTTGGACCTCGCGCCCAGGGAAACATTCACACGGTAAGAttttgaggataggccctggcaggt  >  2:159132/1‑230 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
AGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAG  >  NC_000913/2632493‑2632945

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