Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,572,526 T→C D234G (GAC→GGC)  pqqL ← putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,572,5260TC100.0% 295.6 / NA 93D234G (GAC→GGC) pqqLputative peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (52/41);  total (52/41)

GCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGTCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGC  >  NC_000913/1572436‑1572614
                                                                                          |                                                                                        
gCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCg                                                                                        >  2:1292852/1‑93 (MQ=255)
gCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCg                                                                                        >  2:3856200/1‑93 (MQ=255)
 cTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGa                                                                                       <  2:481732/93‑1 (MQ=255)
 cTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGa                                                                                       <  1:2893430/93‑1 (MQ=255)
  ttttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGat                                                                                      >  1:3802222/1‑93 (MQ=255)
  ttttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGat                                                                                      >  1:3414873/1‑93 (MQ=255)
  ttttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGat                                                                                      >  2:206438/1‑93 (MQ=255)
   tttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGata                                                                                     >  2:3210485/1‑93 (MQ=255)
    ttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGatat                                                                                    <  2:914604/93‑1 (MQ=255)
    ttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGatat                                                                                    <  2:341567/93‑1 (MQ=255)
    ttGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGatat                                                                                    <  1:2414505/93‑1 (MQ=255)
    ttGCCGGCCAGACGCTGTTTTCAGCTGCTTTGTAAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGatat                                                                                    <  1:3627338/93‑1 (MQ=255)
      gTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCg                                                                                  >  1:293881/1‑93 (MQ=255)
         ggCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCg                                                                               >  1:3182434/1‑93 (MQ=255)
         ggCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCg                                                                               >  2:448116/1‑93 (MQ=255)
         ggCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCg                                                                               >  1:4060753/1‑93 (MQ=255)
             aGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCa                                                                           <  1:671047/93‑1 (MQ=255)
               aCGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACg                                                                         >  1:3994791/1‑93 (MQ=255)
                cgcgATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGa                                                                        <  2:1140400/93‑1 (MQ=255)
                  cgATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATa                                                                      <  2:230588/93‑1 (MQ=255)
                  cgATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATa                                                                      >  1:1991221/1‑93 (MQ=255)
                    aTTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATaaa                                                                    >  2:1888305/1‑93 (MQ=255)
                      tttCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAgg                                                                  >  2:902653/1‑93 (MQ=255)
                      tttCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAgg                                                                  >  2:1738195/1‑93 (MQ=255)
                        tCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTc                                                                <  2:3071502/93‑1 (MQ=255)
                        tCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTc                                                                >  1:3378897/1‑93 (MQ=255)
                        tCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTc                                                                <  2:1527544/93‑1 (MQ=255)
                          aGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCat                                                              >  1:364795/1‑93 (MQ=255)
                              gctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGAGCACGATAAAGGTCAtatta                                                          >  2:3968108/1‑93 (MQ=255)
                               ctTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCAtattat                                                         >  1:825202/1‑93 (MQ=255)
                                tttGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATattatt                                                        <  1:4020711/93‑1 (MQ=255)
                                 ttGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTAttt                                                       <  2:3567682/93‑1 (MQ=255)
                                 ttGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTAttt                                                       <  2:3658574/93‑1 (MQ=255)
                                  tGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTg                                                      <  1:911339/93‑1 (MQ=255)
                                   gTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTgg                                                     <  1:175808/93‑1 (MQ=255)
                                   gTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTgg                                                     <  1:2310090/93‑1 (MQ=255)
                                    ttAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGt                                                    <  2:1133627/93‑1 (MQ=255)
                                    ttAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGt                                                    >  1:1388944/1‑93 (MQ=255)
                                       gCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGa                                                 <  2:125293/93‑1 (MQ=255)
                                         cGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATa                                               <  1:1524518/93‑1 (MQ=255)
                                          ggAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATAc                                              >  2:4661632/1‑93 (MQ=255)
                                           gAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATAcc                                             >  1:1280329/1‑93 (MQ=255)
                                            aaGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCa                                            <  1:985842/93‑1 (MQ=255)
                                            aaGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCa                                            <  1:1457783/93‑1 (MQ=255)
                                              gCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGc                                          >  1:2196356/1‑93 (MQ=255)
                                              gCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGc                                          <  2:3301814/93‑1 (MQ=255)
                                              gCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGc                                          <  2:3974802/93‑1 (MQ=255)
                                              gCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGc                                          >  1:3471276/1‑93 (MQ=255)
                                               cttacttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCg                                         <  1:4234239/93‑1 (MQ=255)
                                               cttacttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCg                                         >  2:2087292/1‑93 (MQ=255)
                                                  acttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTg                                      >  1:2693791/1‑93 (MQ=255)
                                                  acttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTg                                      <  1:1889937/93‑1 (MQ=255)
                                                  acttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTg                                      >  2:2291868/1‑93 (MQ=255)
                                                   cttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGa                                     >  2:976874/1‑93 (MQ=255)
                                                   cttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGa                                     <  2:1777601/93‑1 (MQ=255)
                                                    ttaAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGAt                                    <  2:630688/93‑1 (MQ=255)
                                                     taAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATa                                   <  1:2728797/93‑1 (MQ=255)
                                                     taAATTATCCTTTATCAGCGCCAGCGATTCTTTACTGCCGCTATCGCCGACCCCGCTAAAGGTCATATTATTTGGTTGCTACCAGCGTTGATa                                   >  1:1684789/1‑93 (MQ=255)
                                                      aaaTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaa                                  >  1:2439056/1‑93 (MQ=255)
                                                      aaaTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaa                                  >  1:4226463/1‑93 (MQ=255)
                                                      aaaTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaa                                  <  2:1121178/93‑1 (MQ=255)
                                                       aaTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaaa                                 >  1:1455600/1‑93 (MQ=255)
                                                        aTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaaaa                                <  1:353510/93‑1 (MQ=255)
                                                        aTTATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATaaaa                                <  1:3123116/93‑1 (MQ=255)
                                                          tATCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAAtt                              >  1:3591175/1‑93 (MQ=255)
                                                           aTCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTg                             >  2:2960011/1‑93 (MQ=255)
                                                            tCCTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTgg                            <  1:4320699/93‑1 (MQ=255)
                                                              cTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgcg                          >  1:3624073/1‑93 (MQ=255)
                                                              cTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgcg                          >  1:940221/1‑93 (MQ=255)
                                                              cTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgcg                          >  1:1018004/1‑93 (MQ=255)
                                                              cTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgcg                          >  1:1996028/1‑93 (MQ=255)
                                                              cTTTATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgcg                          >  2:1703783/1‑93 (MQ=255)
                                                               tttATCAGCGCCAGCGCTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGgcgc                         >  1:3883104/1‑93 (MQ=255)
                                                                     agcgccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATtgt                   <  2:15741/93‑1 (MQ=255)
                                                                      gcgccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATtgtg                  <  2:1958726/93‑1 (MQ=255)
                                                                      gcgccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATtgtg                  >  2:2728155/1‑93 (MQ=255)
                                                                       cgccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGc                 <  2:180367/93‑1 (MQ=255)
                                                                        gccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGcc                >  1:839421/1‑93 (MQ=255)
                                                                        gccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGcc                >  2:1409601/1‑93 (MQ=255)
                                                                         ccagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCg               >  1:3087835/1‑93 (MQ=255)
                                                                          cagcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCgg              >  2:212509/1‑93 (MQ=255)
                                                                           agcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgt             >  2:526436/1‑93 (MQ=255)
                                                                            gcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgtg            >  1:2560875/1‑93 (MQ=255)
                                                                            gcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgtg            >  2:1695451/1‑93 (MQ=255)
                                                                            gcgcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGgtg            >  1:404096/1‑93 (MQ=255)
                                                                              gcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGCTAAAATTGGCGCAATTGTGCCGGtgtg          >  1:3659347/1‑93 (MQ=255)
                                                                              gcTTCTTTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGtgtg          <  1:2122046/93‑1 (MQ=255)
                                                                               cttcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGa         <  2:2368708/93‑1 (MQ=255)
                                                                                ttcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGAc        <  2:3519887/93‑1 (MQ=255)
                                                                                 tcttTACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGAcc       >  2:3606073/1‑93 (MQ=255)
                                                                                   tttACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGt     >  2:2575999/1‑93 (MQ=255)
                                                                                    ttgcTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATTCCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGGCCGTg    <  2:2084366/90‑1 (MQ=255)
                                                                                     tACTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTgg   >  2:2675412/1‑93 (MQ=255)
                                                                                      tcTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGc  <  1:4620595/92‑1 (MQ=255)
                                                                                      aCTGCCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGc  <  1:3604845/93‑1 (MQ=255)
                                                                                          |                                                                                        
GCTTTTGTCGGCCAGACGCGATTTTCAGCTGCTTTGTTAGCCGGAAGCTTACTTAAATTATCCTTTATCAGCGCCAGCGCTTCTTTACTGTCGATATCGCCGACCACGATAAAGGTCATATTATTTGGTTGATACCAGCGTTGATAAAATTGGCGCAATTGTGCCGGTGTGACCGTGGC  >  NC_000913/1572436‑1572614

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