Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,280 Δ1 bp 6.5% coding (54/324 nt) 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,2800C. 6.5% 248.1 / 9.9 77coding (54/324 nt)iscAcopper(I) binding iron‑sulfur cluster assembly protein
Reads supporting (aligned to +/- strand):  ref base C (38/34);  new base . (3/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

GGTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATCCC  >  NC_002947/980132‑980424
                                                                                                                                                      |                                                                                                                                                
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                                                                                                                                     >  2:239378/1‑113 (MQ=38)
                                                gaagaaAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCGCCCGTCGAGGGAGCGGCgaat                                                                                                                                     <  1:239378/114‑2 (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                                                                          >  8:341140/1‑148 (MQ=255)
                                                                         tGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTg                                                                          >  2:327412/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            <  8:125724/148‑1 (MQ=255)
                                                                                                                                        gccgcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGat            >  6:314165/1‑148 (MQ=255)
                                                                                                                                           gcTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGatcatc        <  5:314165/149‑1 (MQ=255)
                                                                                                                                            cTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCCGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCg       >  8:283206/1‑149 (MQ=255)
                                                                                                                                               ccTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATc    >  4:15932/1‑149 (MQ=255)
                                                                                                                                                cTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATcc   >  4:21867/1‑149 (MQ=255)
                                                                                                                                                  tCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATcc   >  5:380766/1‑148 (MQ=255)
                                                                                                                                                    gACGGGCGCGTCAAAGGTGAGGGCATTCTCCTGTGCGTGCTCACCACCGGTTGCTCTGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATccc  <  3:281800/147‑1 (MQ=255)
                                                                                                                                                    gACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATccc  <  5:165552/147‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GGTACTCGCCGAAGATGCCATCAAGGCAGCCGTACGCGACTACAAGCAGAAGAAAGGCTTGATCTAAGTCGCCTGTTGCAAGGTAAGGAGTTCTGATGGCTATCAGCATGACAGAAGCCGCCGCCAACCACATTCGCCGCTCCCTCGACGGGCGCGGCAAAGGTGAGGGCATTCGCCTGGGCGTGCGCACCACCGGTTGCTCGGGTCTGGCCTACGTGCTGGAATTCGTCGACGAACTGGCGGACGAAGACCAGGTGTTCGAAAACCATGGCGTCAAGGTGATCATCGATCCC  >  NC_002947/980132‑980424

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: