Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,803,319 G→A 5.3% D658N (GAT→AAT)  nrdE → ribonucleoside‑diphosphate reductase 2, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,803,3190GA 5.3% 252.9 / 3.4 76D658N (GAT→AAT) nrdEribonucleoside‑diphosphate reductase 2, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (28/44);  new base A (2/2);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

GTAGAGATACGCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATT  >  NC_000913/2803172‑2803465
                                                                                                                                                   |                                                                                                                                                  
gTAGAGATACGCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGa                                                                                                                                                   >  3:502121/1‑149 (MQ=255)
gTAGAGATACGCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGa                                                                                                                                                   <  3:771299/149‑1 (MQ=255)
 tAGAGATACGCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGAt                                                                                                                                                  <  1:270827/149‑1 (MQ=255)
     gaTACGCAAAGAGGGGAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAgg                                                                                                                                              <  1:723184/149‑1 (MQ=255)
     gaTACGCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAgg                                                                                                                                              <  2:421534/149‑1 (MQ=255)
        aCGCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGc                                                                                                                                           >  1:41506/1‑149 (MQ=255)
          gCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGa                                                                                                                                                   >  1:360283/1‑139 (MQ=255)
          gCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGa                                                                                                                                                   <  2:360283/139‑1 (MQ=255)
          gCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTg                                                                                                                                         >  2:239742/1‑149 (MQ=255)
            aaaGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTc                                                                                                                                       <  1:580273/149‑1 (MQ=255)
              agagGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATGTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTATGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGc                                                                                                                                     >  3:759695/1‑149 (MQ=255)
               gagGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCt                                                                                                                                    <  1:631565/149‑1 (MQ=255)
                 gggCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGa                                                                                                                                  <  4:759695/149‑1 (MQ=255)
                  ggCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGAc                                                                                                                                 >  4:65980/1‑149 (MQ=255)
                     aaaaCAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCt                                                                                                                              >  1:54658/1‑149 (MQ=255)
                     aaaaCAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCt                                                                                                                              >  3:406648/1‑149 (MQ=255)
                         cAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttttt                                                                                                                          <  3:191180/149‑1 (MQ=255)
                          aGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTTGCCATGTCGATCAGGGGCTGTCGCTGACGCtttttt                                                                                                                         >  2:173153/1‑149 (MQ=255)
                             aCGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTcc                                                                                                                      <  3:605586/149‑1 (MQ=255)
                             aCGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTcc                                                                                                                      <  1:98202/149‑1 (MQ=255)
                               gcgTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTcccc                                                                                                                    <  4:777640/149‑1 (MQ=255)
                                     actacCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATAcc                                                                                                              >  1:397433/1‑149 (MQ=255)
                                       tacCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGc                                                                                                            >  2:280551/1‑149 (MQ=255)
                                         cccTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGcca                                                                                                          <  1:239742/149‑1 (MQ=255)
                                         cccTGCCCCGTTTATGACTAACGAAAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGcca                                                                                                          <  1:337508/149‑1 (MQ=255)
                                          ccTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGAGGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTccc                                                                                                                     <  4:304325/137‑1 (MQ=255)
                                          ccTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGAGGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTccc                                                                                                                     >  3:304325/1‑137 (MQ=255)
                                          ccTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACTCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTccc                                                                                                                     <  4:679623/137‑1 (MQ=255)
                                          ccTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTccc                                                                                                                     >  3:679623/1‑137 (MQ=255)
                                          ccTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGccac                                                                                                         <  4:406648/149‑1 (MQ=255)
                                          ccTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGccac                                                                                                         <  2:54658/149‑1 (MQ=255)
                                                 cGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTcgc                                                                                                  <  1:190558/149‑1 (MQ=255)
                                                 cGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTc                                                                                                    <  3:660872/147‑1 (MQ=255)
                                                 cGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTc                                                                                                    >  4:660872/1‑147 (MQ=255)
                                                  gTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTcgcg                                                                                                 <  4:570625/149‑1 (MQ=255)
                                                  gTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTcgcg                                                                                                 <  4:7489/149‑1 (MQ=255)
                                                     tATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGata                                                                                              <  3:232682/149‑1 (MQ=255)
                                                      aTGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGatat                                                                                             <  3:65980/149‑1 (MQ=255)
                                                       tGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATc                                                                                            <  3:312499/149‑1 (MQ=255)
                                                       tGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATc                                                                                            <  1:173153/149‑1 (MQ=255)
                                                         aCTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATcaa                                                                                          <  2:41506/149‑1 (MQ=255)
                                                          cTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTc                                                                                                                       <  4:568534/119‑1 (MQ=255)
                                                          cTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTc                                                                                                                       >  3:568534/1‑119 (MQ=255)
                                                           tAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATcaaca                                                                                        >  4:597585/1‑149 (MQ=255)
                                                               gagaATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgc                                                                                    <  4:484972/149‑1 (MQ=255)
                                                                          gcTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGata                                                                                              >  4:532943/1‑128 (MQ=255)
                                                                          gcTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGata                                                                                              <  3:532943/128‑1 (MQ=255)
                                                                               aTCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTgg                                                                    >  4:244048/1‑149 (MQ=255)
                                                                                  aGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgcgc                                                                 >  1:174815/1‑149 (MQ=255)
                                                                                     aCGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAg                                                              <  4:365547/149‑1 (MQ=255)
                                                                                        cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGt                                                           <  2:174815/149‑1 (MQ=255)
                                                                                        cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGt                                                           >  3:234336/1‑149 (MQ=255)
                                                                                         ttACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAg                                                              <  2:43711/145‑1 (MQ=255)
                                                                                         ttACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAg                                                              >  1:43711/1‑145 (MQ=255)
                                                                                         ttACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTa                                                          >  2:800935/1‑149 (MQ=255)
                                                                                           aCGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATc                                                        <  3:104873/149‑1 (MQ=255)
                                                                                                aTTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCaaacc                                                   <  1:382971/149‑3 (MQ=255)
                                                                                                aTTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAAc                                                   <  1:744476/149‑1 (MQ=255)
                                                                                                      gcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTa                                             >  3:616381/1‑149 (MQ=255)
                                                                                                        aGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                           <  2:642904/149‑1 (MQ=255)
                                                                                                        aGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                           <  3:244048/149‑1 (MQ=255)
                                                                                                         gAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTa                                          <  4:502121/149‑1 (MQ=255)
                                                                                                            aaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACAt                                       <  2:397433/149‑1 (MQ=255)
                                                                                                                   tcGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTg                                >  1:675338/1‑149 (MQ=255)
                                                                                                                    cGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGc                               <  4:234336/149‑1 (MQ=255)
                                                                                                                           tACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGt                                                           >  3:796925/1‑114 (MQ=255)
                                                                                                                           tACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGt                                                           <  4:796925/114‑1 (MQ=255)
                                                                                                                            aCGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATg                       >  3:449849/1‑149 (MQ=255)
                                                                                                                             cgcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAATACATCCGCCTGCGTCAGATgg                      <  1:491841/149‑1 (MQ=255)
                                                                                                                              gcgGAAGCGACTCGCCATGTCAATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                           >  4:424576/1‑127 (MQ=255)
                                                                                                                              gcgGAAGCGACTCGCCATGTCAATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                           >  4:551153/1‑127 (MQ=255)
                                                                                                                              gcgGAAGCGACTCGCCATGTCAATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                           <  3:424576/127‑1 (MQ=255)
                                                                                                                              gcgGAAGCGACTCGCCATGTCAATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                           <  3:551153/127‑1 (MQ=255)
                                                                                                                                                gTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAAt   <  1:280551/149‑1 (MQ=255)
                                                                                                                                                gTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAAt   >  3:62607/1‑149 (MQ=255)
                                                                                                                                                 tCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAAtt  >  3:351084/1‑149 (MQ=255)
                                                                                                                                                   gATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCt                                 <  4:302766/116‑1 (MQ=255)
                                                                                                                                                   gATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCt                                 >  3:302766/1‑116 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTAGAGATACGCAAAGAGGGCAAAACAGGACGCGTTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATT  >  NC_000913/2803172‑2803465

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: