Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,272 T→G 6.4% S16A (TCC→GCC)  iscA → copper(I) binding iron‑sulfur cluster assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,2720TG 6.4% 240.4 / 6.7 78S16A (TCC→GCC) iscAcopper(I) binding iron‑sulfur cluster assembly protein
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base G (3/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

ACTGCTCGGTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCG  >  NC_002947/980125‑980419
                                                                                                                                                    |                                                                                                                                                    
aCTGCTCGGTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTc                                                                                                                                                     >  8:222943/1‑149 (MQ=255)
 cTGCTCGGTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTcc                                                                                                                                                    >  1:353731/1‑149 (MQ=255)
    cTCGGTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTccc                                                                                                                                                   <  3:184656/147‑1 (MQ=255)
       ggTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGa                                                                                                                                              <  5:212555/148‑1 (MQ=255)
       ggTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGa                                                                                                                                              <  7:163925/148‑1 (MQ=255)
       ggTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGAc                                                                                                                                             <  3:57083/149‑1 (MQ=255)
               ccGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCGAACCACATTCGCCGCTCCCTCGACGGGCGCgg                                                                                                                                     <  1:235769/149‑1 (MQ=255)
               ccGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCgg                                                                                                                                     <  7:222943/149‑1 (MQ=255)
                    gATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCGCCCGTCGAGGGAGCGGCgaat                                                                                                                                <  2:302017/149‑2 (MQ=255)
                      tGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGt                                                                                                                              >  8:196673/1‑149 (MQ=255)
                         cATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAg                                                                                                                           <  3:39343/149‑1 (MQ=255)
                          aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGcgcg                                                                                                                                      <  1:359255/137‑1 (MQ=255)
                          aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGcgcg                                                                                                                                      >  2:359255/1‑137 (MQ=255)
                          aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAgg                                                                                                                          >  3:37407/1‑149 (MQ=255)
                          aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAgg                                                                                                                          <  2:88738/149‑1 (MQ=255)
                          aTCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAgg                                                                                                                          <  7:123104/149‑1 (MQ=255)
                           tCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAggg                                                                                                                         >  6:73464/1‑149 (MQ=255)
                            cAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGcgcg                                                                                                                                      <  3:48727/135‑1 (MQ=255)
                            cAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGcgcg                                                                                                                                      >  4:48727/1‑135 (MQ=255)
                                gCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTc                                                                                                                    >  5:76065/1‑149 (MQ=255)
                                     cttaggcgAGTACTAGCAGAAGATAGGCTTGATCTAAGTCGCCTGTTGCAAGGTATGGAGTTCTGTTGGCTCTCAGGATGACAGAAGCCGACGCCAACCACATTCGCCTCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCt                                                                                                                <  3:280101/143‑1 (MQ=255)
                                      gTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTgg                                                                                                              <  7:290984/149‑1 (MQ=255)
                                      gTACGCGACTACAAGCAGAAGAAAGCCTTGATCTACGTCGCCTGTTGCAAGGTAAGGAGGTCTGATGGCTACCAGCATGACAGAAGCCGCCCCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTgg                                                                                                              <  4:107131/149‑1 (MQ=255)
                                          gcgACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGt                                                                                                          <  1:56767/149‑1 (MQ=255)
                                            gACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTgc                                                                                                        <  7:202360/149‑1 (MQ=255)
                                                  aaGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGcacca                                                                                                  >  5:191883/1‑149 (MQ=255)
                                                     cAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCg                                                                                               >  3:288505/1‑149 (MQ=255)
                                                       gaagaaAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCGCCCGTCGAGGGAGCGGCgaat                                                                                                                                <  1:239378/114‑2 (MQ=38)
                                                       gaagaaAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCGCCCGTCGAGGGAGCGGCgaat                                                                                                                                >  2:239378/1‑113 (MQ=38)
                                                          gaaAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGc                                                                                          >  6:165552/1‑149 (MQ=255)
                                                                cTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGt                                                                                    >  6:327636/1‑149 (MQ=255)
                                                                   gATGTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGGGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCTCCACCGGTTGCTCgg                                                                                      <  3:174779/144‑1 (MQ=255)
                                                                   caTCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCCCCGCTCACTCGACGGGCGAGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTAGCTCgggagaa                                                                                 >  4:174779/2‑145 (MQ=255)
                                                                    aTCTAAGTCGCCTGTTTCAAGGTAAGGCGTTCTGTTGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTg                                                                                 <  4:372943/148‑1 (MQ=255)
                                                                    aTCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGAAAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTgg                                                                                >  7:125724/1‑149 (MQ=255)
                                                                     tCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGc                                                                               >  5:113219/1‑149 (MQ=255)
                                                                         aGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGc                                                                               <  4:37407/145‑1 (MQ=255)
                                                                         aGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTa                                                                            <  6:76065/148‑1 (MQ=255)
                                                                         aGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTAc                                                                           <  2:353731/149‑1 (MQ=255)
                                                                           tCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGt                                                                         >  1:46232/1‑149 (MQ=255)
                                                                              ccTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCt                                                                      <  5:288273/149‑1 (MQ=255)
                                                                                tGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTg                                                                     >  2:327412/1‑148 (MQ=255)
                                                                                tGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTg                                                                     >  8:341140/1‑148 (MQ=255)
                                                                                  ttGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGaa                                                                  >  4:200123/1‑149 (MQ=255)
                                                                                    gCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAAtt                                                                >  7:18821/1‑149 (MQ=255)
                                                                                     cAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATtc                                                               <  7:196673/149‑1 (MQ=255)
                                                                                     cAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATtc                                                               >  3:112204/1‑149 (MQ=255)
                                                                                        ggTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATtcgtc                                                            <  7:341140/149‑1 (MQ=255)
                                                                                         gTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATtcgtc                                                            >  6:90676/1‑148 (MQ=255)
                                                                                          tAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTcga                                                          >  5:33100/1‑149 (MQ=255)
                                                                                                gTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACt                                                    >  7:168018/1‑149 (MQ=255)
                                                                                                 ttCTGATGGCTATCAGCATGACAGAAGCCGCCGCGCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACt                                                    >  1:247201/1‑149 (MQ=255)
                                                                                                    tGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCg                                                >  4:181407/1‑149 (MQ=255)
                                                                                                    tGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCg                                                >  4:175468/1‑149 (MQ=255)
                                                                                                      aTGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGc                                                 <  8:168018/146‑1 (MQ=255)
                                                                                                      aTGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGa                                              >  8:244869/1‑149 (MQ=255)
                                                                                                            aTCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCgg                                               <  8:140114/142‑1 (MQ=255)
                                                                                                            aTCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCgg                                               >  7:140114/1‑142 (MQ=255)
                                                                                                             tCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGAc                                       <  3:200123/149‑1 (MQ=255)
                                                                                                                gCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAg                                    <  2:46232/149‑1 (MQ=255)
                                                                                                                 cATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAgg                                   >  6:32334/1‑149 (MQ=255)
                                                                                                                  aTGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGgt                                  >  5:342015/1‑149 (MQ=255)
                                                                                                                   tGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGgtg                                 <  1:327412/149‑1 (MQ=255)
                                                                                                                      cAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTc                              >  4:151504/1‑149 (MQ=255)
                                                                                                                                  ccAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGa                                        >  3:230421/1‑127 (MQ=255)
                                                                                                                                  ccAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGa                                        <  4:230421/127‑1 (MQ=255)
                                                                                                                                     aaCCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGc                                                                               <  7:49287/86‑1 (MQ=255)
                                                                                                                                     aaCCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGc                                                                               >  8:49287/1‑86 (MQ=255)
                                                                                                                                      ctcaccTTTGCCGCGCCCGTCGAGGGAGCGGCGAATGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGt                >  3:274118/3‑148 (MQ=37)
                                                                                                                                       tcaccTTTGCCGCGCCCGTCGAGGGAGCGGCGAATGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTc               >  7:111572/2‑148 (MQ=38)
                                                                                                                                        cacaTTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTc               >  1:295022/1‑147 (MQ=255)
                                                                                                                                          caTTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAgg           <  6:191883/149‑1 (MQ=255)
                                                                                                                                          caTTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAgg           <  5:90676/149‑1 (MQ=255)
                                                                                                                                              cgccgcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGa        <  6:113219/148‑1 (MQ=255)
                                                                                                                                               gccgcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGat       >  6:314165/1‑148 (MQ=255)
                                                                                                                                               gccgcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGat       <  8:125724/148‑1 (MQ=255)
                                                                                                                                                  gcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGatcatc   <  5:314165/149‑1 (MQ=255)
                                                                                                                                                   cTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCCGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCg  >  8:283206/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTGCTCGGTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCG  >  NC_002947/980125‑980419

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: