Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,511 Δ1 bp intergenic (+20/+39) 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,5110T.100.0% 915.8 / NA 200intergenic (+20/+39)yjgN/valShypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (109/91);  total (109/91)

Only 100 of 219 total aligned reads displayed.

AGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCT  >  NC_000913/4478422‑4478598
                                                                                         |                                                                                       
agaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGgagcatg                                                                                    <  2:3121455/93‑1 (MQ=255)
agaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTGCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGgagcatg                                                                                    <  2:5814335/93‑1 (MQ=255)
 gaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGgagcatgc                                                                                   >  2:5638347/1‑93 (MQ=255)
    caacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTcc                                                                                >  1:6305223/1‑93 (MQ=255)
    caacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTcc                                                                                <  2:2050550/93‑1 (MQ=255)
    caacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTcc                                                                                <  2:1298469/93‑1 (MQ=255)
    caacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTcc                                                                                <  1:5704434/93‑1 (MQ=255)
     aacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCg                                                                               >  1:4257011/1‑93 (MQ=255)
     aacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCg                                                                               >  1:5811288/1‑93 (MQ=255)
     aacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCg                                                                               >  1:5733990/1‑93 (MQ=255)
       ccAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGc                                                                             >  1:6289206/1‑93 (MQ=255)
       ccAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGc                                                                             >  1:195586/1‑93 (MQ=255)
        cAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGcc                                                                            >  1:87073/1‑93 (MQ=255)
        cAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGcc                                                                            <  2:1636209/93‑1 (MQ=255)
        cAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGcc                                                                            >  2:4846319/1‑93 (MQ=255)
          gATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                                          >  2:3435890/1‑93 (MQ=255)
          gATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                                          <  2:5335028/93‑1 (MQ=255)
          gATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCtt                                                                          >  1:4008917/1‑93 (MQ=255)
            tAAAGCCTTCTTAGCCAGTATTTCACGGTGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcg                                                                        <  1:2399275/93‑3 (MQ=255)
             aaaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcgt                                                                       >  1:1626280/1‑90 (MQ=255)
              aaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcgtt                                                                      >  2:6248537/1‑89 (MQ=255)
              aaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcgtt                                                                      <  1:6112021/93‑5 (MQ=255)
              aaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCttcgtt                                                                      <  1:4420222/93‑5 (MQ=255)
               aGGCTTCTTAGCCAGTATTTCACGGGTAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGttt                                                                    <  2:3825628/93‑1 (MQ=255)
               aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGttt                                                                    >  1:4813808/1‑93 (MQ=255)
               aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGttt                                                                    <  2:4646110/93‑1 (MQ=255)
                ggCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTTTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGtttt                                                                   >  2:2628728/1‑93 (MQ=255)
                   ttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcat                                                                >  2:3608844/1‑93 (MQ=255)
                    tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTtcatc                                                               >  2:4783841/1‑93 (MQ=255)
                        aGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtg                                                           >  1:3534506/1‑93 (MQ=255)
                         gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                                          >  1:1381202/1‑93 (MQ=255)
                         gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACtgt                                                          <  1:2962062/93‑1 (MQ=255)
                            aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                                       >  1:1666099/1‑93 (MQ=255)
                            aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                                       <  1:5126311/93‑1 (MQ=255)
                            aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                                       <  2:4952738/93‑1 (MQ=255)
                            aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtt                                                       >  2:1576657/1‑93 (MQ=255)
                             gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGGGATAAAAAAGGGCGGAGCATGCGCCGGACCTCGTTTTCCTCACTctgtct                                                      >  1:3384476/1‑91 (MQ=255)
                             gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGttt                                                      <  1:1271262/93‑1 (MQ=255)
                              tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                                     >  2:287573/1‑93 (MQ=255)
                              tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGtttt                                                     <  1:183097/93‑1 (MQ=255)
                               aTTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTg                                                    <  2:292621/93‑1 (MQ=255)
                                tttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGGATGCTCCGGCCTTCGTTTTCATCACTGTGtttttt                                                   >  1:2324666/1‑91 (MQ=255)
                                    aCGGGGAGTTATGCCTTCTTTCCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTAc                                               <  1:2645652/93‑1 (MQ=255)
                                    aCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTAc                                               <  2:1672795/93‑1 (MQ=255)
                                      ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGCCGGAGCACGCTCCCGCCTTCGTTTTCATCACTGTGTTTTTAttttag                                             >  2:2290983/1‑89 (MQ=255)
                                      ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                             >  2:4174434/1‑93 (MQ=255)
                                      ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                             <  1:943687/93‑1 (MQ=255)
                                      ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAg                                             >  2:1646558/1‑93 (MQ=255)
                                       gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                            >  2:5113872/1‑93 (MQ=255)
                                       gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAgc                                            <  2:5130838/93‑1 (MQ=255)
                                           gTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcgg                                        <  2:2941902/93‑1 (MQ=255)
                                            ttATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                       >  1:5943837/1‑93 (MQ=255)
                                            ttATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                       <  2:5187577/93‑1 (MQ=255)
                                            ttATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                       >  2:2399373/1‑93 (MQ=255)
                                            ttATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGCGTTTTGATTACAGCgcggc                                       <  2:5288824/93‑1 (MQ=255)
                                            ttATGCCTCCTTTGCCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                       <  1:5574536/93‑1 (MQ=255)
                                                gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                                   <  2:1283198/93‑1 (MQ=255)
                                                gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                                   >  2:1438133/1‑93 (MQ=255)
                                                gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATa                                   <  1:2662394/93‑1 (MQ=255)
                                                gCCTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCCGCCTTCGTTTTTATTCCTGTGTTTTTATTTAAGAGCGGCGCTa                                   >  2:5397587/1‑93 (MQ=255)
                                                 cctcctTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                  <  1:5427323/93‑1 (MQ=255)
                                                 ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                  >  2:4066487/1‑93 (MQ=255)
                                                 ccTTCTTTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                  >  1:3481991/1‑93 (MQ=255)
                                                  cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                                 <  2:1881536/93‑1 (MQ=255)
                                                  cttcttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                                 >  1:5191356/1‑93 (MQ=255)
                                                    tcttTACTATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAg                               <  1:3109961/93‑1 (MQ=255)
                                                    tcttTACCATTTCTGTAAGAGATCAAACAGGCCGGAGCATGCTCCGGCCCTCGGTTTCACCACTGTGTTTTGCTTACAGCGCGGCGATCACAg                               >  1:4333447/1‑93 (MQ=255)
                                                     cttTACCATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGACTACAGCGCGGCGATAACAGc                              <  2:3588329/93‑1 (MQ=255)
                                                          ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                         <  1:3256146/93‑1 (MQ=255)
                                                          ccATTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgct                         <  2:4121392/93‑1 (MQ=255)
                                                            aTTTCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGt                       <  1:1695050/93‑1 (MQ=255)
                                                             tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTTCGGGCTTCGTTTTTATCACTGTGTTTTGATTTCCACGCGGGGATAACAGCCTTCTGtt                      >  2:214538/1‑93 (MQ=255)
                                                             tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt                      >  1:5283562/1‑93 (MQ=255)
                                                             tttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt                      >  1:2096007/1‑93 (MQ=255)
                                                              ttCTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                     <  2:2534234/93‑1 (MQ=255)
                                                                cTGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCaa                   >  1:249688/1‑93 (MQ=255)
                                                                 tGTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt                  >  2:5887846/1‑93 (MQ=255)
                                                                  gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc                 <  2:729798/93‑1 (MQ=255)
                                                                  gTAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATc                 >  1:3013827/1‑93 (MQ=255)
                                                                   tAAGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCa                <  2:2167961/93‑1 (MQ=255)
                                                                    aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAg               <  2:5520033/93‑1 (MQ=255)
                                                                    aaGAGATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAg               <  2:4671825/93‑1 (MQ=255)
                                                                     agagATAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGt              >  1:1024095/1‑93 (MQ=255)
                                                                      gagaTAAAAAAGGCCGGAGCATGCTTCGGCCTTCGTTTTTTTCACTGTGTTTTGATTTCAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt             >  2:282747/1‑93 (MQ=255)
                                                                      gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt             <  1:370065/93‑1 (MQ=255)
                                                                      gagaTAAAAAAGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt             <  2:4365531/93‑1 (MQ=255)
                                                                        agataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt            <  1:3342218/89‑1 (MQ=255)
                                                                        agaaaaaaaaGCCCGGAGCATGCCCGGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt            <  2:1099406/90‑1 (MQ=255)
                                                                         gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttc           >  2:5055688/4‑93 (MQ=255)
                                                                         gataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGAGCGGCGATAACAGCCTGCTGTTCAATCAGtttc           >  1:5619570/4‑93 (MQ=255)
                                                                          ataaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcg          >  1:3665966/3‑93 (MQ=255)
                                                                           taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc         <  1:4623977/92‑1 (MQ=255)
                                                                           taaaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc         <  2:6007282/92‑1 (MQ=255)
                                                                             aaaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctt       <  1:4947532/93‑1 (MQ=255)
                                                                              aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgcttt      <  1:6192969/93‑1 (MQ=255)
                                                                              aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgcttt      >  1:2142186/1‑93 (MQ=255)
                                                                              aaaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgcttt      <  2:3092568/93‑1 (MQ=255)
                                                                               aaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttc     >  1:3185260/1‑93 (MQ=255)
                                                                               aaaGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttc     >  1:2048322/1‑93 (MQ=255)
                                                                                 aGGCCGGAGCATGCTCCGGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttcgc   >  2:2486573/1‑93 (MQ=255)
                                                                                         |                                                                                       
AGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCT  >  NC_000913/4478422‑4478598

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