Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,696,452 (A)5→4 intergenic (‑381/+10) uidA ← / ← uidR beta‑D‑glucuronidase/transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,696,4480A.97.4% 300.6 / 2.0 76intergenic (‑377/+14)uidA/uidRbeta‑D‑glucuronidase/transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (37/37);  total (38/38)
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.97e-01

ATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGCCAGTACACCGTAGGTCAGATCCAGCAGTA  >  NC_000913/1696323‑1696580
                                                                                                                             |                                                                                                                                    
aTAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAAt                                                                                                                        >  1:1366463/1‑139 (MQ=255)
                  tGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATa                                                                                                      >  1:407065/1‑139 (MQ=255)
                  tGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATa                                                                                                      >  2:445215/1‑139 (MQ=255)
                        cAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCa                                                                                                >  2:114957/1‑139 (MQ=255)
                           tGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAatc                                                                                             >  2:1507599/1‑139 (MQ=255)
                             aaTATTTTGGTGGGAATTCTCCGAAGCCGGGAAGTATAATTTTTGCACAGTGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGaa                                                                                                                                     <  2:1063843/98‑1 (MQ=255)
                             aaTATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGaa                                                                                                                                     >  1:1063843/1‑98 (MQ=255)
                             aaTATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATa                                                                                                      >  2:267322/1‑128 (MQ=255)
                             aaTATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATa                                                                                                      <  1:267322/128‑1 (MQ=255)
                             aaTATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAatcat                                                                                           <  2:1354724/139‑1 (MQ=255)
                             aaTATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAatcat                                                                                           >  2:347942/1‑139 (MQ=255)
                                 ttttGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATcgcg                                                                                       <  1:347942/139‑1 (MQ=255)
                                   ttGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATcgcgcg                                                                                     <  1:1507599/139‑1 (MQ=255)
                                    tGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGt                                                                                    <  2:407065/139‑1 (MQ=255)
                                     ggggggAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAgg                                                                                                                    >  1:139441/1‑106 (MQ=255)
                                     ggggggAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAgg                                                                                                                    <  2:139441/106‑1 (MQ=255)
                                       ggggAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGa                                                                                                                                      >  2:404493/1‑87 (MQ=255)
                                       ggggAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGa                                                                                                                                      <  1:404493/87‑1 (MQ=255)
                                        gggAATTGTCCGAAGGCGGGAAGAATAATTTTTGCACAGAGAGATTACGATATGAAGCGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATc                                                                                <  1:1206298/139‑1 (MQ=255)
                                           aaTTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAgaa                                                                                                                   >  1:233578/1‑99 (MQ=255)
                                           aaTTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAgaa                                                                                                                   <  2:233578/101‑3 (MQ=255)
                                                tcCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCa                                                                        >  2:1259506/1‑139 (MQ=255)
                                                  cGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGa                                                                                                                   >  1:18305/1‑94 (MQ=255)
                                                  cGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGa                                                                                                                   <  2:18305/94‑1 (MQ=255)
                                                  cGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCaaa                                                                      >  2:480137/1‑139 (MQ=255)
                                                          ggAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGaaaagaga                                                                                                                               <  2:1358313/75‑5 (MQ=255)
                                                          ggAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGaaaagaga                                                                                                                               >  1:1358313/1‑71 (MQ=255)
                                                          ggAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATa                                                                                                      >  2:893608/1‑99 (MQ=255)
                                                          ggAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATa                                                                                                      <  1:893608/99‑1 (MQ=255)
                                                          ggAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCtt                                                              >  2:1381585/1‑139 (MQ=255)
                                                                aTAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCAcgcg                                                        <  2:1366463/139‑1 (MQ=255)
                                                                 tAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAacga                                                                   <  2:1411570/127‑1 (MQ=255)
                                                                 tAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAacga                                                                   >  1:1411570/1‑127 (MQ=255)
                                                                 tAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACgcgc                                                       >  1:310952/1‑139 (MQ=255)
                                                                  aaTTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGcc                                                      <  2:310952/139‑1 (MQ=255)
                                                                  aaTTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGcc                                                      <  1:1259506/139‑1 (MQ=255)
                                                                         gCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAatcat                                                                                           <  1:729764/95‑1 (MQ=255)
                                                                         gCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAatcat                                                                                           >  2:729764/1‑95 (MQ=255)
                                                                         gCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATcc                                                                               >  2:842149/1‑107 (MQ=255)
                                                                         gCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATcc                                                                               <  1:842149/107‑1 (MQ=255)
                                                                              gagaAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGa                                                                                                        <  1:917654/77‑1 (MQ=255)
                                                                              gagaAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGa                                                                                                        >  2:917654/1‑77 (MQ=255)
                                                                              gagaAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCt                                          >  1:258180/1‑139 (MQ=255)
                                                                                 aaaTTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATg                                                                                                                 <  2:924354/65‑1 (MQ=255)
                                                                                 aaaTTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATg                                                                                                                 >  1:924354/1‑65 (MQ=255)
                                                                                    ttACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGAGGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGata                                    <  1:428643/139‑1 (MQ=255)
                                                                                         atatGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGc                               <  2:258180/139‑1 (MQ=255)
                                                                                         atatGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGc                               <  1:480137/139‑1 (MQ=255)
                                                                                           atGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGAt                                                                                                                  >  1:5903/1‑54 (MQ=255)
                                                                                           atGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGAt                                                                                                                  <  2:5903/54‑1 (MQ=255)
                                                                                           atGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGc                                                                                                                <  2:1043410/56‑1 (MQ=255)
                                                                                           atGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGc                                                                                                                >  1:1043410/1‑56 (MQ=255)
                                                                                              aaaCGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGCCAGTa                          <  1:1381585/139‑1 (MQ=255)
                                                                                              aaaCGCGGGCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGCCAGTa                          <  1:445215/139‑1 (MQ=255)
                                                                                                    gggCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGt                                                                                                             >  2:1355566/1‑50 (MQ=255)
                                                                                                    gggCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGt                                                                                                             <  1:1355566/50‑1 (MQ=255)
                                                                                                    gggCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAatcat                                                                                           >  1:642377/1‑68 (MQ=255)
                                                                                                    gggCACAGTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAatcat                                                                                           <  2:642377/68‑1 (MQ=255)
                                                                                                           gTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGc                                                                                                                >  2:464657/1‑40 (MQ=255)
                                                                                                           gTTATCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGc                                                                                                                <  1:464657/40‑1 (MQ=38)
                                                                                                              aTCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATcg                                                                                         >  2:234009/1‑60 (MQ=255)
                                                                                                              aTCACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATcg                                                                                         <  1:234009/60‑1 (MQ=255)
                                                                                                                cACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCa                                                            <  1:58116/87‑1 (MQ=255)
                                                                                                                cACCAGCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCa                                                            >  2:58116/1‑87 (MQ=255)
                                                                                                                     gCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCtt                                               <  1:783818/96‑1 (MQ=255)
                                                                                                                     gCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCtt                                               >  2:783818/1‑96 (MQ=255)
                                                                                                                     gCCCGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGCCAGTACACCGTAGGTCAGATCCAGCAGt   >  1:200386/1‑139 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAatcat                                                                                           <  2:1593785/50‑1 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAatcat                                                                                           >  1:1593785/1‑50 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCAcgc                                                         >  2:1379845/1‑84 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCAcgc                                                         <  1:1379845/84‑1 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGcc                                                      <  1:1208510/87‑1 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGcc                                                      >  2:1208510/1‑87 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCtt                                         >  1:367059/1‑100 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCtt                                         <  2:367059/100‑1 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGata                                    >  1:415777/1‑105 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGata                                    >  1:1168781/1‑105 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGata                                    <  2:415777/105‑1 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGata                                    <  2:1168781/105‑1 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGCCAGTACACCGTAGGTc               <  2:529934/126‑1 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGCCAGTACACCGTAGGTc               >  1:529934/1‑126 (MQ=255)
                                                                                                                      cccGCCGAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGCCAGTACACCGTAGGTCAGATCCAGCAGTa  >  1:986305/1‑139 (MQ=255)
                                                                                                                             |                                                                                                                                    
ATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTGAGCAAACGACGCTTCACGCGCCAGGTCTTCCGCTTCGATATCGGCCAGTACACCGTAGGTCAGATCCAGCAGTA  >  NC_000913/1696323‑1696580

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