Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,444,629 G→A 5.6% A87A (GCC→GCT sulA ← cell division inhibitor SulA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,444,6290GA 5.6% 271.7 / 3.0 88A87A (GCC→GCTsulAcell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base G (37/46);  new base A (3/2);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

GCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCGGGCTGG  >  NC_002947/2444483‑2444776
                                                                                                                                                   |                                                                                                                                                    
gCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGc                                                                                                                                                    >  3:73826/1‑149 (MQ=255)
    cATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTgg                                                                                                                                                <  8:310370/149‑1 (MQ=255)
    cATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTgg                                                                                                                                                <  3:145550/149‑1 (MQ=255)
    cATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTgg                                                                                                                                                <  1:263666/149‑1 (MQ=255)
     aTTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTggg                                                                                                                                              >  5:65596/1‑149 (MQ=255)
     aTTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTggg                                                                                                                                              >  8:29158/1‑149 (MQ=255)
      ttGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGt                                                                                                                                             <  3:243286/149‑1 (MQ=255)
         aCGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACATGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCa                                                                                                                                           <  6:224975/148‑1 (MQ=255)
         aCGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAg                                                                                                                                          >  8:254043/1‑149 (MQ=255)
              gCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTg                                                                                                                                     >  3:142769/1‑149 (MQ=255)
                ccagccagCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGc                                                                                                                                   >  5:66868/1‑149 (MQ=255)
                ccagccagCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGc                                                                                                                                   >  7:240188/1‑149 (MQ=255)
                ccagccagCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAAGCCTGACCCAGGCCTGGc                                                                                                                                   <  5:215511/149‑1 (MQ=255)
                 cagccagcTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCg                                                                                                                                  <  2:383543/149‑1 (MQ=255)
                  agccagcTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCgg                                                                                                                                 >  1:2786/1‑149 (MQ=255)
                   gccagcTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGt                                                                                                                                <  1:95947/149‑1 (MQ=255)
                     cagcTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTgg                                                                                                                              <  7:254043/149‑1 (MQ=255)
                       gcTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTgggg                                                                                                                            <  2:2786/149‑1 (MQ=255)
                         tCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGGGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGc                                                                                                                           >  6:382629/1‑148 (MQ=255)
                         tCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGc                                                                                                                           >  7:263620/1‑148 (MQ=255)
                             cacGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGAt                                                                                                                        >  8:18374/1‑147 (MQ=255)
                             cacGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGAt                                                                                                                        >  7:344937/1‑147 (MQ=255)
                             cacGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGAt                                                                                                                        >  2:367106/1‑147 (MQ=255)
                             cacGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGAt                                                                                                                        >  5:58146/1‑147 (MQ=255)
                               cGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAgg                                                                                                                    <  2:129651/149‑1 (MQ=255)
                                   gtgGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGGTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTca                                                                                                                <  3:244581/149‑1 (MQ=255)
                                      gCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTcagc                                                                                                              <  3:388490/148‑1 (MQ=255)
                                       cTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTcagc                                                                                                              >  1:206020/1‑147 (MQ=255)
                                        tGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCtt                                                                                                                                      >  6:68893/1‑122 (MQ=255)
                                             cccAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGc                                                                                                      <  4:73826/149‑1 (MQ=255)
                                              ccAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGcc                                                                                                     >  8:362267/1‑149 (MQ=255)
                                              ccAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGcc                                                                                                     <  6:58146/149‑1 (MQ=255)
                                               cAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAAGCCTGACCCAGGCCTgg                                                                                                                                    <  2:381089/117‑1 (MQ=255)
                                               cAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAAGCCTGACCCAGGCCTgg                                                                                                                                    >  1:381089/1‑117 (MQ=255)
                                                aGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGTGTCAGCCAGCGGGCCTc                                                                                                   <  6:184850/149‑1 (MQ=255)
                                                   tGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGTCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTc                                                                                                <  1:96186/149‑1 (MQ=255)
                                                   tGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTc                                                                                                <  1:101271/149‑1 (MQ=255)
                                                    gCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGAATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTc                                                                                                >  1:257468/1‑149 (MQ=255)
                                                     cAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTcc                                                                                               >  7:358979/1‑148 (MQ=255)
                                                     cAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCct                                                                                              >  8:145510/1‑149 (MQ=255)
                                                      aaaGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCctc                                                                                             >  6:111967/1‑149 (MQ=255)
                                                      aaaGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCctc                                                                                             <  8:344937/149‑1 (MQ=255)
                                                        aGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCtt                                                                                           >  4:107026/1‑149 (MQ=255)
                                                        aGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCtt                                                                                           <  7:318870/149‑1 (MQ=255)
                                                        aGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCtt                                                                                           <  4:242774/149‑1 (MQ=255)
                                                                aGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGt                                                                                   <  2:206020/149‑1 (MQ=255)
                                                                aGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGt                                                                                   >  1:249364/1‑149 (MQ=255)
                                                                  gCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTc                                                                                 >  6:229420/1‑149 (MQ=255)
                                                                   ccAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTcg                                                                                <  7:362267/149‑1 (MQ=255)
                                                                    cAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTcgc                                                                               >  1:355946/1‑149 (MQ=255)
                                                                      gTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTGCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCgcgc                                                                             <  3:287998/149‑1 (MQ=255)
                                                                          cAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCgcagca                                                                         >  6:31943/1‑149 (MQ=255)
                                                                               ggcTTTGATTGCTGCGGGTTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACgggg                                                                    <  7:56433/148‑1 (MQ=255)
                                                                                 cTTTGATTGCTGCGGGGTTGGAGCAGCAGAATAGGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGGGGTGCGTTCAGGGTCAGCCAGCGGGCCTCGTCCTCTTGGGTGAGTTCGCGCAGCACgtgtgc                                                                  <  4:262933/149‑6 (MQ=255)
                                                                                 cTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGc                                                                  <  6:65596/149‑1 (MQ=255)
                                                                                  tttGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGc                                                                  >  5:145225/1‑148 (MQ=255)
                                                                                       aTTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCcagca                                                             <  2:257468/149‑1 (MQ=255)
                                                                                           ctgcGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGc                                                          <  8:263620/148‑1 (MQ=255)
                                                                                            tgcGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCtt                                                        <  7:145510/149‑1 (MQ=255)
                                                                                                gggTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGc                                                    <  2:355946/149‑1 (MQ=255)
                                                                                                   ttGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGt                                                 <  6:211710/149‑1 (MQ=255)
                                                                                                   ttGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGt                                                 >  3:366724/1‑149 (MQ=255)
                                                                                                       agcagcagAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGc                                              >  6:284168/1‑148 (MQ=255)
                                                                                                         cagcagAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAg                                           >  3:188811/1‑149 (MQ=255)
                                                                                                              gAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGgcg                                      >  4:95212/1‑149 (MQ=255)
                                                                                                                aTACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGgcgc                                     <  8:358979/148‑1 (MQ=255)
                                                                                                                  aCGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGcgccgc                                  >  2:85469/1‑149 (MQ=255)
                                                                                                                       cACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGTAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAg                              <  1:85469/148‑1 (MQ=255)
                                                                                                                         cGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGa                           <  7:29158/149‑1 (MQ=255)
                                                                                                                         cGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGa                           <  4:230522/149‑1 (MQ=255)
                                                                                                                          ggTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGAc                          <  1:266434/149‑1 (MQ=255)
                                                                                                                           gTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACa                         <  7:18374/149‑1 (MQ=255)
                                                                                                                           gTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACa                         <  2:249364/149‑1 (MQ=255)
                                                                                                                                gCCCGGCATCGCGCAGCCAGGCCTGGCGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTc                     >  1:204819/1‑149 (MQ=255)
                                                                                                                                gCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCg                    <  4:142769/149‑1 (MQ=255)
                                                                                                                                 cccGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGc                   >  1:184471/1‑149 (MQ=255)
                                                                                                                                  ccGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCt                  <  4:188811/149‑1 (MQ=255)
                                                                                                                                    ggCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGa                <  3:95212/149‑1 (MQ=255)
                                                                                                                                     gCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGaa               <  8:240188/149‑1 (MQ=255)
                                                                                                                                     gCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGaa               <  6:307685/149‑1 (MQ=255)
                                                                                                                                         cgcgCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGc           >  5:138460/1‑149 (MQ=255)
                                                                                                                                          ccaccgCCAAGCCTGACCCAGGCCTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCt          >  1:168779/6‑149 (MQ=255)
                                                                                                                                          ccaccgCCAAGCCTGACCCAGGCCTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCt          >  7:22836/6‑149 (MQ=255)
                                                                                                                                            gcAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCg        >  6:201861/1‑149 (MQ=255)
                                                                                                                                             cAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCgg       <  5:215919/149‑1 (MQ=255)
                                                                                                                                               gCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCGggc     >  4:12463/1‑149 (MQ=255)
                                                                                                                                                  aGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCGggctgg  <  1:367106/149‑1 (MQ=255)
                                                                                                                                                  aGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCGggctgg  <  5:229420/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTTTGATTGCTGCGGGGTTGCAGCAGCAGAATACGTTCACGGTTGAGCCCGGCATCGCGCAGCCAGGCCTGGGTCAGGCTTGGCGGTGGGGCGATCAGGGTCAGCCAGCGGGCCTCGTCCTCTTCGCTCAGTTCGCGCAGCACGGGGGCCAGCAGGCTTTGGCAGTGGCCAGGGGCGCCGCGCAGCGACAGTTCGCTGAACAGCTCGGGCTGG  >  NC_002947/2444483‑2444776

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