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.6% 235.5 / ‑3.9 69intergenic (+354/‑257)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (37/31);  total (38/31)
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.98e-01

GAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  NC_000913/1712374‑1712645
                                                                                                                                           |                                                                                                                                      
gAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         <  2:177126/139‑1 (MQ=255)
gAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         <  1:250653/139‑1 (MQ=255)
     aGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                                                                   <  1:285041/139‑1 (MQ=255)
     aGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                                                                   <  2:1471536/139‑1 (MQ=255)
                 tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         >  2:558654/1‑122 (MQ=255)
                 tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         <  1:558654/122‑1 (MQ=255)
                 tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcagtgca                                                                                                                       >  2:813486/1‑139 (MQ=255)
                                            gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         >  1:533833/1‑95 (MQ=255)
                                            gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         <  2:533833/95‑1 (MQ=255)
                                            gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                            >  2:626200/1‑139 (MQ=255)
                                            gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                            >  2:268282/1‑139 (MQ=255)
                                              ccTACAAATCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGa                                                                                          >  2:1214065/1‑139 (MQ=255)
                                               cTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         >  1:1257912/1‑139 (MQ=255)
                                               cTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         >  2:1451691/1‑139 (MQ=255)
                                                  cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaaaa                                                                                      >  1:785126/1‑139 (MQ=255)
                                                         ttACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               <  1:1424863/139‑1 (MQ=255)
                                                          tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               >  2:1099321/1‑138 (MQ=255)
                                                          tACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               <  1:1099321/138‑1 (MQ=255)
                                                            cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                                                            >  1:448540/1‑139 (MQ=255)
                                                               tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         >  2:1360002/1‑76 (MQ=255)
                                                               tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         <  1:1360002/76‑1 (MQ=255)
                                                               tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         >  2:412280/1‑76 (MQ=255)
                                                               tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                         <  1:412280/76‑1 (MQ=255)
                                                               tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcagtgc                                                                                                                        >  2:1417605/1‑92 (MQ=255)
                                                               tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcagtgc                                                                                                                        <  1:1417605/92‑1 (MQ=255)
                                                               tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGAt                                                                         >  1:1360803/1‑139 (MQ=255)
                                                                 gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTGAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  1:967045/1‑139 (MQ=255)
                                                                 gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  2:1431150/1‑139 (MQ=255)
                                                                                   gACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                                                     <  1:268282/139‑1 (MQ=255)
                                                                                        tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                                >  1:455209/1‑139 (MQ=255)
                                                                                               aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  1:1423970/139‑1 (MQ=255)
                                                                                               aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  2:455209/139‑1 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGc                                                                                                         <  1:574315/72‑1 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGc                                                                                                         >  2:574315/1‑72 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          >  1:612454/1‑135 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  2:612454/135‑1 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  2:910043/1‑139 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  1:1463495/1‑139 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      <  2:1360803/139‑1 (MQ=255)
                                                                                                  aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      <  2:785126/139‑1 (MQ=255)
                                                                                                    aaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    <  1:1431150/139‑1 (MQ=255)
                                                                                                      aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtat                                  >  2:467625/1‑139 (MQ=255)
                                                                                                        gCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                                                                   >  2:581380/1‑40 (MQ=255)
                                                                                                        gCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccc                                                                                                                                   <  1:581380/40‑1 (MQ=255)
                                                                                                          tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                                >  2:1143481/1‑121 (MQ=255)
                                                                                                          tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                                <  1:1143481/121‑1 (MQ=255)
                                                                                                          tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGt                              >  1:1215367/1‑139 (MQ=255)
                                                                                                           aaaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTaa                                                                                 <  1:645543/87‑1 (MQ=255)
                                                                                                           aaaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTaa                                                                                 >  2:645543/1‑87 (MQ=255)
                                                                                                            aaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  2:1215367/139‑1 (MQ=255)
                                                                                                            aaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  1:55867/139‑1 (MQ=255)
                                                                                                            aaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  1:1451691/139‑1 (MQ=255)
                                                                                                            aaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  2:425128/139‑1 (MQ=255)
                                                                                                              tAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACa                          <  2:1380312/139‑1 (MQ=255)
                                                                                                                  aGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGatt                      >  1:1289781/1‑139 (MQ=255)
                                                                                                                    aaCATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         >  2:475845/1‑70 (MQ=255)
                                                                                                                    aaCATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         <  1:475845/70‑1 (MQ=255)
                                                                                                                    aaCATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                           <  1:561699/100‑1 (MQ=255)
                                                                                                                    aaCATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                           >  2:561699/1‑100 (MQ=255)
                                                                                                                    aaCATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGattat                    >  1:298963/1‑139 (MQ=255)
                                                                                                                      cATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                                                   <  2:958397/106‑1 (MQ=255)
                                                                                                                      cATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                                                   >  1:958397/1‑106 (MQ=255)
                                                                                                                      cATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      <  1:1342950/119‑1 (MQ=255)
                                                                                                                      cATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  2:1342950/1‑119 (MQ=255)
                                                                                                                      cATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGa                        >  2:834911/1‑133 (MQ=255)
                                                                                                                      cATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGa                        <  1:834911/133‑1 (MQ=255)
                                                                                                                      cATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTCATATGTAACAGAttatta                  <  2:967045/139‑1 (MQ=255)
                                                                                                                        tatCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      <  2:465889/117‑1 (MQ=255)
                                                                                                                        tatCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  1:465889/1‑117 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTcc                                                                <  1:289152/90‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTcc                                                                >  2:289152/1‑90 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                                                             >  1:185261/1‑93 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                                                             <  2:185261/93‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  1:712929/1‑111 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  2:712929/111‑1 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACaa               >  1:126430/1‑139 (MQ=255)
                                                                                                                         atCGCCGTTAAGCCATAACCCCCCATCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAAAGATGGATAAAATTTCCCGGAATGGTTGAATTACCCGCCTAAGTTATAAGGAAAAGATTATTACaa               >  1:1096941/1‑139 (MQ=255)
                                                                                                                                ttAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         <  1:305115/58‑1 (MQ=255)
                                                                                                                                ttAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         >  2:305115/1‑58 (MQ=255)
                                                                                                                                ttAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTcc                                                                <  2:1311820/83‑1 (MQ=255)
                                                                                                                                ttAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTcc                                                                >  1:1311820/1‑83 (MQ=255)
                                                                                                                                  aaGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGt      >  2:7089/1‑139 (MQ=255)
                                                                                                                                      cATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGa  <  2:1463495/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                      
GAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  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: