Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,503 (A)7→6 intergenic (+12/+47) 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,4970A.97.9% 854.0 / 10.5 186intergenic (+6/+53)yjgN/valShypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base . (105/77);  total (109/77)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 216 total aligned reads displayed.

CCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTT  >  NC_000913/4478405‑4478587
                                                                                            |                                                                                          
ccttTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATa                                                                                            <  2:3383186/93‑1 (MQ=255)
 cttTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaa                                                                                           >  1:5833399/1‑93 (MQ=255)
   ttAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaa                                                                                         <  1:6190629/93‑1 (MQ=255)
    tAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaa                                                                                        >  1:4760321/1‑93 (MQ=255)
    tAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaa                                                                                        >  2:4778725/1‑93 (MQ=255)
    tAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaa                                                                                        >  1:2138156/1‑93 (MQ=255)
     aaTCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaa                                                                                       <  1:3281192/93‑1 (MQ=255)
     aaTCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaa                                                                                       >  2:1642108/1‑93 (MQ=255)
      aTCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaag                                                                                      <  2:2169406/93‑2 (MQ=255)
       tCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaagg                                                                                     <  2:3059620/93‑3 (MQ=255)
       tCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaagg                                                                                     <  2:4360469/93‑3 (MQ=255)
        cAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaaggc                                                                                    >  2:3105723/1‑90 (MQ=255)
          aTAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCg                                                                                 >  1:255285/1‑93 (MQ=255)
           tAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCgg                                                                                <  1:4997735/93‑1 (MQ=255)
             aagaagAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGGCCgag                                                                               >  2:3354756/1‑91 (MQ=255)
              agaagaACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagc                                                                             >  1:2505288/1‑91 (MQ=255)
              agaagaACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagc                                                                             >  1:5438866/1‑91 (MQ=255)
              agaagaACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagc                                                                             <  1:941209/93‑3 (MQ=255)
               gaagaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagca                                                                            >  1:1280794/1‑90 (MQ=255)
               gaagaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagca                                                                            <  1:1182003/93‑4 (MQ=255)
               gaagaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGagca                                                                            >  2:1569301/1‑90 (MQ=255)
                aagaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGGCgggggct                                                                            >  1:5543843/1‑89 (MQ=255)
                  gaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGgagcatgc                                                                        <  1:3185273/93‑1 (MQ=255)
                  gaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGgagcatgc                                                                        >  1:3765243/1‑93 (MQ=255)
                   aacaacCAGATAAAAGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCt                                                                       <  1:2670553/93‑1 (MQ=255)
                    acaacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTc                                                                      >  1:4162357/1‑93 (MQ=255)
                    acaacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTc                                                                      >  1:4726529/1‑93 (MQ=255)
                      aacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCg                                                                    <  2:3186823/93‑1 (MQ=255)
                      aacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCg                                                                    <  2:1171664/93‑1 (MQ=255)
                      aacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTAGGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCg                                                                    >  2:4536023/1‑93 (MQ=255)
                       acCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCgg                                                                   >  1:4200337/1‑93 (MQ=255)
                       acCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCgg                                                                   >  1:2825899/1‑93 (MQ=255)
                       acCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCgg                                                                   >  2:5090452/1‑93 (MQ=255)
                        ccAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGc                                                                  >  1:4064977/1‑93 (MQ=255)
                        ccAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGc                                                                  >  2:5288469/1‑93 (MQ=255)
                         cAGATAAAGTCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGcc                                                                 <  2:2517187/93‑1 (MQ=255)
                         cAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAACAAAGGCCGGAGCATGCTCCGGcc                                                                 >  1:396957/1‑93 (MQ=255)
                         cAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGcc                                                                 <  1:4100306/93‑1 (MQ=255)
                           gATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                               <  1:4525504/93‑1 (MQ=255)
                           gATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                               <  1:773641/93‑1 (MQ=255)
                           gATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCATTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                               <  2:3563906/93‑1 (MQ=255)
                            aTAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttc                                                              >  2:2712131/1‑92 (MQ=255)
                            aTAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttc                                                              >  2:3894649/1‑92 (MQ=255)
                             tAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcg                                                             >  2:5774671/1‑91 (MQ=255)
                              aaaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcgt                                                            >  1:1477508/1‑90 (MQ=255)
                              aaaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcgt                                                            >  2:5232858/1‑90 (MQ=255)
                                aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGttt                                                         >  1:1722957/1‑93 (MQ=255)
                                 ggCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGtttt                                                        >  1:2475672/1‑93 (MQ=255)
                                    ttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcat                                                     >  1:3524990/1‑93 (MQ=255)
                                    ttcttAGCCAGTATGTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCGTCGTTTtcat                                                     >  1:3012581/1‑93 (MQ=255)
                                     tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                                    <  1:3863944/93‑1 (MQ=255)
                                     tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                                    >  1:1679930/1‑93 (MQ=255)
                                       ttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGGTCCGGCCTTCGTTTTTATCAc                                                  >  1:160400/1‑93 (MQ=255)
                                        tAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACt                                                 >  2:5508457/1‑93 (MQ=255)
                                         aGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtg                                                >  2:4895512/1‑93 (MQ=255)
                                          gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                               >  2:5652338/1‑93 (MQ=255)
                                          gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                               >  2:4195419/1‑93 (MQ=255)
                                          gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                               >  2:2660758/1‑93 (MQ=255)
                                           ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgtg                                              <  2:3230335/93‑1 (MQ=255)
                                             aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                            >  2:5867338/1‑93 (MQ=255)
                                               tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                          >  1:5167038/1‑93 (MQ=255)
                                               tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                          <  1:4467131/93‑1 (MQ=255)
                                                aTTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTg                                         >  2:3011556/1‑93 (MQ=255)
                                                 tttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGa                                        >  2:553917/1‑93 (MQ=255)
                                                 tttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGa                                        <  1:1954871/93‑1 (MQ=255)
                                                   tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGAtt                                      >  2:554861/1‑93 (MQ=255)
                                                    cACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTa                                     >  2:5949496/1‑93 (MQ=255)
                                                     aCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTTATCAATGTGGTTttattaa                                    >  1:2482315/1‑92 (MQ=255)
                                                     aCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTAc                                    >  1:3938084/1‑93 (MQ=255)
                                                     aCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTAc                                    <  2:3151073/93‑1 (MQ=255)
                                                      cGGGGAGTTATGGCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATCCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACa                                   <  2:1267343/93‑1 (MQ=255)
                                                       ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                  <  2:1753290/93‑1 (MQ=255)
                                                        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGACAGGAGCATGCTCCAGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                 >  2:3876547/1‑93 (MQ=255)
                                                           aGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGcgcg                              <  2:4635917/93‑1 (MQ=255)
                                                               aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                          >  2:926696/1‑93 (MQ=255)
                                                               aTGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                          <  2:4588379/93‑1 (MQ=255)
                                                                 gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                        >  1:4945172/1‑93 (MQ=255)
                                                                 gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                        >  2:6157573/1‑93 (MQ=255)
                                                                 gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                        >  2:1897732/1‑93 (MQ=255)
                                                                 gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                        <  2:5244425/93‑1 (MQ=255)
                                                                  ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                       <  2:5372558/93‑1 (MQ=255)
                                                                   cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGGGCCTGCTCCCGCCTTCGTTTTTATCCCTGTGTTTTTATTACAGCGCGGCGATaaa                      >  2:71641/1‑92 (MQ=255)
                                                                      cttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGc                   <  2:6083989/93‑1 (MQ=255)
                                                                          acaaTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc               <  2:2389685/90‑1 (MQ=255)
                                                                            cATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgctg             <  2:143009/93‑1 (MQ=255)
                                                                             aTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGt            <  1:657430/93‑1 (MQ=255)
                                                                             aTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGt            <  1:3320716/93‑1 (MQ=255)
                                                                             aTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGt            >  2:1684187/1‑93 (MQ=255)
                                                                              tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt           >  1:5357865/1‑93 (MQ=255)
                                                                              tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt           <  1:4629687/93‑1 (MQ=255)
                                                                                   gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc      <  2:6471787/93‑1 (MQ=255)
                                                                                     aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGACTGCTGTTCAATCAg    >  2:2779075/1‑93 (MQ=255)
                                                                                      agagATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGt   >  1:5900462/1‑93 (MQ=255)
                                                                                      agagATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGt   >  1:6164919/1‑93 (MQ=255)
                                                                                       gagaTAAAAAAGGCGGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt  <  1:5217431/93‑1 (MQ=255)
                                                                                       gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt  <  1:5655282/93‑1 (MQ=255)
                                                                                       gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt  >  1:6237265/1‑93 (MQ=255)
                                                                                       gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt  >  1:3065912/1‑93 (MQ=255)
                                                                                       gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt  <  1:2754606/93‑1 (MQ=255)
                                                                                       gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt  <  2:852344/93‑1 (MQ=255)
                                                                                            |                                                                                          
CCTTTAATCAATAAAGAAGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTT  >  NC_000913/4478405‑4478587

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