Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,572,647 A→G S852S (AGT→AGC pqqL ← putative periplasmic M16 family zinc metalloendopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,572,6470AG100.0% 392.3 / NA 107S852S (AGT→AGCpqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (51/56);  total (51/56)

Only 100 of 107 total aligned reads displayed.

TGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGA  >  NC_000913/1572513‑1572774
                                                                                                                                      |                                                                                                                                
tGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              <  1:382754/139‑1 (MQ=255)
       tCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATcc                                                                                                                           >  1:1342573/1‑135 (MQ=255)
                 aGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGc                                                                                                             >  2:152834/1‑139 (MQ=255)
                   gATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAACCAGACGCtt                                                                                                           <  1:378535/139‑1 (MQ=255)
                    aTCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTa                                                                                                          <  1:707994/139‑1 (MQ=255)
                             ttGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTa                                                                                                          <  2:859483/130‑1 (MQ=255)
                             ttGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTa                                                                                                          >  1:859483/1‑130 (MQ=255)
                                  tAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCAtt                                                                                            >  1:751978/1‑139 (MQ=255)
                                    aGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCg                                                                                          <  2:751978/139‑1 (MQ=255)
                                          tttACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              <  1:1295452/97‑1 (MQ=255)
                                          tttACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              >  2:1295452/1‑97 (MQ=255)
                                          tttACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCt                                                                                        >  1:1258112/1‑135 (MQ=255)
                                          tttACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCt                                                                                        <  2:1258112/135‑1 (MQ=255)
                                            tACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAAcc                                                                                                       <  2:32123/118‑1 (MQ=255)
                                            tACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAAcc                                                                                                       >  1:32123/1‑118 (MQ=255)
                                              cAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTaa                                                                                <  2:1269456/139‑1 (MQ=255)
                                               aaTAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACt                                                                                                 >  2:179186/1‑121 (MQ=255)
                                               aaTAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACt                                                                                                 <  1:179186/121‑1 (MQ=255)
                                                 tAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGc                                                                                                             <  2:1614692/107‑1 (MQ=255)
                                                    tgtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              <  1:419627/87‑1 (MQ=255)
                                                    tgtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              >  2:419627/1‑87 (MQ=255)
                                                    tgtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCAt                                                                                                     <  1:182511/112‑1 (MQ=255)
                                                    tgtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCAt                                                                                                     >  2:182511/1‑112 (MQ=255)
                                                    tgtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAAc                                                                                     >  2:851210/1‑128 (MQ=255)
                                                    tgtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAAc                                                                                     <  1:851210/128‑1 (MQ=255)
                                                    tgtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTaa                                                                                >  2:608690/1‑133 (MQ=255)
                                                    tgtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTaa                                                                                <  1:608690/133‑1 (MQ=255)
                                                    tgtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtc                                                                          >  1:608618/1‑139 (MQ=255)
                                                     gtTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                                                                         >  1:884393/1‑139 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGa                                                                                                                >  1:622992/1‑99 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGAc                                                                                                               <  1:408882/100‑1 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGAc                                                                                                               <  2:706101/100‑1 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGAc                                                                                                               >  2:408882/1‑100 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGAc                                                                                                               >  1:706101/1‑100 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTa                                                                                 >  2:1337722/1‑130 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTa                                                                                 <  1:1337722/130‑1 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACa                                                                              >  2:1556332/1‑133 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACa                                                                              >  2:600035/1‑133 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACa                                                                              <  1:1556332/133‑1 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACa                                                                              <  1:600035/133‑1 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcat                                                                        >  1:667726/1‑139 (MQ=255)
                                                      ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcat                                                                        >  1:953681/1‑139 (MQ=255)
                                                        cGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGa                                                                                                                <  2:314341/97‑1 (MQ=255)
                                                        cGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGa                                                                                                                >  1:314341/1‑97 (MQ=255)
                                                         gCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTcatcat                                                                     <  2:1000970/139‑1 (MQ=255)
                                                           tAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGTTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTaa                                                                                >  2:1378971/1‑126 (MQ=255)
                                                           tAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGTTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTaa                                                                                <  1:1378971/126‑1 (MQ=255)
                                                              ttgttgAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcatca                                                                      <  2:263773/133‑1 (MQ=255)
                                                              ttgttgAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcatca                                                                      >  1:263773/1‑133 (MQ=255)
                                                              ttgttgAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGt                                                                >  2:705740/1‑139 (MQ=255)
                                                                gttgAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              <  2:790467/75‑1 (MQ=255)
                                                                   gAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              >  2:115098/1‑72 (MQ=255)
                                                                   gAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              <  1:115098/72‑1 (MQ=255)
                                                                   gAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCg                                                                                          >  1:499232/1‑108 (MQ=255)
                                                                   gAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCg                                                                                          <  2:499232/108‑1 (MQ=255)
                                                                       gctacgctGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGa                                                       <  2:504452/139‑1 (MQ=255)
                                                                        ctacgctGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGAc                                                      <  2:884393/139‑1 (MQ=255)
                                                                        ctacgctGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGAc                                                      <  2:667726/139‑1 (MQ=255)
                                                                                ttGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCa                                                                                                                  >  2:853014/1‑71 (MQ=255)
                                                                                ttGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCa                                                                                                                  <  1:853014/71‑1 (MQ=255)
                                                                                  ggATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGc                                           <  1:152834/139‑1 (MQ=255)
                                                                                  ggATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGc                                           <  2:608618/139‑1 (MQ=255)
                                                                                        cGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCtt                                                                                                                        >  1:269484/1‑57 (MQ=255)
                                                                                        cGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCtt                                                                                                                        <  2:269484/57‑1 (MQ=255)
                                                                                        cGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTa                                                                                 <  2:672239/96‑1 (MQ=255)
                                                                                        cGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTa                                                                                 >  1:672239/1‑96 (MQ=255)
                                                                                        cGAGGCTGCCCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              >  1:1387492/1‑51 (MQ=38)
                                                                                        cGAGGCTGCCCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGa                                                                                                                              <  2:1387492/51‑1 (MQ=38)
                                                                                          aGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAAcc                                                                                                       >  2:742091/1‑72 (MQ=255)
                                                                                          aGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAAcc                                                                                                       <  1:742091/72‑1 (MQ=255)
                                                                                          aGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCg                                                                                          >  2:23495/1‑85 (MQ=255)
                                                                                          aGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCg                                                                                          <  1:23495/85‑1 (MQ=255)
                                                                                          aGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCt                                                                                        >  2:162169/1‑87 (MQ=255)
                                                                                          aGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCt                                                                                        <  1:162169/87‑1 (MQ=255)
                                                                                             cTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCt                                                                                        <  1:1530786/84‑1 (MQ=255)
                                                                                             cTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCt                                                                                        >  2:1530786/1‑84 (MQ=255)
                                                                                              tGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCGCTTCATTCGCTCACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACt                               >  1:957237/1‑139 (MQ=255)
                                                                                                cgcTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCa                                                                                                                  <  1:1495858/55‑1 (MQ=255)
                                                                                                cgcTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCa                                                                                                                  >  2:1495858/1‑55 (MQ=255)
                                                                                                cgcTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGa                             <  2:1561602/139‑1 (MQ=255)
                                                                                                cgcTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGa                             <  1:609354/139‑1 (MQ=255)
                                                                                                cgcTGAACGTTTTGCTGGTATTCATACAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGa                             >  1:1212416/1‑139 (MQ=255)
                                                                                                               tGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACa                                                                              >  1:1185216/1‑76 (MQ=255)
                                                                                                               tGGTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACa                                                                              <  2:1185216/76‑1 (MQ=255)
                                                                                                                 gTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCGTCACTTCATTCg                                                                                          <  1:1368172/62‑1 (MQ=255)
                                                                                                                 gTATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCg                                                                                          >  2:1368172/1‑62 (MQ=255)
                                                                                                                  tATTCATTCAGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGAt           <  2:953681/139‑1 (MQ=255)
                                                                                                                    ttcatacaGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCa         <  2:1212416/139‑1 (MQ=255)
                                                                                                                         tcaGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCga    <  1:705740/139‑1 (MQ=255)
                                                                                                                           aGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAAGCCAGCAAATGACTGa                             <  1:796755/113‑1 (MQ=255)
                                                                                                                           aGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAAGCCAGCAAATGACTGa                             >  2:796755/1‑113 (MQ=255)
                                                                                                                           aGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCaaa                                    >  2:35975/1‑105 (MQ=255)
                                                                                                                           aGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCaaa                                    <  1:35975/105‑1 (MQ=255)
                                                                                                                           aGTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCgaga  >  2:447620/1‑139 (MQ=255)
                                                                                                                            gTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGAc                                                      <  1:1438081/87‑1 (MQ=255)
                                                                                                                            gTTCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGAc                                                      >  2:1438081/1‑87 (MQ=255)
                                                                                                                             ttCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTcact                                                                                        >  1:749994/1‑49 (MQ=255)
                                                                                                                             ttCTTGCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTcact                                                                                        <  2:749994/52‑4 (MQ=255)
                                                                                                                                  gCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGc                                           >  1:511637/1‑91 (MQ=255)
                                                                                                                                  gCTCGCTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGc                                           <  2:511637/91‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
TGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGA  >  NC_000913/1572513‑1572774

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