Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,457 G→A 13.5% intergenic (+22/+49) hutF → / ← PP_5037 formiminoglutamate deiminase/lipocalin family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,740,4570GA13.5% 137.4 / 10.3 59intergenic (+22/+49)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (21/30);  new base A (2/6);  total (23/36)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

GTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTC  >  NC_002947/5740311‑5740604
                                                                                                                                                    |                                                                                                                                                   
gTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTTGGTTGGAGcc                                                                                                                                                   >  4:360534/1‑149 (MQ=255)
    cGATCGCCAGGTGCGCGATGTGATTGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTc                                                                                                                                               <  3:110302/149‑1 (MQ=255)
    cGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTc                                                                                                                                               <  7:120223/149‑1 (MQ=255)
    cGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTc                                                                                                                                               <  7:125903/149‑1 (MQ=255)
        cGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTcgcgg                                                                                                                                           >  2:195893/1‑148 (MQ=255)
         gCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGt                                                                                                                                           <  3:72266/148‑1 (MQ=255)
          ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTg                                                                                                                                          >  5:374808/1‑148 (MQ=255)
          ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGc                                                                                                                                         <  4:7895/149‑1 (MQ=255)
           cAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGct                                                                                                                                        <  1:195893/149‑2 (MQ=255)
             ggTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGcttg                                                                                                                                      <  2:207965/149‑4 (MQ=255)
               tGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGccc                                                                                                                                                  <  1:13380/135‑1 (MQ=255)
               tGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGccc                                                                                                                                                  >  2:13380/1‑135 (MQ=255)
                     aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGaa                                                                                                                              <  2:257603/149‑1 (MQ=255)
                      tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                             <  2:100575/149‑1 (MQ=255)
                       gtgAGGGTGACTGGGCACTGGGTGGTTCGCCAGGGGCGGCATGCCCAGGAAGAGGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTTAGCGTACCGGGGTGGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                             <  2:247755/148‑1 (MQ=37)
                        tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                           <  3:54414/149‑1 (MQ=255)
                        tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg                                                                                                                            <  2:244165/148‑1 (MQ=255)
                          aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTTGGGTTGAGCCCCTAGCCTGGAAGCCCGCCAAAGGtt                                                                                                                         >  2:356442/1‑149 (MQ=255)
                          aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggc                                                                                                                          <  2:191283/148‑2 (MQ=38)
                          aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggc                                                                                                                          >  1:191283/1‑147 (MQ=38)
                           ttttGAATGGGGAATGGGTGTTGCGTCAGGGGCGGCATGCGCAGGTAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGTCTGAGCGTACCGGGGTTGCGTTGCAACGCTTCGCGGGGTTGCACGCGAAGGGGTg                                                                                                                        <  8:373868/146‑1 (MQ=255)
                           tGGTGAATGGGCAATGGGGGGTGCGCCAGGGGCGGCAAGCCCATGAAGTTGTAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTTCTGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                         <  1:304224/149‑1 (MQ=255)
                            ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTTCGTTTGAGCCCATACCCTGGAAGCCCCGGAAAGGTTg                                                                                                                        >  8:261094/1‑148 (MQ=255)
                            ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                           >  2:349803/1‑145 (MQ=255)
                            ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                           <  1:349803/145‑1 (MQ=255)
                              tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTTGAGCCCCTGGCGGGGCAGCCCGCGCAGGGGTggca                                                                                                                     >  5:345102/1‑146 (MQ=255)
                              tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                              >  7:356195/1‑140 (MQ=255)
                              tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                              <  8:356195/140‑1 (MQ=255)
                                aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCATAGCGTGGAAGCCCGCGAAGGGTTtgaacg                                                                                                                   >  2:21489/1‑149 (MQ=255)
                                aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCCTTGCAGCCCATAGCGTGCCACCCCGCGCAGGGTTgccacg                                                                                                                   >  6:173765/1‑149 (MQ=255)
                                            gtggggCGCCAGGGGCGGCATGCCCAGGAAGAGGAACGCGGGCGGGCGTTTGCCAAGGTGGTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTGGCAGCCCTGCGCGTGCAACCCCGCGAAGGGTTGCAACGCAACCCCaaaaa                                                                                                       <  7:93242/149‑1 (MQ=255)
                                             tcgtccgcCAGGGGCGGCATGCCCAGGAAGATGAACGCGGCCGGGCGTTTGCCAAGGTGTTGCCCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCaaaaaa                                                                                                      <  5:198105/144‑1 (MQ=255)
                                                   gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCATAGCGTGGAAGGCCGCGCAGGGGTGCAACGCCACCCCACAAAACTCAGt                                                                                                >  6:360659/1‑149 (MQ=255)
                                                   gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                             >  1:135291/1‑120 (MQ=38)
                                                   gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                             <  2:135291/120‑1 (MQ=38)
                                                    ccAGGGGCCGCATGCCCACGAAGATGAAAGCGGGCGGGCGTTTGCCAAGCTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTAACCGGGCAACGCCCCCCAGGGTTGCAAAGCCACCCCAAAAAACtagctg                                                                                               >  8:283920/1‑144 (MQ=255)
                                                             gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCCATTCAGCCCATACCGTGGCACCCCCCCAAGGGTTGCAACCCAACCCCAAAAAACAGAGAGCAtttttt                                                                                       >  6:183937/1‑147 (MQ=255)
                                                                  cccAGGAAGATGAAAGGGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAtt                                                                                  <  7:261094/148‑1 (MQ=255)
                                                                  cccAGGAAGATGAAAGCGGGCGGTCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCTCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                 <  5:183937/149‑1 (MQ=255)
                                                                  cccAGGAAGATGAAAGCGGGCGGGCTTTTGCCAAGGTGTTGCGCCTGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                 <  5:173765/149‑1 (MQ=255)
                                                                  cccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                 <  1:21489/149‑1 (MQ=255)
                                                                   ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                                <  6:377057/149‑1 (MQ=255)
                                                                    cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCCGCCCCTCGCCGGGCAGCCCCCGCCCGGGTGCCCCGCCACCCCAAAAAACACAGTGCCATTTTAATAttttc                                                                               >  7:4128/1‑147 (MQ=255)
                                                                    cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                               <  5:360659/149‑1 (MQ=255)
                                                                      ggAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTg                                                                             <  6:374808/149‑1 (MQ=255)
                                                                       taaGATCAAAGCGGGGGGGCGTTTGCCATGGTGTTGCGCCTGGTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGa                                                                            <  6:131083/148‑1 (MQ=255)
                                                                              aaaGCGGGCGGGCGTTTGCAAAGGTGTTGCGCCAGCTTTTGGGCTTAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                      <  1:8980/148‑1 (MQ=255)
                                                                                    ggcgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCAAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgcc                                                               <  1:356442/149‑1 (MQ=255)
                                                                                    ggcgggcgTTTGCCAAGGTGTGGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGGCCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgcc                                                               <  3:360534/149‑1 (MQ=255)
                                                                                                     gtgtTGCGCCAGCTTGTGGGCTGCGCGTACCGGGGTTGGGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                               <  8:4128/148‑1 (MQ=255)
                                                                                                          gcgcCAGCTTTTGGGCTGCGCGTACCGGGGTGGCGTGGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCAt                                          <  7:283920/148‑1 (MQ=255)
                                                                                                           cgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                        <  3:40293/149‑1 (MQ=255)
                                                                                                                                        gttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCACTGGACAGCAGGt            >  2:34171/1‑149 (MQ=255)
                                                                                                                                            gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtct         >  3:362277/1‑149 (MQ=255)
                                                                                                                                              gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcg                             >  7:360020/1‑127 (MQ=255)
                                                                                                                                              gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcg                             <  8:360020/127‑1 (MQ=255)
                                                                                                                                                  cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt   >  1:224756/1‑149 (MQ=255)
                                                                                                                                                  cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt   >  7:283667/1‑149 (MQ=255)
                                                                                                                                                   aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc  >  3:63730/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTC  >  NC_002947/5740311‑5740604

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: