Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 268,746 G→C 8.6% intergenic (‑77/+23) PP_0216 ← / ← desA sensory box/GGDEF family protein/delta‑9 fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947268,7460GC 8.6% 174.9 / 7.7 58intergenic (‑77/+23)PP_0216/desAsensory box/GGDEF family protein/delta‑9 fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base G (23/30);  new base C (3/2);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGTTG  >  NC_002947/268599‑268895
                                                                                                                                                   |                                                                                                                                                     
gggAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgggg                                                                                                                                                     <  1:541959/150‑1 (MQ=255)
gggAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgggg                                                                                                                                                     >  1:444523/1‑150 (MQ=255)
 ggAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTTATCTTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggggg                                                                                                                                                    <  1:1152540/150‑1 (MQ=255)
      aTCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCCCAGCGGGGCGTg                                                                                                                                               <  1:1707230/150‑1 (MQ=255)
       tCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGc                                                                                                                                              <  1:951918/150‑1 (MQ=255)
       tCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGc                                                                                                                                              <  1:952119/150‑1 (MQ=255)
         gTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCg                                                                                                                                            >  1:146265/1‑150 (MQ=255)
               aTGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAc                                                                                                                                      <  1:469804/150‑1 (MQ=255)
                   tCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGCAACCTTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcg                                                                                                                                  <  1:1841578/150‑1 (MQ=255)
                   tCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcg                                                                                                                                  <  2:717238/150‑1 (MQ=255)
                   tCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGACGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcg                                                                                                                                  <  2:1001378/150‑1 (MQ=255)
                      tGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTc                                                                                                                               <  2:1673977/150‑1 (MQ=255)
                        gCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAg                                                                                                                             >  2:623506/1‑150 (MQ=255)
                         ccATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAgg                                                                                                                            <  2:902552/150‑1 (MQ=255)
                             gtgCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCTGTCAGGCGcc                                                                                                                        <  1:1046500/150‑1 (MQ=255)
                                    gAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCggg                                                                                                                 >  2:424355/1‑150 (MQ=255)
                                     aGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGgc                                                                                                                <  1:1601029/150‑1 (MQ=255)
                                            ctcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGa                                                                                                         <  1:323613/150‑1 (MQ=255)
                                             tcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCCGACCGgcgcctgacc                                                                                                        <  2:73226/150‑7 (MQ=255)
                                               aaCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTa                                                                                                      >  1:433916/1‑150 (MQ=255)
                                                     aaaGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCtt                                                                                                <  1:1473312/150‑1 (MQ=255)
                                                     aaaGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATGCCGGATTCGAGTGGGGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCtt                                                                                                <  1:415283/150‑1 (MQ=255)
                                                       aGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                              >  1:420579/1‑150 (MQ=255)
                                                       aGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                              <  2:420579/150‑1 (MQ=255)
                                                          gTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACTGCGGTCAGGCGCCGGTCGGGCGCATCGATTAAGTCTTGAGtt                                                                                           <  2:1261447/150‑1 (MQ=255)
                                                          gTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCg                                                                                                                       >  1:966421/1‑122 (MQ=255)
                                                          gTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCg                                                                                                                       <  2:966421/122‑1 (MQ=255)
                                                               gCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgct                                                                                      >  1:1880403/1‑150 (MQ=255)
                                                                 cAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgctgc                                                                                    >  2:525623/1‑150 (MQ=255)
                                                                   gCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCCACGGCTGGCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGGTGTGCTGCgg                                                                                  >  1:1058461/1‑149 (MQ=255)
                                                                     aTCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAg                                                                                <  2:433916/150‑1 (MQ=255)
                                                                         ggCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTcgc                                                                            >  2:1474269/1‑150 (MQ=255)
                                                                         ggCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGCTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTcgc                                                                            >  1:1896224/1‑150 (MQ=255)
                                                                            ccGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGTCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCa                                                                         >  2:257966/1‑150 (MQ=255)
                                                                              gCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGt                                                                       <  1:1474269/150‑1 (MQ=255)
                                                                                cACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGc                                                                     >  1:690625/1‑150 (MQ=255)
                                                                                cACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGc                                                                     <  1:1853386/150‑1 (MQ=255)
                                                                                     aGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGa                                                                <  1:525623/150‑1 (MQ=255)
                                                                                          atatTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCa                                                                                                             <  1:761004/100‑1 (MQ=255)
                                                                                          atatTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCa                                                                                                             >  2:761004/1‑100 (MQ=255)
                                                                                              tCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTgg                                                       <  1:416636/150‑1 (MQ=255)
                                                                                                   aTTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCg                                                  <  1:1024823/150‑1 (MQ=255)
                                                                                                     tCGAGTGGAGAGGTGGGGGCGAAGTTCCGTATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGCGTTGTGCTGCGAAGTCGCGCAGTGTCTGGATGCCGCTGGTTTCGGc                                                <  2:676789/150‑1 (MQ=255)
                                                                                                     tCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGc                                                <  2:146265/150‑1 (MQ=255)
                                                                                                       gAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCt                                              >  1:1826099/1‑150 (MQ=255)
                                                                                                          tGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCAt                                           <  1:623506/150‑1 (MQ=255)
                                                                                                               agGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTAcc                                      >  1:945510/1‑150 (MQ=255)
                                                                                                                    ggggCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCTTGGCTGCCGCTGGTTTCGGCCTCATGTACCCAGTc                                 <  1:1026858/150‑1 (MQ=255)
                                                                                                                    ggggCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTc                                 >  1:830896/1‑150 (MQ=255)
                                                                                                                            gTTCCGAATGGCACGCACCGCTGGGGGGCGGGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGGCTAATTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATgg                         <  1:1772536/150‑1 (MQ=255)
                                                                                                                             ttCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGc                        >  1:738873/1‑150 (MQ=255)
                                                                                                                               ttgatcggCACGCCCCCCAGCGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTg                      >  2:480231/7‑150 (MQ=255)
                                                                                                                                   tcggCACGCCCCCCAGCGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCttt                              <  2:1509510/136‑1 (MQ=255)
                                                                                                                                   tcggCACGCCCCCCAGCGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCttt                              >  1:1509510/3‑138 (MQ=255)
                                                                                                                                    tGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCTGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGc                 <  1:840949/150‑1 (MQ=255)
                                                                                                                                        aCGCCCCCCAGCGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGCTGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCCGCATGt             >  2:86858/1‑150 (MQ=255)
                                                                                                                                         cGCCCCGCTGGGGGTTGTTCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTAGAGTTGTGCTGCGAAGTCGCACATTGCTTGGATGTCGCTGGTTTCGGCCTCATGTACCCAGTCTTTTATTGCTGCCAGCATGTc            >  2:594511/1‑150 (MQ=255)
                                                                                                                                            cccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTg         >  1:1915158/1‑150 (MQ=255)
                                                                                                                                             ccGCTGGGGGGCGTGTCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTgg        >  1:543774/1‑150 (MQ=255)
                                                                                                                                             ccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGGGTAAGTCGTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTgg        >  1:69338/1‑150 (MQ=255)
                                                                                                                                              cGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGCCGCGACGACCATGCCAGCTGTCTTTTTCTGTCTCATGTACCCAGTCTTTCATGGCTTCCAGCATGTCGTGGc       >  1:223789/1‑150 (MQ=255)
                                                                                                                                                   gggggCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGTTg  <  2:1880403/150‑1 (MQ=255)
                                                                                                                                                   gggggCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGTTg  >  1:70073/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGTTG  >  NC_002947/268599‑268895

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