Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,409,994 G→C 6.4% F499L (TTC→TTG acnA‑I ← aconitate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,409,9940GC 6.4% 193.1 / 4.7 62F499L (TTC→TTGacnA‑Iaconitate hydratase
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base C (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

AGCACCGAGGCCACGGTGAGGTCGGCGCTGCCAATGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACGTGTTGGTGCATGAGGTGATC  >  NC_002947/2409845‑2410141
                                                                                                                                                     |                                                                                                                                                   
ctcACCGTGGCCTCGGTGAGGTCGGCGCTGCCAATGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTg                                                                                                                                                     >  2:1878627/3‑150 (MQ=255)
  cACCGAGGCCACGGTGAGGTCGGCGCTGCCAATGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGaa                                                                                                                                                   <  1:2216053/150‑1 (MQ=255)
       aGGCCACGGTGAGGTCGGCGCTGCCAATGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAagtagt                                                                                                                                              <  2:725812/150‑1 (MQ=255)
       aGGCCACGGTGAGGTCGGCGCTGCCAATGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAagtagt                                                                                                                                              <  1:1824732/150‑1 (MQ=255)
          ccACGGCGAGGTCGGCGCTGCCAATGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCAGCTGCCTTGAAGTAGTCAg                                                                                                                                           <  1:852826/150‑1 (MQ=255)
                          gcTGCCAATGACTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGGAGTAGTCAGTGACGACTTTGGAGcc                                                                                                                           >  2:642014/1‑150 (MQ=255)
                            tGCCAATGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCg                                                                                                                         <  1:1395324/150‑1 (MQ=255)
                              ccAATGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGACGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGCGCCCTgc                                                                                                                       >  2:168613/1‑150 (MQ=255)
                                aaTGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGgcgc                                                                                                                     <  2:1144220/150‑1 (MQ=255)
                                  tGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCa                                                                                                                   <  1:2356513/150‑1 (MQ=255)
                                            gATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCtt                                                                                                         >  1:363058/1‑150 (MQ=255)
                                              tCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGa                                                                                                       <  2:1440873/150‑1 (MQ=255)
                                                gCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCGGCCATAGCCGACCCGGTCACAACCCAGTTGGACGAGGTAAGGCGTGATCCCGGCTGCCTTGAAGTAGTCAGTGGCGAGTTTGAAGGCCGGCGCCAGCGAAGATTTGAcc                                                                                                     >  2:139875/1‑150 (MQ=255)
                                                 ccTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGTCCGGCGCCAGCGAACTCTTGAccc                                                                                                    >  2:1548036/1‑150 (MQ=255)
                                                        cAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGTCCGGCGCCAGCGAACTCTTGACCCAAGGttt                                                                                             <  1:1548036/150‑1 (MQ=255)
                                                        cAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGttt                                                                                             >  2:1204811/1‑150 (MQ=255)
                                                        cAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGttt                                                                                             >  1:1004220/1‑150 (MQ=255)
                                                          gCGGGCCGGAGGTGCGGATGCAGGTGGTGAAGCCCACATCGTTCAGACGACAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGCCTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGc                                                                                           <  2:193587/150‑1 (MQ=255)
                                                            gggCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGt                                                                                         <  2:1004220/150‑1 (MQ=255)
                                                            gggCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGt                                                                                         <  1:131971/150‑1 (MQ=255)
                                                                     ttGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGccc                                                                                <  2:363058/150‑1 (MQ=255)
                                                                          gATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTc                                                                           >  1:939254/1‑150 (MQ=255)
                                                                               aGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGCAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTg                                                                      >  1:1036827/1‑150 (MQ=255)
                                                                                ggtggtGCAGCCATAGCGGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGc                                                                     >  2:1237368/1‑150 (MQ=255)
                                                                                    gtgCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGGAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCtttt                                                                 <  1:642014/150‑1 (MQ=255)
                                                                                    gtgCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCtttt                                                                 <  1:1878627/150‑1 (MQ=255)
                                                                                            aTAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGccac                                                         <  1:1659647/150‑1 (MQ=255)
                                                                                                 cGACCAGATCAAAACCCAGTTGGTCGAGGGAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGcc                                                    <  1:2344965/150‑1 (MQ=255)
                                                                                                     cAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGc                                                                                           >  2:2219357/1‑107 (MQ=255)
                                                                                                     cAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGc                                                                                           <  1:2219357/107‑1 (MQ=255)
                                                                                                      aGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCa                                               >  1:1122559/1‑150 (MQ=255)
                                                                                                        aTCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTg                                                                                                                                <  1:809828/67‑1 (MQ=255)
                                                                                                        aTCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTg                                                                                                                                >  2:809828/1‑67 (MQ=255)
                                                                                                        aTCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGc                                             >  1:506841/1‑150 (MQ=255)
                                                                                                          cAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCca                                           >  2:2358392/1‑150 (MQ=255)
                                                                                                           aaaaCCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCcat                                          <  1:1204811/150‑1 (MQ=255)
                                                                                                               cccAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGccacca                                                       <  1:459610/133‑1 (MQ=255)
                                                                                                               cccAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGccacca                                                       >  2:459610/1‑133 (MQ=255)
                                                                                                                     tGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACCTTGGAGCCCCGCGCCAGCGAACCCTTGACACCAGGCTTGAGTTGCAGGTCCTTTTCCAGTGCGCTTTTCTCCTCCAGCCCGGCAGCCATCATCACACTg                                <  2:144426/150‑1 (MQ=255)
                                                                                                                        tCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAgg                                                                                   <  2:2112490/96‑1 (MQ=255)
                                                                                                                        tCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAgg                                                                                   >  1:2112490/1‑96 (MQ=255)
                                                                                                                          acgctAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGtt                           >  2:2256856/5‑150 (MQ=255)
                                                                                                                             gTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACtt                                                                                                                                  <  2:952917/44‑1 (MQ=255)
                                                                                                                             gTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACtt                                                                                                                                  >  1:952917/1‑44 (MQ=255)
                                                                                                                               aaGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACg                      >  1:1337807/1‑150 (MQ=255)
                                                                                                                               aaGGCGTGAGCCCGGCTGCCTTCAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGAc                       >  2:1003111/1‑149 (MQ=255)
                                                                                                                               aaGGCGTGAGCCCGGCTGCCTTCAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGAc                       <  1:1002937/149‑1 (MQ=255)
                                                                                                                               aaGGCGTGAGCCCGGCTGCCTTCAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGAc                       <  1:1003111/149‑1 (MQ=255)
                                                                                                                               aaGGCGTGAGCCCGGCTGCCTTCAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGAc                       >  2:1002937/1‑149 (MQ=255)
                                                                                                                                aGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCAAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACgt                     >  1:444212/1‑150 (MQ=255)
                                                                                                                                  gCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCtt                                                                                                         <  1:2375815/64‑1 (MQ=255)
                                                                                                                                  gCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCtt                                                                                                         >  2:2375815/1‑64 (MQ=255)
                                                                                                                                  gCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACgtgt                   <  1:2358392/150‑1 (MQ=255)
                                                                                                                                  gCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACgtgt                   >  1:2288791/1‑150 (MQ=255)
                                                                                                                                   cGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCa                                                                                     <  1:2340113/83‑1 (MQ=255)
                                                                                                                                   cGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCa                                                                                     >  2:2340113/1‑83 (MQ=255)
                                                                                                                                    gTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACGTGTTg                 >  1:1450077/1‑150 (MQ=255)
                                                                                                                                         cccGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACGTGTTGGTGCa            >  1:1008946/1‑150 (MQ=255)
                                                                                                                                          ccGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACGTGTTGGTGCAt           <  2:444212/150‑1 (MQ=255)
                                                                                                                                          ccGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGACCTCGGCTCTATGGAAGTCTTGACCCAAGGTTTGCGTTTCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACGTGTTGGTGCAt           <  1:540859/150‑1 (MQ=255)
                                                                                                                                           cGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACGTGTTGGTGCATg          >  1:1876089/1‑150 (MQ=255)
                                                                                                                                              cTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACGTGTTGGTGCATGAgg       <  1:1237368/150‑1 (MQ=255)
                                                                                                                                                   ttGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACGTGTTGGTGCATGAGGTGATc  <  2:1008946/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGCACCGAGGCCACGGTGAGGTCGGCGCTGCCAATGGCTTTTTCGATCGCCTCGTCCAGCGGGCCGGAGTTGCCGATGCAGGTGGTGCAGCCATAGCCGACCAGATCAAAACCCAGTTGGTCGAGGTAAGGCGTGAGCCCGGCTGCCTTGAAGTAGTCAGTGACGACTTTGGAGCCCGGCGCCAGCGAACTCTTGACCCAAGGTTTGCGTTGCAGGCCCTTTTCCAGTGCTTTTTTCGCCACCAGCCCGGCAGCCATCATCACACTGGGGTTGGACGTGTTGGTGCATGAGGTGATC  >  NC_002947/2409845‑2410141

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