Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,483,939 G→A 6.8% V269I (GTC→ATC)  puuA‑I → glutamate‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,483,9390GA 6.8% 180.3 / 2.3 61V269I (GTC→ATC) puuA‑Iglutamate‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base A (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

GCACAGCGACGATGCCCTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACG  >  NC_002947/2483794‑2484078
                                                                                                                                                     |                                                                                                                                           
gCACAGCGACGATGCCCTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTc                                                                                                                                           <  4:164308/148‑1 (MQ=255)
  acaGCGACGATGCCCTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGc                                                                                                                                        <  6:182671/149‑1 (MQ=255)
      cgacgaTGCCCTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTg                                                                                                                                     >  3:131273/1‑148 (MQ=255)
          gaTGCCCTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTacga                                                                                                                                >  1:90007/1‑149 (MQ=255)
             gCCCTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACATCAGCCTGTACGAcgc                                                                                                                             >  3:249142/1‑149 (MQ=255)
              cccTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCa                                                                                                                                          <  1:229938/135‑1 (MQ=255)
              cccTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCa                                                                                                                                          >  2:229938/1‑135 (MQ=255)
               ccTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGAcgccg                                                                                                                           <  5:369978/149‑1 (MQ=255)
               ccTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACAAcgccg                                                                                                                           >  7:220443/1‑149 (MQ=255)
                cTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGAcgccgc                                                                                                                          >  4:237107/1‑149 (MQ=255)
                    gTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCCAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCgg                                                                                                                       >  7:267384/1‑148 (MQ=255)
                       cggcggACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACATCAGCCTGTACGACGCCGCCGGCAAt                                                                                                                   <  4:249142/149‑1 (MQ=255)
                        ggcggACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATa                                                                                                                  <  2:90007/149‑1 (MQ=255)
                              cTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATcc                                                                                                            <  4:23970/149‑1 (MQ=255)
                                 gTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGCCGACGCCGCCAATAACATCCTgg                                                                                                         >  6:289449/1‑149 (MQ=255)
                                    ggCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCgg                                                                                                      >  3:71110/1‑149 (MQ=255)
                                       aCTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCggg                                                                                                     >  2:86931/1‑147 (MQ=255)
                                       aCTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCggg                                                                                                     >  4:377740/1‑147 (MQ=255)
                                       aCTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGaa                                                                                                   <  5:314210/149‑1 (MQ=255)
                                           ttCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCg                                                                                                <  3:377740/148‑1 (MQ=255)
                                                  gCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACATCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAgcg                                                                                         >  4:30280/1‑148 (MQ=255)
                                                     gTACCCGCGGTATCGCGTTGAAGCACCGCTACCGTGCCAGCTTCATGAGCAAACCCTACTTGGACGCGCCAGGCACTTAGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAgcg                                                                                         <  1:126556/149‑1 (MQ=255)
                                                        cccGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCATGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTg                                                                                  >  3:339580/1‑149 (MQ=255)
                                                           gcgGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATc                                                                                                             >  3:67027/1‑119 (MQ=255)
                                                           gcgGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATc                                                                                                             <  4:67027/119‑1 (MQ=255)
                                                           gcgGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTgcgc                                                                               <  3:39137/149‑1 (MQ=255)
                                                           gcgGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTgcgc                                                                               <  4:71110/149‑1 (MQ=255)
                                                           gcgGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACAACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTgcgc                                                                               <  8:220443/149‑1 (MQ=255)
                                                                    cgATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTg                                                                      <  1:86931/149‑1 (MQ=255)
                                                                     gATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTgg                                                                     >  4:396731/1‑149 (MQ=255)
                                                                             aCGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGc                                                             <  1:243446/149‑1 (MQ=255)
                                                                             aCGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACCTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGACCAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGc                                                             <  4:131273/149‑1 (MQ=255)
                                                                              cGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCggg                                                                                                     >  6:112236/1‑108 (MQ=255)
                                                                              cGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCggg                                                                                                     <  5:112236/108‑1 (MQ=255)
                                                                               ggCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCt                                                           >  1:155645/1‑149 (MQ=255)
                                                                                    ccGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAg                                                                                           >  5:140890/1‑112 (MQ=255)
                                                                                    ccGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAg                                                                                           <  6:140890/112‑1 (MQ=255)
                                                                                    ccGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGaa                                                       <  5:289449/148‑1 (MQ=255)
                                                                                          cAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACtgctgc                                                <  3:237107/149‑1 (MQ=255)
                                                                                          cAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACtgctgc                                                >  3:40032/1‑149 (MQ=255)
                                                                                            gCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTa                                                                                                                                   >  6:69014/1‑64 (MQ=255)
                                                                                                 aTGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCATGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACt                                          <  4:339580/148‑1 (MQ=255)
                                                                                                       aaaCCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTgg                                                         >  6:396442/1‑127 (MQ=255)
                                                                                                       aaaCCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTgg                                                         <  5:396442/127‑1 (MQ=255)
                                                                                                               ccTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCttt                            <  4:281779/148‑1 (MQ=255)
                                                                                                               ccTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCttt                            <  3:396731/148‑1 (MQ=255)
                                                                                                               ccTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTg                           >  8:117173/1‑149 (MQ=255)
                                                                                                               ccTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTg                           >  4:4644/1‑149 (MQ=255)
                                                                                                                cTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGc                          <  2:155645/149‑1 (MQ=255)
                                                                                                                cTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACATCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGc                          <  3:30280/149‑1 (MQ=255)
                                                                                                                       cgcCAGGC‑‑‑‑AGCTGCATGCATGTGAAAGTCAGCCGGTACGACGCCTCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAAc                   <  2:80910/149‑1 (MQ=255)
                                                                                                                         ccAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACca                 <  8:267384/149‑1 (MQ=255)
                                                                                                                           aGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACcac                >  4:171715/1‑148 (MQ=255)
                                                                                                                                       gcatgcatGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGAGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTAc       >  8:111968/1‑149 (MQ=255)
                                                                                                                                            catgTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACg  >  6:268727/1‑149 (MQ=255)
                                                                                                                                             atgTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGccccc                      >  8:6325/1‑128 (MQ=255)
                                                                                                                                             atgTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGccccc                      <  7:6325/128‑1 (MQ=255)
                                                                                                                                               gtgCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCttt                            >  4:39028/1‑120 (MQ=255)
                                                                                                                                               gtgCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCttt                            <  3:39028/120‑1 (MQ=255)
                                                                                                                                                tgCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTgg                                                                     <  6:216612/78‑1 (MQ=255)
                                                                                                                                                tgCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTgg                                                                     >  5:216612/1‑78 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GCACAGCGACGATGCCCTCAGTGCGGCGGACTGGTCGGCACTGTTCTGCCGCAGTACCCGCGGTATCGCGATGAAGCACGGCTACCGTGCCAGCTTCATGAGCAAACCCTACCTGGACGCGCCAGGC‑‑‑‑AGCGGCATGCATGTGCACGTCAGCCTGTACGACGCCGCCGGCAATAACATCCTGGCGGGAAGCGAACAGCGCAAGCTGCGCCATGCGGTGGCAGGCTGCCTGGAACTGCTGCCGCACTGCATGCCGATCTTTGCCCCCAACCACAATGCCTACCGACG  >  NC_002947/2483794‑2484078

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