Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,814,962 C→T 7.3% P366S (CCC→TCC)  katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,814,9620CT 7.3% 168.5 / 2.8 55P366S (CCC→TCC) katEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base C (33/18);  new base T (2/2);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTA  >  NC_000913/1814814‑1815110
                                                                                                                                                    |                                                                                                                                                    
tttgcACCGCCGCGAGGTGTTTGAAGCCATTGAAGCAGGCGATTTTCCTGAATACGAACTGGGCTTCCAGTTGATTCCTGATGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGTAATAACTGGTGc                                                                                                                                                      <  2:911441/145‑1 (MQ=255)
cTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCTGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGc                                                                                                                                                      <  2:335006/149‑1 (MQ=255)
 ttCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGcc                                                                                                                                                     <  1:767619/149‑1 (MQ=255)
                ttGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCgg                                                                                                                                      >  1:1524179/1‑149 (MQ=255)
                     ggAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCaaaa                                                                                                                                 >  1:1152252/1‑149 (MQ=255)
                     ggAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCaaaa                                                                                                                                 >  2:1087377/1‑149 (MQ=255)
                          ccATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTg                                                                                                                            <  2:1278810/149‑1 (MQ=255)
                           cATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGc                                                                                                                           >  2:30586/1‑149 (MQ=255)
                            aTTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCt                                                                                                                          >  2:764416/1‑149 (MQ=255)
                            aTTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGAGGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCt                                                                                                                          <  2:1178532/149‑1 (MQ=255)
                             ttGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTc                                                                                                                         >  2:962286/1‑149 (MQ=255)
                             ttGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTc                                                                                                                         >  2:1103246/1‑149 (MQ=255)
                                  gCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGTTGTCCGTTCAGAGTGTCGGCAAAATGGTGCTCAATCg                                                                                                                    >  2:580873/1‑149 (MQ=255)
                                  gCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGTCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCg                                                                                                                    >  2:140749/1‑149 (MQ=255)
                                          ttttCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAATCTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCgg                                                                                                            >  1:445592/1‑149 (MQ=255)
                                                  aaTACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAActtc                                                                                                    <  2:896457/149‑1 (MQ=255)
                                                     aCGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAATTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCttt                                                                                                 <  2:1121060/149‑1 (MQ=255)
                                                      cGATCTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTg                                                                                                <  2:1238217/149‑1 (MQ=255)
                                                      cGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTg                                                                                                <  1:746015/149‑1 (MQ=255)
                                                       gAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACACGGATAACTTCTTTGc                                                                                               >  2:1214253/1‑149 (MQ=255)
                                                           tGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTg                                                                                                <  1:1357467/144‑1 (MQ=255)
                                                           tGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTg                                                                                                >  2:1357467/1‑144 (MQ=255)
                                                             ggCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCaaaa                                                                                                                                 >  2:268596/1‑109 (MQ=255)
                                                             ggCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCaaaa                                                                                                                                 <  1:268596/109‑1 (MQ=255)
                                                               cTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACg                                                                                       >  1:833148/1‑149 (MQ=255)
                                                                ttCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGa                                                                                      >  1:577927/1‑149 (MQ=255)
                                                                   cAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACa                                                                                   >  1:1584127/1‑149 (MQ=255)
                                                                   cAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGATCAGCGTGTCGGCAAAATGGTGCTCAATCGCTACCCGGATAACTTCTTTTGTGAAAACGaaaa                                                                                   >  2:110768/1‑147 (MQ=255)
                                                                     gTTGATTCCTGAAGTAGATGAATTCAAGTTCGACATCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAgg                                                                                 <  2:815624/149‑1 (MQ=255)
                                                                     gTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAgg                                                                                 >  2:906766/1‑149 (MQ=255)
                                                                             cTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTc                                                                         >  2:337071/1‑149 (MQ=255)
                                                                                aagaagATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCAt                                                                      >  1:192963/1‑149 (MQ=255)
                                                                                aagaagATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCAt                                                                      >  1:1182676/1‑149 (MQ=255)
                                                                                  caagaTGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATcc                                                                    <  1:1321211/148‑1 (MQ=255)
                                                                                   aagaTGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCt                                                                   <  1:1392228/149‑1 (MQ=255)
                                                                                     gaTGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTgg                                                                 >  2:197280/1‑149 (MQ=255)
                                                                                      aTGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGACCAGGCGGCTTTCCATCCTggg                                                                >  2:1144151/1‑149 (MQ=255)
                                                                                           ttCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCatat                                                           >  2:731186/1‑149 (MQ=255)
                                                                                           ttCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCatat                                                           >  1:748107/1‑149 (MQ=255)
                                                                                                          gATCATCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGCAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGggactcga                                            >  1:717558/1‑149 (MQ=255)
                                                                                                           aTCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGTCCGTTCAGCGTGTCGTCAAAATGGTGCTAAATCGCAACCCGGATAACTTCTTTGCTGAAAACTAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGggactggac                                           <  1:580873/149‑1 (MQ=255)
                                                                                                           aTCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGTCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGggactggac                                           <  1:140749/149‑1 (MQ=255)
                                                                                                           aTCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGggactggac                                           >  2:1287134/1‑149 (MQ=255)
                                                                                                           aTCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGggactggac                                           <  1:1214253/149‑1 (MQ=255)
                                                                                                            tcttctCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGggactggact                                          >  2:1506425/1‑149 (MQ=255)
                                                                                                              ttctCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTc                                        <  2:642925/149‑1 (MQ=255)
                                                                                                                ctcGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCAc                                      >  1:1437018/1‑149 (MQ=255)
                                                                                                                ctcGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCAc                                      >  1:953790/1‑149 (MQ=255)
                                                                                                                ctcGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCGGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCAc                                      <  1:1103246/149‑1 (MQ=255)
                                                                                                                    aTCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAAc                                  <  2:598966/149‑1 (MQ=255)
                                                                                                                       caaccaaACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCTTATCGTGCCGGGACTGGACTTCACCAACGAt                               <  2:577927/149‑1 (MQ=255)
                                                                                                                         accaaACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATcc                             >  2:1237892/1‑149 (MQ=255)
                                                                                                                          ccaaACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCg                            >  2:596686/1‑149 (MQ=255)
                                                                                                                            aaaCTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGCTCCGCt                          >  2:163254/1‑149 (MQ=255)
                                                                                                                                ttATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTg                      >  1:1376634/1‑149 (MQ=255)
                                                                                                                                  aTCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTGCa                    <  1:764416/149‑1 (MQ=255)
                                                                                                                                       ggAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGAc               >  1:53798/1‑149 (MQ=255)
                                                                                                                                           caaCTGGTGCCCGTTCAGCGTGGCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTATGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGttt           <  1:1018660/148‑1 (MQ=255)
                                                                                                                                                   ccccGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCt   <  2:1152252/148‑1 (MQ=255)
                                                                                                                                                    cccGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCta  >  1:1085328/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTCAAGTTCGACTTCGATCTTCTCGATCCAACCAAACTTATCCCGGAAGAACTGGTGCCCGTTCAGCGTGTCGGCAAAATGGTGCTCAATCGCAACCCGGATAACTTCTTTGCTGAAAACGAACAGGCGGCTTTCCATCCTGGGCATATCGTGCCGGGACTGGACTTCACCAACGATCCGCTGTTGCAGGGACGTTTGTTCTCCTA  >  NC_000913/1814814‑1815110

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: