Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,426,288 G→A intergenic (+206/‑166) ydaY → / → tmpR pseudogene, Rac prophage,Phage or Prophage Related/Rac prophage, pseudogene, tail protein family,Phage or Prophage Related, putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,2880GA100.0% 249.7 / NA 77intergenic (+206/‑166)ydaY/tmpRpseudogene, Rac prophage,Phage or Prophage Related/Rac prophage, pseudogene, tail protein family,Phage or Prophage Related, putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/42);  total (35/42)

GGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGG  >  NC_000913/1426154‑1426423
                                                                                                                                      |                                                                                                                                       
ggggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAattgt                                                                                                                                     <  1:958548/139‑1 (MQ=255)
  ggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTa                                                                                                                                   >  2:1458649/1‑139 (MQ=255)
    tgcctttgccttTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAAt                                                                                                                                 <  1:593641/139‑1 (MQ=255)
            cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTaca                                                                                                                         >  1:663332/1‑139 (MQ=255)
              tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       <  1:678113/139‑1 (MQ=255)
              tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       <  2:935668/139‑1 (MQ=255)
              tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       <  2:1530021/139‑1 (MQ=255)
              tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       <  2:890542/139‑1 (MQ=255)
              tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       <  1:874225/139‑1 (MQ=255)
              tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       <  2:761462/139‑1 (MQ=255)
                    tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTa                                                                                                                 <  2:663332/139‑1 (MQ=255)
                        aCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACaga                                                                                                             <  1:339461/139‑1 (MQ=255)
                            agaTAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGAtg                                                                                                         >  2:1596679/1‑139 (MQ=255)
                                  aTGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       <  2:692010/119‑1 (MQ=255)
                                  aTGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       >  1:692010/1‑119 (MQ=255)
                                      atatTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTc                                                                                               <  1:1458649/139‑1 (MQ=255)
                                        atTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCtt                                                                                             >  2:1170503/1‑139 (MQ=255)
                                            aGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTa                                                                                                                           <  2:1300126/105‑1 (MQ=255)
                                            aGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTa                                                                                                                           >  1:1300126/1‑105 (MQ=255)
                                                 cTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAAt                                                                                    >  1:644323/1‑139 (MQ=255)
                                                    aTCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCga                                                                                 >  1:522596/1‑139 (MQ=255)
                                                       gggTATTGATTGAGGGGAAGGAATTTAAAGAAGTGTTTTTTGCATTGGATAACAAGAATTTCAGGAGCAGTGTGGTTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgatgat                                                                             <  1:686690/138‑1 (MQ=255)
                                                         gTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGatgatg                                                                            >  2:1402228/1‑139 (MQ=255)
                                                            ttgattgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTc                                                                         <  1:1596679/139‑1 (MQ=255)
                                                              gattgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCt                                                                       >  1:779581/1‑139 (MQ=255)
                                                               attgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTaca                                                                                                                         <  2:1176845/88‑1 (MQ=255)
                                                               attgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTaca                                                                                                                         >  1:1176845/1‑88 (MQ=255)
                                                                   aGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAAt                                                                                    >  2:679436/1‑121 (MQ=255)
                                                                   aGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAAt                                                                                    <  1:679436/121‑1 (MQ=255)
                                                                    gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTa                                                                                                                                   <  2:1057848/73‑1 (MQ=255)
                                                                    gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTa                                                                                                                                   >  1:1057848/1‑73 (MQ=255)
                                                                    gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       >  1:559419/1‑85 (MQ=255)
                                                                    gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACAct                                                                                                                       <  2:559419/85‑1 (MQ=255)
                                                                    gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACag                                                                                                              >  1:645542/1‑94 (MQ=255)
                                                                    gcgcGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACag                                                                                                              <  2:645542/94‑1 (MQ=255)
                                                                         aTGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACa                                                                                                               >  1:144578/1‑88 (MQ=255)
                                                                         aTGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACa                                                                                                               <  2:144578/88‑1 (MQ=255)
                                                                                      aGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCg                                               >  1:29003/1‑139 (MQ=255)
                                                                                        tGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTa                                                                                                                                   >  1:1380512/1‑53 (MQ=255)
                                                                                        tGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTa                                                                                                                                   <  2:1380512/53‑1 (MQ=255)
                                                                                        tGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgatga                                                                              >  1:1284188/1‑106 (MQ=255)
                                                                                        tGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgatga                                                                              <  2:1284188/106‑1 (MQ=255)
                                                                                        tGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGatgatg                                                                            >  1:461709/1‑108 (MQ=255)
                                                                                        tGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGatgatg                                                                            <  2:461709/108‑1 (MQ=255)
                                                                                           tCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATc                                          <  2:29003/139‑1 (MQ=255)
                                                                                           tCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATc                                          <  1:172839/139‑1 (MQ=255)
                                                                                           tCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATc                                          <  2:644323/139‑1 (MQ=255)
                                                                                                  aCTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATg                                                                                                    <  2:888184/74‑1 (MQ=255)
                                                                                                  aCTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATg                                                                                                    >  1:888184/1‑74 (MQ=255)
                                                                                                      gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTa                                                                                                                 >  1:23874/1‑57 (MQ=255)
                                                                                                      gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTa                                                                                                                 <  2:23874/57‑1 (MQ=255)
                                                                                                      gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATg                                                                                                    <  2:676929/70‑1 (MQ=255)
                                                                                                      gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATg                                                                                                    >  1:676929/1‑70 (MQ=255)
                                                                                                      gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                       >  1:153473/1‑131 (MQ=255)
                                                                                                      gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                       <  2:153473/131‑1 (MQ=255)
                                                                                                      gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATc                                  >  2:549028/1‑136 (MQ=255)
                                                                                                      gATGACCAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATc                                  <  1:549028/136‑1 (MQ=255)
                                                                                                           ccAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCa                                                            >  1:1131200/1‑105 (MQ=255)
                                                                                                           ccAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCa                                                            <  2:1131200/105‑1 (MQ=255)
                                                                                                            cAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATgga                         >  1:1170503‑M2/1‑136 (MQ=255)
                                                                                                            cAGTATTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATgga                         >  1:1402228‑M2/1‑136 (MQ=255)
                                                                                                                aTTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATgg                                                                                                   >  1:877014/1‑61 (MQ=255)
                                                                                                                aTTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATgg                                                                                                   <  2:877014/61‑1 (MQ=255)
                                                                                                                aTTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGc                                        >  1:916668/1‑120 (MQ=255)
                                                                                                                aTTTTCAGGAGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGc                                        <  2:916668/120‑1 (MQ=255)
                                                                                                                         aGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgatga                                                                              >  2:973547/1‑73 (MQ=255)
                                                                                                                         aGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgatga                                                                              >  2:1289857/1‑73 (MQ=255)
                                                                                                                         aGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgatga                                                                              <  1:973547/73‑1 (MQ=255)
                                                                                                                         aGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCgatga                                                                              <  1:1289857/73‑1 (MQ=255)
                                                                                                                         aGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGgaga                                                                  >  2:1366420/1‑85 (MQ=255)
                                                                                                                         aGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGgaga                                                                  <  1:1366420/85‑1 (MQ=255)
                                                                                                                         aGCAGTGTGTGTAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATggaaggtgcgaataag            <  2:1086822‑M2/139‑17 (MQ=255)
                                                                                                                             gtgtgtgtAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                       >  1:649060/1‑108 (MQ=255)
                                                                                                                             gtgtgtgtAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                       <  2:649060/108‑1 (MQ=255)
                                                                                                                              tgtgtgtAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGaa                                                                 <  1:1571726/81‑1 (MQ=255)
                                                                                                                              tgtgtgtAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGaa                                                                 >  2:1571726/1‑81 (MQ=255)
                                                                                                                              tgtgtgtAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATggaaggtgcgaataagcgggg       <  1:166467‑M2/139‑22 (MQ=255)
                                                                                                                                   gtAATTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATggaaggtgcgaataagcggggaaatt  <  2:82301‑M2/139‑27 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGG  >  NC_000913/1426154‑1426423

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