Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,712,518 (C)6→7 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.C98.9% 313.6 / ‑3.8 91intergenic (+354/‑257)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (46/44);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 109 total aligned reads displayed.

GAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  NC_000913/1712374‑1712645
                                                                                                                                           |                                                                                                                                      
gAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         <  2:798238/139‑1 (MQ=255)
  aTCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑cc                                                                                                                                      <  1:330752/139‑1 (MQ=255)
     aGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         <  1:645771/134‑1 (MQ=255)
     aGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         >  2:645771/1‑134 (MQ=255)
          aaaaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccccctct                                                                                                                              <  2:723360/139‑5 (MQ=255)
            aaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTca                                                                                                                            >  1:1813911/1‑139 (MQ=255)
            aaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                   <  2:128211/132‑1 (MQ=255)
            aaTCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                   >  1:128211/1‑132 (MQ=255)
               cATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTcagtg                                                                                                                         <  2:1813911/139‑1 (MQ=255)
                           cacacGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAAt                                                                                                             >  2:826012/1‑139 (MQ=255)
                            acacGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATc                                                                                                            <  2:883414/139‑1 (MQ=255)
                                            gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                            >  2:1700799/1‑139 (MQ=255)
                                            gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                            >  1:812509/1‑139 (MQ=255)
                                            gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                   >  2:312994/1‑100 (MQ=255)
                                            gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                   <  1:312994/100‑1 (MQ=255)
                                                  cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaaaa                                                                                      >  2:816277/1‑139 (MQ=255)
                                                         ttACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               <  2:812509/139‑1 (MQ=255)
                                                          tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                                                              >  2:1988913/1‑139 (MQ=255)
                                                          tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                                                              >  2:1911878/1‑139 (MQ=255)
                                                            cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                                                            >  1:1789754/1‑139 (MQ=255)
                                                            cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                                                            >  1:483428/1‑139 (MQ=255)
                                                            cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccccc                                                                                                                                 <  2:130445/86‑2 (MQ=255)
                                                            cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccccc                                                                                                                                 >  1:130445/1‑85 (MQ=255)
                                                              cTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                                                          <  1:1700799/139‑1 (MQ=255)
                                                               tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccccctct                                                                                                                              <  1:1996321/86‑5 (MQ=255)
                                                               tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccccctct                                                                                                                              >  2:1996321/1‑82 (MQ=255)
                                                                 tatatcagacataaaCTGGACTATCTCACTAGTAAAAACGATAAATACCAGAACCTATCGCCGTCAAGCCATCACCCCCCC‑TCTCAGTGCAGTCAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       <  1:1228979/127‑1 (MQ=255)
                                                                 gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  1:407466/1‑139 (MQ=255)
                                                                 gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  2:289001/1‑139 (MQ=255)
                                                                      ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         >  2:531217/1‑69 (MQ=255)
                                                                      ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         <  1:531217/69‑1 (MQ=255)
                                                                          cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCg                                                              <  2:726510/139‑1 (MQ=255)
                                                                           aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                                                             <  1:1945815/139‑1 (MQ=255)
                                                                                        tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGACTgtctcttata                                                >  2:1228979/1‑130 (MQ=255)
                                                                                          tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  2:407466/139‑1 (MQ=255)
                                                                                           cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTcccc                                             >  2:941215/1‑139 (MQ=255)
                                                                                             cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  1:1988913/139‑1 (MQ=255)
                                                                                             cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  2:495801/1‑139 (MQ=255)
                                                                                             cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                   >  1:431019/1‑51 (MQ=255)
                                                                                             cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                   <  2:431019/51‑1 (MQ=255)
                                                                                              tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  1:1528263/139‑1 (MQ=255)
                                                                                               aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  2:1757265/139‑1 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                                                          <  2:826053/103‑1 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                                                          >  1:826053/1‑103 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  1:130590/131‑1 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  1:795792/1‑139 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      <  2:795792/139‑1 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  2:478242/1‑139 (MQ=255)
                                                                                                   aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGt                                     >  2:310453/1‑139 (MQ=255)
                                                                                                      aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAtt                                                                                    <  1:1235476/89‑1 (MQ=255)
                                                                                                      aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAtt                                                                                    >  2:1235476/1‑89 (MQ=255)
                                                                                                      aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  2:274979/131‑1 (MQ=255)
                                                                                                      aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtat                                  >  2:1101478/1‑139 (MQ=255)
                                                                                                        gCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         <  1:186863/82‑1 (MQ=255)
                                                                                                        gCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         >  2:186863/1‑82 (MQ=255)
                                                                                                        gCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          >  1:354483/1‑129 (MQ=255)
                                                                                                        gCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  2:354483/129‑1 (MQ=255)
                                                                                                          tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                                >  2:8518/1‑121 (MQ=255)
                                                                                                          tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                                <  1:8518/121‑1 (MQ=255)
                                                                                                          tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  2:633685/128‑1 (MQ=255)
                                                                                                          tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         >  1:633685/1‑128 (MQ=255)
                                                                                                           aaaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         >  1:533055/1‑127 (MQ=255)
                                                                                                           aaaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  2:533055/127‑1 (MQ=255)
                                                                                                            ataagagacagatatCGCCGTTAAGCCATCA‑CCCCCCTTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  1:1309004/128‑1 (MQ=255)
                                                                                                            aaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAGTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  1:478242/139‑1 (MQ=255)
                                                                                                                 cAGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGat                       <  1:941215/139‑1 (MQ=255)
                                                                                                                  aGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  2:911321/115‑1 (MQ=255)
                                                                                                                  aGAACATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  1:911321/1‑115 (MQ=255)
                                                                                                                    aaCATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAGAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  1:2005333/131‑1 (MQ=255)
                                                                                                                    aaCATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAGAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            >  2:2005333/1‑131 (MQ=255)
                                                                                                                    aaCATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  1:1905212/118‑1 (MQ=255)
                                                                                                                    aaCATATCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         >  2:1905212/1‑118 (MQ=255)
                                                                                                                        tatCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         >  1:400957/1‑114 (MQ=255)
                                                                                                                        tatCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACa                <  1:826012/139‑1 (MQ=255)
                                                                                                                        tatCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACa                <  1:816277/139‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTaa                                                                                 >  1:223978/1‑73 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTaa                                                                                 <  2:223978/73‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGa                                                            >  2:96519/1‑94 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGa                                                            <  1:96519/94‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  1:101127/1‑108 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  2:101127/108‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  1:357026/113‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         >  2:357026/1‑113 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    >  1:1343853/1‑118 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    <  2:1343853/118‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  1:726703/126‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACaa               >  2:1704963/1‑139 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACa                >  2:67352/1‑138 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACa                <  1:67352/138‑1 (MQ=255)
                                                                                                                           cGCCGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAg             >  2:1302353/1‑139 (MQ=255)
                                                                                                                              cGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc          <  1:310453/139‑1 (MQ=255)
                                                                                                                              cGTTAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc          <  1:1911878/139‑1 (MQ=255)
                                                                                                                                ttAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    <  2:1306638/111‑1 (MQ=255)
                                                                                                                                ttAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    >  1:1306638/1‑111 (MQ=255)
                                                                                                                                ttAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACtt        >  2:1349051/1‑139 (MQ=255)
                                                                                                                                 tAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  1:1994401/1‑100 (MQ=255)
                                                                                                                                 tAAGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  2:1994401/100‑1 (MQ=255)
                                                                                                                                  aaGCCATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATGTGCCGTTACGGTTTTTGAAAAATTTACCGCTGGATAACATTTCCCGGAATAGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTAcctgtctcttat      >  2:1309004/1‑128 (MQ=255)
                                                                                                                                      cATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGa  <  1:289001/139‑1 (MQ=255)
                                                                                                                                      cATCACCCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGa  <  1:1101478/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                      
GAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCC‑TCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  NC_000913/1712374‑1712645

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