Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,710,773 A→G intergenic (+591/‑20) nth → / → ydgR endonuclease III/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,710,7730AG99.6% 832.7 / NA 255intergenic (+591/‑20)nth/ydgRendonuclease III/putative transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (137/117);  total (137/118)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 263 total aligned reads displayed.

TGAACATTACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCT  >  NC_000913/1710683‑1710862
                                                                                          |                                                                                         
tGAACATTACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTg                                                                                         <  2:717868/93‑1 (MQ=255)
tGAACATTACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTg                                                                                         >  2:6193416/1‑93 (MQ=255)
  aaCATTACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTa                                                                                       >  2:799046/1‑93 (MQ=255)
     aTTACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGAACCCCCGTTAATATGGGGTGAaaaa                                                                                    >  2:2960918/1‑93 (MQ=255)
       tACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTaaaaaa                                                                                  <  2:2061115/93‑1 (MQ=255)
       tACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTaaaaaa                                                                                  >  2:1113766/1‑93 (MQ=255)
        aCCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTaaaaaaa                                                                                 <  1:6158337/93‑1 (MQ=255)
         ccTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAag                                                                                >  2:2869792/1‑93 (MQ=255)
          cTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAaga                                                                               >  1:1091973/1‑93 (MQ=255)
          cTGCCGTTTCCCCTCCCACTATAACAATGGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAaga                                                                               <  1:1445946/93‑1 (MQ=255)
           tGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAagag                                                                              >  1:4357374/1‑93 (MQ=255)
           tGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAagag                                                                              <  2:2698117/93‑1 (MQ=255)
           cgCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAagag                                                                              <  2:1866936/92‑1 (MQ=255)
             ccGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGt                                                                            <  1:6466904/93‑1 (MQ=255)
             ccGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGt                                                                            <  1:320908/93‑1 (MQ=255)
              cGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTa                                                                           <  2:3805112/93‑1 (MQ=255)
               gTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaa                                                                          <  1:6117988/93‑1 (MQ=255)
               gTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaa                                                                          >  1:3320278/1‑93 (MQ=255)
                tttCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaaa                                                                         >  2:3171303/1‑93 (MQ=255)
                tttCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaaa                                                                         >  2:6353420/1‑93 (MQ=255)
                 ttCCCCTCCCACTATAGCAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaaaa                                                                        <  1:2517954/93‑1 (MQ=255)
                 ttCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaaaa                                                                        <  1:3901465/93‑1 (MQ=255)
                 ttCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaaaa                                                                        <  2:4939770/93‑1 (MQ=255)
                 ttCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaaaa                                                                        >  1:4597206/1‑93 (MQ=255)
                 ttCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaaaa                                                                        <  2:5987710/93‑1 (MQ=255)
                 ttCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTaaaa                                                                        <  1:2908538/93‑1 (MQ=255)
                  tcgccTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCGCCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAg                                                                       <  2:2182324/90‑1 (MQ=255)
                   ccccTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAgt                                                                      >  1:772020/1‑93 (MQ=255)
                     ccTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAgtgt                                                                    >  1:1000833/1‑93 (MQ=255)
                     ccTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAgtgt                                                                    >  1:5716376/1‑93 (MQ=255)
                     ccTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAgtgt                                                                    >  1:3145667/1‑93 (MQ=255)
                      cTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTc                                                                   <  2:3081219/93‑1 (MQ=255)
                       tCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTcc                                                                  >  2:4536635/1‑93 (MQ=255)
                         ccACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCCCCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCAc                                                                <  2:4300353/93‑1 (MQ=255)
                         ccACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCAc                                                                <  2:4149630/93‑1 (MQ=255)
                         ccACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCAc                                                                >  1:6725614/1‑93 (MQ=255)
                           aCTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTg                                                              >  1:22626/1‑93 (MQ=255)
                              ataACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaa                                                           >  1:5634907/1‑93 (MQ=255)
                               taACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaa                                                          <  1:5575117/93‑1 (MQ=255)
                                aaCAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaac                                                         >  2:2346736/1‑93 (MQ=255)
                                 aCAATTGCGCGTATGTTTTTTATACATAACGCGATAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaacc                                                        >  2:3423618/1‑93 (MQ=255)
                                 aCAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaacc                                                        <  1:1571985/93‑1 (MQ=255)
                                   aaTTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaaccaa                                                      <  2:5052450/93‑1 (MQ=255)
                                    aTTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaaccaaa                                                     <  2:5615041/93‑1 (MQ=255)
                                     ttGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaaccaaaa                                                    <  1:2940231/93‑1 (MQ=255)
                                      tGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaaccaaaaa                                                   >  2:3690460/1‑93 (MQ=255)
                                       gcgcGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaaccaaaaac                                                  <  2:2598930/93‑1 (MQ=255)
                                        cgcgTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaaccaaaaacc                                                 <  1:3917325/93‑1 (MQ=255)
                                         gcgTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaaccaaaaacca                                                <  2:4374031/93‑1 (MQ=255)
                                          cgTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaaccaaaaaccaa                                               <  2:3779419/93‑1 (MQ=255)
                                          cgTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCaaaccaaaaaccaa                                               <  2:1169922/93‑1 (MQ=255)
                                           gTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAAc                                              <  2:1628026/93‑1 (MQ=255)
                                           gTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAAc                                              <  2:1805431/93‑1 (MQ=255)
                                           gTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAAc                                              >  2:5717787/1‑93 (MQ=255)
                                           gTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAAc                                              <  1:5042913/93‑1 (MQ=255)
                                            tATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACt                                             <  1:1895385/93‑1 (MQ=255)
                                              tGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGa                                           >  2:6414958/1‑93 (MQ=255)
                                              tGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGa                                           >  1:2720118/1‑93 (MQ=255)
                                                ttttttATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGaaa                                         <  2:1706806/93‑1 (MQ=255)
                                                 tttttATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAg                                        <  2:483891/93‑1 (MQ=255)
                                                  ttttATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGc                                       <  2:3932179/93‑1 (MQ=255)
                                                    ttATACATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGt                                     >  1:4751707/1‑93 (MQ=255)
                                                     tataCATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTc                                    >  2:3613769/1‑93 (MQ=255)
                                                      ataCATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCa                                   <  1:1795254/93‑1 (MQ=255)
                                                      ataCATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCa                                   <  1:5970314/93‑1 (MQ=255)
                                                       taCATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAg                                  >  2:271542/1‑93 (MQ=255)
                                                        aCATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGt                                 >  1:4877675/1‑93 (MQ=255)
                                                         cATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGtt                                <  2:6284772/93‑1 (MQ=255)
                                                         cATAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGtt                                <  1:5026298/93‑1 (MQ=255)
                                                          aTAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGttt                               >  2:5129608/1‑93 (MQ=255)
                                                          aTAACGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGttt                               <  2:2593325/93‑1 (MQ=255)
                                                            aaCGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGa                             <  1:4281063/93‑1 (MQ=255)
                                                            aaCGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGa                             >  1:3050292/1‑93 (MQ=255)
                                                             aCGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGaa                            <  2:4413020/93‑1 (MQ=255)
                                                             aCGCGAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGaa                            >  2:3021315/1‑93 (MQ=255)
                                                               gcgAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACg                          >  2:3979014/1‑93 (MQ=255)
                                                                cgAGAAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGc                         <  2:4026775/93‑1 (MQ=255)
                                                                   gaAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCttt                      >  2:5387931/1‑93 (MQ=255)
                                                                   gaAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCttt                      >  1:4026519/1‑93 (MQ=255)
                                                                   gaAAGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCttt                      >  2:2678225/1‑93 (MQ=255)
                                                                    aaaGCACCCCCGTTAATATGGGGAGTAAAAAAAGAGGTAAAAGTGTACACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTc                     >  1:5373322/1‑93 (MQ=255)
                                                                     aaGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCa                    <  1:4813664/93‑1 (MQ=255)
                                                                      aGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCaa                   >  1:2873290/1‑93 (MQ=255)
                                                                      aGCACCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCaa                   >  2:1404312/1‑93 (MQ=255)
                                                                         aCCCCCGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAaaca                >  2:6452304/1‑93 (MQ=255)
                                                                          cccccGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACAGCAAATCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAaacaa               >  2:6722953/1‑93 (MQ=255)
                                                                           ccccGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAaacaac              <  1:6641887/93‑1 (MQ=255)
                                                                           ccccGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAaacaac              >  2:1533577/1‑93 (MQ=255)
                                                                           ccccGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAaacaac              >  2:6525930/1‑93 (MQ=255)
                                                                            cccGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAAcc             >  1:255963/1‑93 (MQ=255)
                                                                            cccGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAAcc             >  2:6267506/1‑93 (MQ=255)
                                                                             ccGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCg            <  2:5199200/93‑1 (MQ=255)
                                                                              cGTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGa           <  1:964287/93‑1 (MQ=255)
                                                                               gTTAATATGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGaa          <  1:450344/93‑1 (MQ=255)
                                                                                ttAATATGGTGTGTAAAAAAAGAGGTAAAACTGTCCGCTGCAAACCAAAACCCAACTGAAACCGTCAGTTTGAACGCTTTCAAACAACCGAAg         <  1:5803403/93‑1 (MQ=255)
                                                                                    tatGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGt     >  2:174036/1‑93 (MQ=255)
                                                                                    tatGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGt     >  1:6405065/1‑93 (MQ=255)
                                                                                     atGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGtt    >  2:3582163/1‑93 (MQ=255)
                                                                                      tGGGGTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTc   >  2:1981343/1‑93 (MQ=255)
                                                                                       ggggTGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCt  <  2:5758920/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGAACATTACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCT  >  NC_000913/1710683‑1710862

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