Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,793,103 C→G 7.0% G194A (GGC→GCC)  yibB ← YibB family protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,793,1030CG 7.0% 454.6 / 14.7 135G194A (GGC→GCC) yibBYibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base C (64/61);  new base G (6/4);  total (70/65)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 141 total aligned reads displayed.

GTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATAGCCAAAATCTAT  >  NC_000913/3792955‑3793251
                                                                                                                                                    |                                                                                                                                                    
  aaTAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCa                                                                                                                                                    <  2:369670/149‑1 (MQ=255)
  aaTAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCa                                                                                                                                                    >  2:709021/1‑149 (MQ=255)
    tAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAAt                                                                                                                                                  >  2:1344213/1‑149 (MQ=255)
        acaTGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATg                                                                                                                                              >  1:1093729/1‑149 (MQ=255)
        acaTGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATg                                                                                                                                              <  1:482828/149‑1 (MQ=255)
        acaTGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATg                                                                                                                                              <  1:1463691/149‑1 (MQ=255)
        acaTGACAAAAATGCCCTGGTCATCGTCAACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATg                                                                                                                                              <  1:58394/149‑1 (MQ=255)
           tGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGaaa                                                                                                                                           >  1:1356223/1‑149 (MQ=255)
             aCAATAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATc                                                                                                                                         <  1:255658/149‑1 (MQ=255)
             aCAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATc                                                                                                                                         <  2:1093729/149‑1 (MQ=255)
              cAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGGCAATCATgtatata                                                                                                                                        <  1:1539938/149‑7 (MQ=255)
               aaaaaTGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCaa                                                                                                                                       >  2:643184/1‑149 (MQ=255)
                  aaTGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCaaaaa                                                                                                                                    <  2:914431/149‑1 (MQ=255)
                  aaTGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTACAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCaaaaa                                                                                                                                    <  1:408588/149‑1 (MQ=255)
                   aTGCCCTGGTCATCTTCCACAATCTTATTATTTATTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGATGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAAc                                                                                                                                   <  1:65038/149‑1 (MQ=255)
                    tGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAAcc                                                                                                                                  <  1:1354031/149‑1 (MQ=255)
                      cccTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCtg                                                                                                                                <  2:427019/149‑1 (MQ=255)
                            tCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGttgttg                                                                                                                          <  1:699732/149‑1 (MQ=255)
                            tCATCGGACACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGttgttg                                                                                                                          <  2:1386803/149‑1 (MQ=255)
                             cATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTga                                                                                                                         <  2:713066/149‑1 (MQ=255)
                             cATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTACAGAACTAATTTATAGAATTCTTTACATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTga                                                                                                                         <  1:1160941/149‑1 (MQ=255)
                              atacaacacaATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTgag                                                                                                                        <  1:551298/143‑1 (MQ=255)
                              aTCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTgag                                                                                                                        >  2:1270843/1‑149 (MQ=255)
                              aTCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTgag                                                                                                                        <  1:1490266/149‑1 (MQ=255)
                                 gTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGAtg                                                                                                                     >  1:608585/1‑149 (MQ=255)
                                  tCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCATTAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGAtgt                                                                                                                    <  1:885877/149‑1 (MQ=255)
                                  tCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGAtgt                                                                                                                    >  1:569287/1‑149 (MQ=255)
                                  tCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACAAACAATTGCGCCCCCAATAATATATACATTATTGCCAATCATGAAATCAAAAAACTGTTGTTGAGAtgt                                                                                                                    >  1:825474/1‑149 (MQ=255)
                                   ccACAATGTTATTTATCAGTGTTATTTTTTGACTTTCCAGAACTAACTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAATCCTGTTGTTGAGAtgtg                                                                                                                   <  2:1253138/149‑1 (MQ=255)
                                    cacaATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGa                                                                                                                  <  2:1289473/149‑1 (MQ=255)
                                          gTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTaa                                                                                                            <  1:1114964/149‑1 (MQ=255)
                                           ttattaTTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAg                                                                                                           >  2:703579/1‑149 (MQ=255)
                                           ttattaTTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAg                                                                                                           <  2:1523730/149‑1 (MQ=255)
                                            tattatTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGc                                                                                                          <  2:620354/149‑1 (MQ=255)
                                            tattatTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGc                                                                                                          >  1:716009/1‑149 (MQ=255)
                                             attattTAGTGTTATTTTTTGTCTTTCCAGAAATAATTTATAGAATTCTTTCCATTTATGCTGTTAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                                                                                         <  2:1047105/149‑1 (MQ=255)
                                                 tttAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCttt                                                                                                     >  2:262256/1‑149 (MQ=255)
                                                      tgtTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTta                                                                                                <  1:1523560/149‑1 (MQ=255)
                                                      tgtTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTta                                                                                                <  2:376482/149‑1 (MQ=255)
                                                      tgtTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTAGTTGAGATGAGACAGTAAGAACtttttttt                                                                                                >  2:531884/1‑148 (MQ=255)
                                                        ttATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtata                                                                                              <  1:1530307/149‑1 (MQ=255)
                                                        ttATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtata                                                                                              <  2:525437/149‑1 (MQ=255)
                                                          aTTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGt                                                                                            >  1:737143/1‑149 (MQ=255)
                                                          aTTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGt                                                                                            <  1:709021/149‑1 (MQ=255)
                                                              tttGACTTTCCAGAACTAATTTATAGAAATCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCAAATAATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGACCTATATATATAGTGaaa                                                                                        >  1:1426255/1‑149 (MQ=255)
                                                               ttGACTTTCCAGAACTAATTTATAGAATTATTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGaaaa                                                                                       <  1:984750/149‑1 (MQ=255)
                                                                  aCTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGAt                                                                                    <  2:1357283/149‑1 (MQ=255)
                                                                      tCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGGGACCCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGCTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCAt                                                                                <  2:1115806/149‑1 (MQ=255)
                                                                          gAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTa                                                                            >  1:181180/1‑149 (MQ=255)
                                                                          gAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTa                                                                            >  2:1151340/1‑149 (MQ=255)
                                                                             atAATTTATAGTATTCTTTAGATTTATGCTGTTAACCAAGAATTGCGTCCTCAATAATATATACATGATTGCGAATCATGAAATCAAAAACCTGTTGTTGTGATGTGATAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt                                                                         <  2:799796/148‑1 (MQ=255)
                                                                              tAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACttt                                                                        >  2:1269461/1‑149 (MQ=255)
                                                                               aaTTTATAGAATTCATTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc                                                                       <  1:1167616/149‑1 (MQ=255)
                                                                                aTTTATAGAATTCTTTAAATTTATGCTGTGAACCAACAATTGAGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtca                                                                      <  1:756219/149‑1 (MQ=255)
                                                                                          ttCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGa                                                            >  1:1362421/1‑149 (MQ=255)
                                                                                           tCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGaa                                                           <  2:1356223/149‑1 (MQ=255)
                                                                                             tttCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAAt                                                         >  1:879728/1‑149 (MQ=255)
                                                                                             tttCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAAt                                                         >  2:1005621/1‑149 (MQ=255)
                                                                                               tCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATcc                                                       >  2:291343/1‑149 (MQ=255)
                                                                                                ccATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATccc                                                      >  2:473325/1‑149 (MQ=255)
                                                                                                caaTTTATGCTGAGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGATATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATccc                                                      <  2:576994/147‑1 (MQ=255)
                                                                                                  aTTTATGCTGTGAACCAACAATTGCGCCCACAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGAAAGTAAGACCTTTTTTTATAGTGAAAAGATTCATCTTACTTTCATCAAATGGGAAATCCCaa                                                    >  1:94997/1‑149 (MQ=255)
                                                                                                       tGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAAtttt                                               >  2:1546132/1‑149 (MQ=255)
                                                                                                       tGCTGTGAACCAAAAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAAtttt                                               <  2:1021519/149‑1 (MQ=255)
                                                                                                            tgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATc                                          >  2:496519/1‑149 (MQ=255)
                                                                                                              aaCCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCt                                        >  2:187750/1‑149 (MQ=255)
                                                                                                              aaCCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGCAAATCCCAAATTTTCAATCCt                                        >  1:652855/1‑149 (MQ=255)
                                                                                                                ccAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCg                                      >  2:606682/1‑149 (MQ=255)
                                                                                                                ccAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCg                                      >  1:249005/1‑149 (MQ=255)
                                                                                                                ccAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCg                                      >  1:40655/1‑149 (MQ=255)
                                                                                                                ccAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCAAAATTTTCAATCCTCg                                      >  1:1052642/1‑149 (MQ=255)
                                                                                                                 caacaaTTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGa                                     <  2:1295251/149‑1 (MQ=255)
                                                                                                                  aacaaTTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTATTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAg                                    >  2:853196/1‑149 (MQ=255)
                                                                                                                     aaTTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTa                                 >  2:494291/1‑149 (MQ=255)
                                                                                                                      aTTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTAc                                <  1:606682/149‑1 (MQ=255)
                                                                                                                            cccccAATTATATTTACATGATTGGCAATCATTAAATCAAAAACCTGGTGTTGAGAGGTGACAGTAAGGCCTTTTTTTATAGTGTTAAGATGCATCTTACTTTAATCAAATGGGGAATCCCAAATTTTCAATCCTCGAGTTACATTTgg                          <  2:762807/149‑1 (MQ=255)
                                                                                                                            cccccAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTgg                          <  2:223064/149‑1 (MQ=255)
                                                                                                                            cccccAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTgg                          >  1:677339/1‑149 (MQ=255)
                                                                                                                              cccAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAACATGAATATTACTTTCATCAAATGGAAAATCCCAAATTTTCAATCCAAGATTTACATTTggat                        >  1:175353/1‑147 (MQ=255)
                                                                                                                               ccAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCtt                       >  2:168608/1‑149 (MQ=255)
                                                                                                                                cAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAATTCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTa                      >  1:1457244/1‑149 (MQ=255)
                                                                                                                                 aataatATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTAt                     >  1:1279438/1‑149 (MQ=255)
                                                                                                                                 aataatATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTAt                     >  1:472877/1‑149 (MQ=255)
                                                                                                                                  ataataTATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTAtt                    >  2:1092101/1‑148 (MQ=255)
                                                                                                                                   taataTATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATgg                   <  1:119731/149‑1 (MQ=255)
                                                                                                                                    aataTATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGc                  <  1:863397/149‑1 (MQ=255)
                                                                                                                                       atataGATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAAt               <  2:1426255/149‑1 (MQ=255)
                                                                                                                                       atataCATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAAt               <  1:712277/149‑1 (MQ=255)
                                                                                                                                        tataCATGATTGCCAATCGTGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATa              >  1:476083/1‑149 (MQ=255)
                                                                                                                                        tataCATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATa              >  1:1312753/1‑149 (MQ=255)
                                                                                                                                        atttcATGATTGGCAATCATGCAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTTTAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAATTTTTCAATCCTCGAGTTACATTTGGCTTATGGCAAtt              >  1:1063595/5‑148 (MQ=255)
                                                                                                                                        atttcATGATTGGCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATa              >  1:524186/5‑149 (MQ=255)
                                                                                                                                        atttcATGATTGGCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATa              >  1:672704/5‑149 (MQ=255)
                                                                                                                                         tttcATGATTGGCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATAg             >  2:24139/4‑149 (MQ=255)
                                                                                                                                          ttcATGATTGGCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATAGc            >  2:513030/3‑149 (MQ=255)
                                                                                                                                           aCATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATAGcc           >  1:674547/1‑149 (MQ=255)
                                                                                                                                                    ccAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATAGCCAAAATCTAt  >  1:1062399/1‑149 (MQ=255)
                                                                                                                                                    ccAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATAGCCAAAATCTAt  >  2:575627/1‑149 (MQ=255)
                                                                                                                                                    ccAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATAGCCAAAATCTAt  >  1:52515/1‑149 (MQ=255)
                                                                                                                                                    ccAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACATTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATAGCCAAAATCTAt  <  1:262256/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAATTTTCAATCCTCGAGTTACATTTGGCTTATGGCAATAGCCAAAATCTAT  >  NC_000913/3792955‑3793251

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