Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,324,955 A→G 5.8% intergenic (‑44/+46) folP ← / ← ftsH 7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,324,9550AG 5.8% 199.7 / 2.7 69intergenic (‑44/+46)folP/ftsH7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base G (2/2);  total (38/31)
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.32e-01

TTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGT  >  NC_000913/3324808‑3325092
                                                                                                                                                   |                                                                                                                                         
ttATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCAc                                                                                                                                          >  4:193534/1‑149 (MQ=255)
ttATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCAc                                                                                                                                          >  3:131838/1‑149 (MQ=255)
      gTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCg                                                                                                                                    <  5:235641/149‑1 (MQ=255)
      gTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCg                                                                                                                                    <  3:45789/149‑1 (MQ=255)
      gTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCg                                                                                                                                    >  7:129603/1‑149 (MQ=255)
       tGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgg                                                                                                                                   <  1:5203/149‑1 (MQ=255)
             caTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGtttt                                                                                                                             >  8:134809/1‑149 (MQ=255)
             caTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGtttt                                                                                                                             <  3:111764/149‑1 (MQ=255)
                 cGAAAAGGAATCAGGCGGGACGGTGAGGAGCCCCATTACGTGACGATGGCTAAGGGCCAGGGAAGTACCCTGGGCAAAGCGTGTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttttc                                                                                                                          <  5:65482/148‑4 (MQ=255)
                  gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTa                                                                                                                        >  7:167016/1‑149 (MQ=255)
                    aaaGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTAc                                                                                                                       >  8:152656/1‑148 (MQ=255)
                      aGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGTCGTGCTCCGGGGttttttc                                                                                                                          >  6:34862/1‑140 (MQ=255)
                      aGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttttc                                                                                                                          <  5:34862/143‑4 (MQ=255)
                        gAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATa                                                                                                                  <  3:193534/149‑1 (MQ=255)
                          aTCAGGCGGCACGGTGAGCATCCCCATTCCGTGAGGATGGCTAAGGTCCAGTGACGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAACTGATAACTAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACa                                                                                                                <  2:186860/149‑1 (MQ=255)
                           tCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAg                                                                                                               <  7:152656/149‑1 (MQ=255)
                                cGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATc                                                                                                          <  5:59409/149‑1 (MQ=255)
                                cGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATc                                                                                                          <  8:129603/149‑1 (MQ=255)
                                 gTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCt                                                                                                         <  4:184713/149‑1 (MQ=255)
                                        tGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCgg                                                                                                  >  3:148401/1‑149 (MQ=255)
                                               ccccAAAACCTGAGGACGGCTAAGGTCCAGTGAACTACCCTGGGCAAAGAGTTTCATGCTGTAATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATCCAGTCATCTGATGCGGGAactta                                                                                           <  2:236499/149‑1 (MQ=255)
                                                 ccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATATGAAAAAACCCCGGAGCCGGCCCCGGGGTTTTCGGTACAATTACAGTCATCGGATGCGGGAacttact                                                                                         >  2:166369/1‑149 (MQ=255)
                                                 ccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttact                                                                                         >  4:9348/1‑149 (MQ=255)
                                                 ccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCACGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttact                                                                                         >  2:98110/1‑149 (MQ=255)
                                                  cATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttactt                                                                                        <  6:234271/149‑1 (MQ=255)
                                                        gTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGGCATCTGATGCGGGAACTTACTGGTCGcc                                                                                  >  2:195763/1‑149 (MQ=255)
                                                                  cTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATc                                                                                                          >  8:36804/1‑115 (MQ=255)
                                                                  cTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATc                                                                                                          <  7:36804/115‑1 (MQ=255)
                                                                  cTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctct                                                                        >  1:90118/1‑149 (MQ=255)
                                                                  cTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctct                                                                        >  7:197275/1‑149 (MQ=255)
                                                                  cTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctct                                                                        >  1:86775/1‑149 (MQ=255)
                                                                   tAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctct                                                                        >  4:55967/1‑148 (MQ=255)
                                                                         ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACACCGGAGCACGCCCCGGGGTTTTCGGTACTAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATgg                                                                 >  1:65774/1‑149 (MQ=255)
                                                                            gTGAAGTACCCTTGTCAAAGAGTTTCAGGATGTTATCCCGGGTATGAATTGATAAGAAAACACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtgt                                                              <  6:39886/149‑1 (MQ=255)
                                                                                 gTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTAccc                                                         >  7:101874/1‑149 (MQ=255)
                                                                                  tACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCg                                                        <  4:148401/149‑1 (MQ=255)
                                                                                     ccTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTTACATGGTGTTACCCGGGt                                                     >  6:135332/1‑149 (MQ=255)
                                                                                     ccTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGt                                                     >  1:19064/1‑149 (MQ=255)
                                                                                      cTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGtt                                                    >  5:119200/1‑149 (MQ=255)
                                                                                       tGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGCCATGGTGTTACCCGGGTTc                                                   >  5:150283/1‑149 (MQ=255)
                                                                                             aaGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATTCAGTCATCTGATGCGTGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTa                                             <  2:86775/149‑1 (MQ=255)
                                                                                                agTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTAcgc                                          <  8:167016/149‑1 (MQ=255)
                                                                                                 gTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATAAAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTAcgcg                                         >  4:173980/1‑149 (MQ=255)
                                                                                                  tttCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCgg                                        <  4:131838/149‑1 (MQ=255)
                                                                                                   ttCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTGACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCgg                                        <  7:134809/148‑1 (MQ=255)
                                                                                                       tgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcat                                   <  5:135332/149‑1 (MQ=255)
                                                                                                        gatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTaa                                                                               >  2:237029/1‑104 (MQ=255)
                                                                                                        gatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTaa                                                                               <  1:237029/104‑1 (MQ=255)
                                                                                                                  cTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGAc                          >  8:162039/1‑147 (MQ=255)
                                                                                                                  cTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGAc                          >  4:89974/1‑147 (MQ=255)
                                                                                                                         gAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGa                                                                      >  3:92863/1‑96 (MQ=255)
                                                                                                                         gAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGa                                                                      <  4:92863/96‑1 (MQ=255)
                                                                                                                            ttGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg              <  3:55967/149‑1 (MQ=255)
                                                                                                                             tGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg              >  3:127890/1‑148 (MQ=255)
                                                                                                                              gATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTgg             <  2:90118/148‑1 (MQ=255)
                                                                                                                               aTAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTAccc                                                         <  2:238094/103‑1 (MQ=255)
                                                                                                                               aTAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGa            <  2:19064/148‑1 (MQ=255)
                                                                                                                               aTAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGAc           >  6:40185/1‑149 (MQ=255)
                                                                                                                                tAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGAc           >  6:188111/1‑148 (MQ=255)
                                                                                                                                gtaccgaaaaCCCGGGGGCGTGCTCCGGGGGTGTCGGTACAAATACAGTCATCTGCTGCGGGAACTTACTTGTCGCCTTACTGGTCTGCCATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCTCCCGGGCGAGGAGGCTTTGGAc           >  7:249148/7‑148 (MQ=255)
                                                                                                                                   gAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGGGGAGCCTTTGGACTAc        >  1:180184/1‑148 (MQ=255)
                                                                                                                                    aaaaaaCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGa            <  2:52781/143‑1 (MQ=255)
                                                                                                                                    aaaaaaCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGa            >  1:52781/1‑143 (MQ=255)
                                                                                                                                    aaaaaaCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCAGCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCACCCGGACGAGGAGCCTTTGGACTACCa      >  1:162937/1‑149 (MQ=255)
                                                                                                                                        aaCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGt  >  7:164087/1‑149 (MQ=255)
                                                                                                                                           cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCggg                                                      >  4:15087/1‑94 (MQ=255)
                                                                                                                                           cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCggg                                                      <  3:15087/94‑1 (MQ=255)
                                                                                                                                                 gCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg              <  7:75474/128‑1 (MQ=255)
                                                                                                                                                 gCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTg              >  8:75474/1‑128 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
TTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGT  >  NC_000913/3324808‑3325092

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