Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,373,078 T→C T167A (ACC→GCC)  nanA ← N‑acetylneuraminate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,373,0780TC100.0% 231.8 / NA 61T167A (ACC→GCC) nanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/31);  total (30/31)

GATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTG  >  NC_000913/3372941‑3373225
                                                                                                                                         |                                                                                                                                                   
gatgTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAgcgc                                                                                                                                         >  2:548693/1‑150 (MQ=255)
 atgTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCAc                                                                                                                                       >  1:805887/1‑151 (MQ=255)
      gtAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACg                                                                                                                                  <  1:562209/151‑1 (MQ=255)
            aCTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGc                                                                                                                            >  2:104913/1‑151 (MQ=255)
             cTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAgg                                                                                                                             <  1:829901/149‑1 (MQ=255)
               gCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCa                                                                                                                           <  1:980048/149‑1 (MQ=255)
                ccGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATg                                                                                                                        >  1:649250/1‑151 (MQ=255)
                   aTACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTa                                                                                                                      <  2:891651/150‑1 (MQ=255)
                    tACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTaa                                                                                                                     <  2:89121/150‑1 (MQ=255)
                           aTCAGTGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAgtgtg                                                                                                             <  2:399172/151‑1 (MQ=255)
                             cAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAgtgtgt                                                                                                            <  1:278871/150‑1 (MQ=255)
                             cAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGtt                                                                                                           >  2:957967/1‑151 (MQ=255)
                              aGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCAGCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGtt                                                                                                           <  2:389777/150‑1 (MQ=255)
                                 gcccgccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatc                                                                                                       >  2:777015/1‑151 (MQ=255)
                                  cccgccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatct                                                                                                      <  1:857634/151‑1 (MQ=255)
                                  cccgccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatct                                                                                                      <  1:356476/151‑1 (MQ=255)
                                  cccgccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatct                                                                                                      >  2:303484/1‑151 (MQ=255)
                                   ccgccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatctg                                                                                                     <  2:443647/151‑1 (MQ=255)
                                    cgccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatctga                                                                                                    <  2:945912/151‑1 (MQ=255)
                                      ccAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTtgatctgatc                                                                                                  >  1:175849/1‑151 (MQ=255)
                                       cagcagACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATcc                                                                                                 >  2:67413/1‑151 (MQ=255)
                                        agcagACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATcc                                                                                                 >  2:244363/1‑150 (MQ=255)
                                        agcagACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCa                                                                                                <  1:494967/151‑1 (MQ=255)
                                          cagACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAgg                                                                                              >  1:27453/1‑151 (MQ=255)
                                          cagACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAgg                                                                                              <  1:200541/151‑1 (MQ=255)
                                             aCCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTc                                                                                           <  2:585160/151‑1 (MQ=255)
                                                  agGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGttc                                                                                      >  2:338464/1‑150 (MQ=255)
                                                    gCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGtttt                                                                                     >  1:809155/1‑150 (MQ=255)
                                                    gCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGtttt                                                                                     >  2:899859/1‑150 (MQ=255)
                                                     cGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTAc                                                                                   >  2:1025275/1‑151 (MQ=255)
                                                      gAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTAcc                                                                                  <  2:911870/151‑1 (MQ=255)
                                                        aGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTAcccc                                                                                <  1:397095/151‑1 (MQ=255)
                                                         gATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCa                                                                               <  1:43235/151‑1 (MQ=255)
                                                          aTTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCAc                                                                              >  2:622326/1‑151 (MQ=255)
                                                           tttCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACt                                                                             >  2:388203/1‑151 (MQ=255)
                                                              cgtcgtAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTc                                                                            >  2:1023604/1‑149 (MQ=255)
                                                                tcgtAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAgg                                                                         >  2:977338/1‑150 (MQ=255)
                                                                 cgtAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGc                                                                       >  1:1013843/1‑151 (MQ=255)
                                                                           ataGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAAtgttg                                                             <  1:674705/151‑1 (MQ=255)
                                                                               agCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTac                                                          >  2:11186/1‑150 (MQ=255)
                                                                                  acaAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTAcacca                                                      <  1:440564/151‑1 (MQ=255)
                                                                                   caAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTAcaccac                                                     <  1:167029/151‑1 (MQ=255)
                                                                                     aGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACaccacc                                                    >  1:724373/1‑150 (MQ=255)
                                                                                       aTCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                 >  2:304411/1‑151 (MQ=255)
                                                                                           ggATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCa                                             <  1:853553/151‑1 (MQ=255)
                                                                                                 taacgacgAATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCAt                                       >  2:292991/3‑151 (MQ=255)
                                                                                                    cgacgGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccg                                    >  2:137823/1‑151 (MQ=255)
                                                                                                     gacgGATCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgc                                   <  1:726174/151‑1 (MQ=255)
                                                                                                          aTCTGCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAAt                              >  1:806796/1‑151 (MQ=255)
                                                                                                              gCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                          <  1:338032/151‑1 (MQ=255)
                                                                                                              gCTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                           <  1:444697/150‑1 (MQ=255)
                                                                                                               cTCCATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                          >  2:947840/1‑150 (MQ=255)
                                                                                                                 ccATCTGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                         <  2:406915/149‑1 (MQ=255)
                                                                                                                       gATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                 >  1:311752/1‑151 (MQ=255)
                                                                                                                        aTAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                 >  1:427800/1‑150 (MQ=255)
                                                                                                                        aTAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                 >  2:255522/1‑150 (MQ=255)
                                                                                                                        aTAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGa                <  2:1040793/151‑1 (MQ=255)
                                                                                                                          agagATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa              <  1:387908/151‑1 (MQ=255)
                                                                                                                              attgCCAGAGGCCTGTTTCAGCGCACATACTCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCACACTCAGGGCTGGAATGTTGTACACCACCATCGGCACACCATCCGCCGAATCAATAATTGCCCGATa              <  2:742905/144‑1 (MQ=255)
                                                                                                                                   cAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCa     <  1:700638/151‑1 (MQ=255)
                                                                                                                                      agGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtg  <  2:1064866/151‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
GATGTTGTAGGTACTGCCGATACCACCATCAGCGCCCGCCAGCAGACCAGAGGCGAAGATTTCGTCGTAACCGTTATAGAGCACAAGATCAGGATGTTCACGACGGATCTGCTCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTG  >  NC_000913/3372941‑3373225

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