Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,998,171 G→T 10.5% intergenic (‑71/‑145) PP_4405 ← / → PP_4406 sensory box protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,998,1710GT10.5% 216.0 / 13.5 76intergenic (‑71/‑145)PP_4405/PP_4406sensory box protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base T (4/4);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

GCTGATCAAATCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGAAACA  >  NC_002947/4998026‑4998316
                                                                                                                                                     |                                                                                                                                                 
gCTGATCAAATCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGc                                                                                                                                                <  2:278035/149‑1 (MQ=255)
       aaaTCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCTCGCGGCaatc                                                                                                                                         <  3:344536/149‑3 (MQ=255)
        aaTCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCTCTCGCGGCaatc                                                                                                                                         >  1:371800/1‑146 (MQ=255)
        aaTCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCTCTCGCGGCaatc                                                                                                                                         <  2:371800/148‑3 (MQ=255)
         aTCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAgat                                                                                                                                        >  7:287660/1‑148 (MQ=255)
            gAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGatgatg                                                                                                                                    >  6:155413/1‑149 (MQ=255)
             aaCAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGc                                                                                                                                   <  1:193132/149‑1 (MQ=255)
                     gtgCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACg                                                                                                                           >  8:256819/1‑149 (MQ=255)
                         aaTTGCTGTTCACTGAAGCCGGCTGGTACCGGCCTGGTGTCTTGACTCATGCTATGCTCTCCCTATTGCCCTGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGc                                                                                                                        <  1:132142/149‑1 (MQ=255)
                             gCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGc                                                                                                                    <  1:181168/148‑1 (MQ=255)
                                 ttCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCatat                                                                                                                <  5:410022/148‑1 (MQ=255)
                                 ttCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATg                                                                                                               >  2:367138/1‑149 (MQ=255)
                                  tCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCAG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCatatcc                                                                                                              >  1:64038/1‑147 (MQ=255)
                                   cACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGc                                                                                                              <  2:11952/148‑1 (MQ=255)
                                    aCTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCTCGCGGCaatc                                                                                                                                         <  1:76621/120‑3 (MQ=255)
                                    aCTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCTCGCGGCaatc                                                                                                                                         >  2:76621/1‑118 (MQ=255)
                                    aCTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGccc                                                                                                            >  6:29883/1‑149 (MQ=255)
                                            cGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGt                                                                                                    <  3:183613/149‑1 (MQ=255)
                                            cGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGt                                                                                                    <  2:64038/149‑1 (MQ=255)
                                                       ggCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGAGGAGCACAGTGGTGTTGATTACCGCGCGCGGCAAGATGATGCATATCACGGGAAAGCATATCCCCACCCACGTGAAAcca                                                                                             >  1:126904/1‑147 (MQ=255)
                                                             gtgtCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCaa                                                                                                <  4:365987/136‑1 (MQ=255)
                                                             gtgtCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCaa                                                                                                >  3:365987/1‑136 (MQ=255)
                                                             gtgtCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGc                                                                                    >  6:57956/1‑148 (MQ=255)
                                                                 cTTGACTCATGGACGCTCTCCCTTATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAAc                                                                                <  4:23224/148‑1 (MQ=255)
                                                                  ttGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCaa                                                                                                                      <  4:346542/109‑1 (MQ=255)
                                                                  ttGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCaa                                                                                                                      >  3:346542/1‑109 (MQ=255)
                                                                     aCTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGccc                                                                            <  5:57956/148‑1 (MQ=255)
                                                                     aCTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCt                                                                           <  5:154110/149‑1 (MQ=255)
                                                                         aTGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCTAACGCCCTGAtt                                                                       >  1:5017/1‑149 (MQ=255)
                                                                         aTGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCa                                                                                  <  3:258538/138‑1 (MQ=255)
                                                                         aTGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCa                                                                                  >  4:258538/1‑138 (MQ=255)
                                                                         aTGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGa                                                                         <  5:29883/147‑1 (MQ=255)
                                                                             gACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTc                                                                     <  7:201388/148‑1 (MQ=255)
                                                                              aCGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACgg                                                                                                                          >  2:50302/1‑94 (MQ=255)
                                                                               cGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGAc                                                                  <  8:287660/149‑1 (MQ=255)
                                                                                      ccTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGcaca                                                           >  1:396209/1‑149 (MQ=255)
                                                                                      ccTGATGCCCAGCAACTGCGCGCCGATCTCGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGcac                                                            <  6:349934/148‑1 (MQ=255)
                                                                                        tGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAgc                                                         >  1:383460/1‑149 (MQ=255)
                                                                                              ccAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGTAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTgg                                                    >  6:225509/1‑148 (MQ=255)
                                                                                              ccAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTgg                                                    >  8:111111/1‑148 (MQ=255)
                                                                                               cAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCa                                                  >  4:21621/1‑149 (MQ=255)
                                                                                                aGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGt                                                                                      <  3:54245/112‑1 (MQ=255)
                                                                                                aGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGt                                                                                      >  4:54245/1‑112 (MQ=255)
                                                                                                     cTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGAtt                                            <  2:383460/149‑1 (MQ=255)
                                                                                                      tGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGAttt                                           >  8:185992/1‑149 (MQ=255)
                                                                                                      tGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACACACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGAttt                                           >  4:14905/1‑149 (MQ=255)
                                                                                                           tCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTg                                                     <  7:44537/134‑1 (MQ=255)
                                                                                                           tCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTg                                                     >  8:44537/1‑134 (MQ=255)
                                                                                                           tCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTgcgc                                      >  7:21306/1‑149 (MQ=255)
                                                                                                             cGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCtt                                    >  5:97402/1‑149 (MQ=255)
                                                                                                             cGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCt                                     <  1:367138/148‑1 (MQ=255)
                                                                                                              gATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTg                                   <  5:155413/149‑1 (MQ=255)
                                                                                                              gATCAAGGCGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCTGCACAGCGCTGGCATGGATTTTGCGCTTg                                   >  8:215931/1‑149 (MQ=255)
                                                                                                                  aaGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTa                                <  3:386605/148‑1 (MQ=255)
                                                                                                                     gggCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTg                            <  6:379496/149‑1 (MQ=255)
                                                                                                                       gCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGcc                          <  3:14905/149‑1 (MQ=255)
                                                                                                                            cG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGTAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACga                     <  5:225509/149‑1 (MQ=255)
                                                                                                                            cG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTg                                   >  8:298526/1‑135 (MQ=255)
                                                                                                                            cG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTg                                   <  7:298526/135‑1 (MQ=255)
                                                                                                                            cG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACga                     <  7:256819/149‑1 (MQ=255)
                                                                                                                                 aCAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAAc                       >  1:64778/1‑145 (MQ=255)
                                                                                                                                     tggtgTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGAc                                                                  <  4:279583/98‑1 (MQ=255)
                                                                                                                                     tggtgTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGAc                                                                  >  3:279583/1‑98 (MQ=255)
                                                                                                                                     tggtgTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCGATCGAGATTTc               >  7:271818/1‑149 (MQ=255)
                                                                                                                                     tggtgTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTc               >  3:245488/1‑149 (MQ=255)
                                                                                                                                        atcatcttGCCGCTCGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTcc              >  7:111050/7‑147 (MQ=255)
                                                                                                                                         tcatcttGCCGCTCGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAag           <  8:111050/144‑1 (MQ=255)
                                                                                                                                          ttGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGGTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAagg          <  1:288834/149‑1 (MQ=255)
                                                                                                                                          ttGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAagg          <  6:97402/149‑1 (MQ=255)
                                                                                                                                          ttGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAag           >  4:359278/1‑148 (MQ=255)
                                                                                                                                            tcttGCCGCTCGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAaggag        >  2:86399/3‑149 (MQ=255)
                                                                                                                                             aTTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAaggagg       <  1:169920/149‑1 (MQ=255)
                                                                                                                                              ttGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAaggagg       <  8:21306/148‑1 (MQ=255)
                                                                                                                                                 ccGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGaaa    <  3:21621/148‑1 (MQ=255)
                                                                                                                                                 ccGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGaaa    >  8:151754/1‑148 (MQ=255)
                                                                                                                                                 ccGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGaaa    >  8:292037/1‑148 (MQ=255)
                                                                                                                                                  cgcgAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGAAACa  <  7:215931/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCTGATCAAATCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCG‑‑‑ACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGAAACA  >  NC_002947/4998026‑4998316

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