Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,712,518 (C)6→8 intergenic (+360/‑251) nth → / → dtpA DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,712,5121.C94.6% 317.3 / 7.7 93intergenic (+354/‑257)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base C (49/39);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
*NC_0009131,712,5122.C93.5% 315.0 / 10.9 93intergenic (+354/‑257)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base C (49/38);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 106 total aligned reads displayed.

AATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTG  >  NC_000913/1712375‑1712640
                                                                                                                                            ||                                                                                                                                   
aaTCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑c                                                                                                                                   <  2:933738/139‑1 (MQ=255)
  tCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTAAAAAACATTACACTGGCTATTTTTCAGAAACTGGATTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TAA‑‑‑cc‑‑c                                                                                                                               <  1:299252/139‑1 (MQ=255)
     ggTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑cc‑‑cccc                                                                                                                            <  1:75108/139‑1 (MQ=255)
     ggTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑cc‑‑cccc                                                                                                                            <  2:79338/139‑1 (MQ=255)
      gTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑cc‑‑ccccc                                                                                                                           >  2:1374222/1‑138 (MQ=255)
          aaaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑CC‑‑CCCCcctct                                                                                                                       <  2:186474/139‑2 (MQ=255)
                              cGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTg                                                                                                     <  1:499964/139‑1 (MQ=255)
                                 aacaaATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCg                                                                                                  >  2:1249383/1‑139 (MQ=255)
                                           gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                                                                        >  2:1604931/1‑139 (MQ=255)
                                           gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                                                                        >  2:350944/1‑139 (MQ=255)
                                            aCCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATAGCCGTTAAGCCA‑‑TCA‑‑‑CCCACCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                       >  2:325186/1‑139 (MQ=255)
                                                aCAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaa                                                                                   >  2:1246232/1‑139 (MQ=255)
                                                 cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑cc‑‑cccc                                                                                                                            >  1:1516314/1‑95 (MQ=255)
                                                 cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑cc‑‑cccc                                                                                                                            <  2:1516314/95‑1 (MQ=255)
                                                      cATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                                             >  1:606964/1‑139 (MQ=255)
                                                      cATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                                             >  2:14079/1‑139 (MQ=255)
                                                        ttACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACCCCC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTaa                                                                            <  1:1583504/139‑1 (MQ=255)
                                                        ttACACT‑GCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                          <  2:215224/139‑1 (MQ=255)
                                                         tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑cc‑‑cccc                                                                                                                            >  1:1311919/1‑87 (MQ=255)
                                                         tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑cc‑‑cccc                                                                                                                            <  2:1311919/87‑1 (MQ=255)
                                                                gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                          >  1:182191/1‑132 (MQ=255)
                                                                gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                          <  2:182191/132‑1 (MQ=255)
                                                                gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATa                                                                   >  2:629166/1‑139 (MQ=255)
                                                                gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATa                                                                   >  1:600393/1‑139 (MQ=255)
                                                                   aTTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                                                <  1:325186/139‑1 (MQ=255)
                                                                     ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAtt                                                              >  2:877894/1‑139 (MQ=255)
                                                                     ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAtt                                                              >  1:1361810/1‑139 (MQ=255)
                                                                     ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACAtt                                                              >  2:300920/1‑139 (MQ=255)
                                                                      tttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑CC‑‑CCCCcctctcagtg                                                                                                                  >  1:855318/1‑78 (MQ=255)
                                                                      tttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑CC‑‑CCCCcctctcagtg                                                                                                                  <  2:855318/84‑7 (MQ=255)
                                                                               cTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                                                    <  1:629166/139‑1 (MQ=255)
                                                                                  gACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGtt                                                 >  1:943737/1‑139 (MQ=255)
                                                                                       tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAtt                                                                               <  1:875333/104‑1 (MQ=255)
                                                                                       tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAtt                                                                               >  2:875333/1‑104 (MQ=255)
                                                                                       tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                            >  2:456429/1‑139 (MQ=255)
                                                                                       tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                            >  2:689830/1‑139 (MQ=255)
                                                                                        ctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                           >  1:186726/1‑139 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                                                <  2:1497091/114‑1 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                                                >  1:1497091/1‑114 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCg                                       >  2:656823/1‑139 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCg                                       >  1:817718/1‑139 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCg                                       >  2:1109749/1‑139 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCg                                       >  1:1272654/1‑139 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCg                                       >  1:1391074/1‑139 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCAC‑‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                      <  1:1374222/139‑1 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑c                                                                                                                                   <  1:1124842/47‑1 (MQ=255)
                                                                                            cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑c                                                                                                                                   >  2:1124842/1‑47 (MQ=255)
                                                                                             tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                         >  2:158799/1‑136 (MQ=255)
                                                                                             tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                         <  1:158799/136‑1 (MQ=255)
                                                                                               gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                    <  1:350944/139‑1 (MQ=255)
                                                                                               gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                    <  2:1391074/139‑1 (MQ=255)
                                                                                               gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                    <  1:456429/139‑1 (MQ=255)
                                                                                               gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                    <  2:817718/139‑1 (MQ=255)
                                                                                               gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑c                                                                                                                                   <  1:1590366/44‑1 (MQ=255)
                                                                                               gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑c                                                                                                                                   >  2:1590366/1‑44 (MQ=255)
                                                                                                tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTcagtgc                                                                                                                   <  2:514781/59‑1 (MQ=255)
                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTcagtgc                                                                                                                   <  2:426626/58‑1 (MQ=255)
                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTcagtgc                                                                                                                   >  1:426626/1‑58 (MQ=255)
                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                      <  2:339895/119‑1 (MQ=255)
                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                      >  1:339895/1‑119 (MQ=255)
                                                                                                 aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCa                                  >  2:780634/1‑139 (MQ=255)
                                                                                                  aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                 <  1:14079/139‑1 (MQ=255)
                                                                                                  aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                 <  1:780634/139‑1 (MQ=255)
                                                                                                  aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                 <  1:877894/139‑1 (MQ=255)
                                                                                                  aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                 <  2:1361810/139‑1 (MQ=255)
                                                                                                     aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCTCC‑‑‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAAACTGCCGTTACGGTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTctggttcttaaaca                              >  2:57665/1‑126 (MQ=255)
                                                                                                       gCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                 >  1:511581/1‑134 (MQ=255)
                                                                                                         tAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                         >  2:1201310/1‑124 (MQ=255)
                                                                                                         tAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                         <  1:1201310/124‑1 (MQ=255)
                                                                                                          aaaTAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                         <  2:1491523/123‑1 (MQ=255)
                                                                                                          aaaTAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                         >  1:1491523/1‑123 (MQ=255)
                                                                                                           aaTAGCAGAACATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTa                        >  1:1371736/1‑139 (MQ=255)
                                                                                                                   aaCATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                    <  1:1544226/119‑1 (MQ=255)
                                                                                                                   aaCATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                    >  2:1544226/1‑119 (MQ=255)
                                                                                                                     cATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                                                <  1:33066/89‑1 (MQ=255)
                                                                                                                     cATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                                                >  2:33066/1‑89 (MQ=255)
                                                                                                                     cATATCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGattatt              >  2:1225676/1‑139 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                    >  2:559719/1‑66 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                    <  1:559719/66‑1 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                          >  1:72131/1‑76 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                          <  2:72131/76‑1 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                                                         <  2:1260920/77‑1 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                                                         >  1:1260920/1‑77 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                      <  1:549660/96‑1 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                      >  2:549660/1‑96 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACa           >  2:298010/1‑139 (MQ=255)
                                                                                                                        atCGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACa           >  2:168709/1‑139 (MQ=255)
                                                                                                                          cGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                    <  1:352064/112‑1 (MQ=255)
                                                                                                                          cGCCGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                    >  2:352064/1‑112 (MQ=255)
                                                                                                                             cGTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGa      >  2:1374129/1‑139 (MQ=255)
                                                                                                                              gTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc     <  1:1374129/139‑1 (MQ=255)
                                                                                                                              gTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc     <  1:1109749/139‑1 (MQ=255)
                                                                                                                              gTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc     <  1:656823/139‑1 (MQ=255)
                                                                                                                              gTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc     <  1:298010/139‑1 (MQ=255)
                                                                                                                              gTTAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc     <  1:1249383/139‑1 (MQ=255)
                                                                                                                               ttAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                         >  2:622938/1‑102 (MQ=255)
                                                                                                                               ttAAGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                         <  1:622938/102‑1 (MQ=255)
                                                                                                                                 aaGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTGCAAAGGACTTg  >  2:1052579/1‑139 (MQ=255)
                                                                                                                                  aGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                      <  2:368151/102‑1 (MQ=255)
                                                                                                                                  aGCCA‑‑TCACC‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                      >  1:368151/1‑102 (MQ=255)
                                                                                                                                            ||                                                                                                                                   
AATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCA‑‑TCA‑‑‑CC‑‑CCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTG  >  NC_000913/1712375‑1712640

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