Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 268,744 T→A 8.2% intergenic (‑75/+25) 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,7440TA 8.2% 181.9 / 5.8 61intergenic (‑75/+25)PP_0216/desAsensory box/GGDEF family protein/delta‑9 fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base A (3/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGC  >  NC_002947/268599‑268890
                                                                                                                                                 |                                                                                                                                                  
gggAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgggg                                                                                                                                                >  1:444523/1‑150 (MQ=255)
gggAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgggg                                                                                                                                                <  1:541959/150‑1 (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                                                                                         <  2:420579/150‑1 (MQ=255)
                                                       aGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                         >  1:420579/1‑150 (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                                                                                                        >  2:761004/1‑100 (MQ=255)
                                                                                          atatTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCa                                                                                                        <  1:761004/100‑1 (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)
                                                                                                                                                 |                                                                                                                                                  
GGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGC  >  NC_002947/268599‑268890

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: ‑