Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,527 Δ1 bp intergenic (+36/+23) yjgN → / ← valS hypothetical protein/valine tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,478,5270G.99.5% 985.4 / ‑4.1 209intergenic (+36/+23)yjgN/valShypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (122/86);  total (122/87)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

AGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCA  >  NC_000913/4478437‑4478614
                                                                                          |                                                                                       
    ttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCCGGCTTCGGTTtcat                                                                                >  1:5754827/1‑93 (MQ=255)
     tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                                                               >  1:3252925/1‑93 (MQ=255)
      cttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatca                                                                              >  2:3141841/1‑93 (MQ=255)
      cttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatca                                                                              <  2:3467914/93‑1 (MQ=255)
       ttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCACTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCAc                                                                             >  1:3011201/1‑93 (MQ=255)
        tAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACt                                                                            <  2:1524975/93‑1 (MQ=255)
         aGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtg                                                                           >  2:5982889/1‑93 (MQ=255)
           ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                                                         >  1:6128033/1‑93 (MQ=255)
           ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                                                         <  2:1370140/93‑1 (MQ=255)
           ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                                                         >  1:2040309/1‑93 (MQ=255)
            cAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTgtgt                                                                        >  2:3902474/1‑93 (MQ=255)
              gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGttt                                                                      >  1:2152542/1‑93 (MQ=255)
              gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCCGCCCTCGTTTTCATCACTGTGttt                                                                      >  1:2367899/1‑93 (MQ=255)
                 tttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGGAAGCTGCGGCCTCCGTTTTTACCACCGTGctttta                                                                   >  1:80522/1‑91 (MQ=255)
                 tttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGACCGAGCATGCTCCGGCCCTCGTTTTCCTCACTGTGGTTTgt                                                                   >  2:42358/1‑92 (MQ=255)
                  ttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAt                                                                  <  1:5698015/93‑1 (MQ=255)
                   tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAtt                                                                 >  2:3286806/1‑93 (MQ=255)
                      cGGGTAGTTGTGCCTTCTTTCCCATTTCGGTAAGAGATAAAAAGGGCCGGAGCATGCTCTGGCCTTCGTTTTCATCACTGTGTTTTGATTACa                                                              <  2:4569903/93‑1 (MQ=255)
                       ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                                             <  2:5755625/93‑1 (MQ=255)
                       ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                                             >  2:5312388/1‑93 (MQ=255)
                        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                                            >  1:1978462/1‑93 (MQ=255)
                        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                                            <  2:2995086/93‑1 (MQ=255)
                        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                                            <  1:4010545/93‑1 (MQ=255)
                        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                                            >  2:1838190/1‑93 (MQ=255)
                         ggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcg                                                           >  2:5017749/1‑93 (MQ=255)
                          gAGTTATGCCTTCTTTCCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcgc                                                          <  2:5589828/93‑1 (MQ=255)
                          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcgc                                                          >  1:3856917/1‑93 (MQ=255)
                          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgcgc                                                          <  2:645875/93‑1 (MQ=255)
                            gTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcgg                                                        <  2:3588703/93‑1 (MQ=255)
                               aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                                                     <  2:2088679/93‑1 (MQ=255)
                               aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                                                     >  1:2413961/1‑93 (MQ=255)
                                 gctttcttTCCAATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                                                   <  1:5593551/90‑1 (MQ=255)
                                 gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                                                   <  1:5228294/93‑1 (MQ=255)
                                 gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                                                   >  1:6635873/1‑93 (MQ=255)
                                 gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                                                   >  1:1191425/1‑93 (MQ=255)
                                  ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTTATCACTGTGTTTTTATTACAGCGCGGCGATaa                                                  >  2:3749913/1‑93 (MQ=255)
                                  ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                                  >  1:6008575/1‑93 (MQ=255)
                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                                                 <  1:1930720/93‑1 (MQ=255)
                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCAACACTGTGGTTTGATTACAGCGCGGCGATCAc                                                 >  1:4827113/1‑93 (MQ=255)
                                      cttTCCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGc                                              <  1:4515789/93‑1 (MQ=255)
                                      cttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGc                                              <  2:862425/93‑1 (MQ=255)
                                       tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGcc                                             >  2:3042105/1‑93 (MQ=255)
                                       tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGcc                                             <  1:3305888/93‑1 (MQ=255)
                                       tttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGcc                                             <  2:1341262/93‑1 (MQ=255)
                                        ttACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                                            <  2:406884/93‑1 (MQ=255)
                                          aCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc                                          <  1:2967103/93‑1 (MQ=255)
                                           ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                                         >  2:6247578/1‑93 (MQ=255)
                                           ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                                         >  1:442541/1‑93 (MQ=255)
                                            cATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgctg                                        >  1:638211/1‑93 (MQ=255)
                                             aTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGt                                       >  1:1593881/1‑93 (MQ=255)
                                               ttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                                     >  1:5397302/1‑93 (MQ=255)
                                                tCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCa                                    >  2:793647/1‑93 (MQ=255)
                                                tCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCa                                    >  1:2189206/1‑93 (MQ=255)
                                                  tGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt                                  <  1:2721744/93‑1 (MQ=255)
                                                  tGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt                                  <  2:3110014/93‑1 (MQ=255)
                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCCATc                                 >  1:2855453/1‑93 (MQ=255)
                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc                                 >  1:2859519/1‑93 (MQ=255)
                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc                                 >  1:2794155/1‑93 (MQ=255)
                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc                                 <  2:2219905/93‑1 (MQ=255)
                                                    tAAGAGATAAAAAAGGCCGGGGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTTAATCa                                >  2:153727/1‑93 (MQ=255)
                                                       gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt                             <  1:5758346/93‑1 (MQ=255)
                                                            taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc                         >  1:4286074/2‑93 (MQ=255)
                                                            taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc                         >  2:5805931/2‑93 (MQ=255)
                                                             aaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATTACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgct                        >  1:3082017/1‑93 (MQ=255)
                                                              aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctt                       <  1:50493/93‑1 (MQ=255)
                                                               aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgcttt                      >  1:3870789/1‑93 (MQ=255)
                                                               aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgcttt                      >  2:1450116/1‑93 (MQ=255)
                                                               aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgcttt                      <  1:3396744/93‑1 (MQ=255)
                                                               aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttacgcttt                      <  1:5520622/93‑1 (MQ=255)
                                                                 aaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAATAGCTTGCTGTTCAATTAGtttcgctttcg                    >  1:641204/1‑93 (MQ=255)
                                                                  aGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttcgc                   >  1:1016365/1‑93 (MQ=255)
                                                                  aGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttcgc                   >  2:2349134/1‑93 (MQ=255)
                                                                   ggCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTttcgctttcgct                  <  2:1640903/93‑1 (MQ=255)
                                                                   ggCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTttcgctttcgct                  >  1:1706936/1‑93 (MQ=255)
                                                                   ggCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTttcgctttcgct                  >  1:2774090/1‑93 (MQ=255)
                                                                   ggCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTttcgctttcgct                  >  2:3988872/1‑93 (MQ=255)
                                                                   ggCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTttcgctttcgct                  <  1:4979616/93‑1 (MQ=255)
                                                                       cggagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTcc               <  2:682083/89‑1 (MQ=255)
                                                                        ggagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCg              >  2:1588815/4‑93 (MQ=255)
                                                                        ggagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCg              >  2:1814970/4‑93 (MQ=255)
                                                                        ggagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCg              >  1:1331985/4‑93 (MQ=255)
                                                                        ggagcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCg              <  2:3929901/90‑1 (MQ=255)
                                                                          agcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCa            >  2:2147238/2‑93 (MQ=255)
                                                                           gcatgcTCCGTCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCAt           <  2:5956274/93‑1 (MQ=255)
                                                                           gcatgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCAt           <  2:6039449/93‑1 (MQ=255)
                                                                            catgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATa          <  1:754365/93‑1 (MQ=255)
                                                                            catgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATa          >  1:5152983/1‑93 (MQ=255)
                                                                            catgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATa          >  2:5418614/1‑93 (MQ=255)
                                                                             atgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAg         >  2:5740871/1‑93 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGTATAGc        <  1:4327774/93‑1 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        >  2:1194320/1‑93 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        <  2:2277905/93‑1 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        >  1:696874/1‑93 (MQ=255)
                                                                              tgcTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGc        <  1:1040616/93‑1 (MQ=255)
                                                                                 tCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCtt     >  2:4314042/1‑93 (MQ=255)
                                                                                  ccGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTc    >  2:1041655/1‑93 (MQ=255)
                                                                                  ccGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTc    >  1:3581774/1‑93 (MQ=255)
                                                                                   cGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTcc   >  2:124373/1‑93 (MQ=255)
                                                                                   cGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTcc   <  2:635586/93‑1 (MQ=255)
                                                                                   cGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTcc   <  2:2008/93‑1 (MQ=255)
                                                                                          |                                                                                       
AGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCA  >  NC_000913/4478437‑4478614

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