Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,994,415 G→A pseudogene (34/60 nt) pbl → pseudogene, peptidoglycan‑binding enzyme family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,994,4150GA100.0% 279.1 / NA 86pseudogene (34/60 nt)pblpseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (45/41);  total (45/41)

TAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATT  >  NC_000913/2994303‑2994540
                                                                                                                |                                                                                                                             
tAACGTGGCGATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACa                                                                                                                            <  2:1548601/116‑1 (MQ=255)
tAACGTGGCGATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACa                                                                                                                            >  1:1548601/1‑116 (MQ=255)
   cGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGcc                                                                                                  >  1:779334/1‑139 (MQ=255)
   cGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGcc                                                                                                  >  1:1430777/1‑139 (MQ=255)
     tGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTa                                                                                                >  2:1190958/1‑139 (MQ=255)
     tGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTa                                                                                                >  1:127988/1‑139 (MQ=255)
       gCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGt                                                                                                                       <  2:98055/114‑1 (MQ=255)
       gCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGt                                                                                                                       >  1:98055/1‑114 (MQ=255)
             ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        <  2:779334/139‑1 (MQ=255)
             ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        <  1:1557665/139‑1 (MQ=255)
                              cGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        <  1:972992/122‑1 (MQ=255)
                              cGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        >  2:972992/1‑122 (MQ=255)
                              cGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCa                                                                       >  1:456736/1‑139 (MQ=255)
                                ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCct                                                                                          >  1:1090248/1‑118 (MQ=255)
                                ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCct                                                                                          <  2:1090248/118‑1 (MQ=255)
                                  ttttGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        <  2:658666/117‑1 (MQ=255)
                                  ttttGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        >  1:658666/1‑117 (MQ=255)
                                     tGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACAccc                                                                <  1:1428651/139‑1 (MQ=255)
                                       cAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCAt                                                              <  2:279182/139‑1 (MQ=255)
                                                          gagaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                           <  2:127988/139‑1 (MQ=255)
                                                          gagaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                           <  2:1430777/139‑1 (MQ=255)
                                                           agaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                           >  2:397161/1‑138 (MQ=255)
                                                           agaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                           <  1:397161/138‑1 (MQ=255)
                                                                   atatGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                  <  1:1190958/139‑1 (MQ=255)
                                                                    tatGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATAc                                 >  1:1012822/1‑139 (MQ=255)
                                                                      tGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTa                                                         <  2:1350295/113‑1 (MQ=255)
                                                                      tGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTa                                                         >  1:1350295/1‑113 (MQ=255)
                                                                           cAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        <  1:331987/77‑1 (MQ=255)
                                                                           cAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        >  2:331987/1‑77 (MQ=255)
                                                                           cAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACtgtg                          >  1:1287210/1‑139 (MQ=255)
                                                                                       tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt                                                                                   <  2:1117550/70‑1 (MQ=255)
                                                                                       tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttt                                                                                   >  1:1117550/1‑70 (MQ=255)
                                                                                        gTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                             >  1:749131/1‑75 (MQ=255)
                                                                                        gTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                             <  2:749131/75‑1 (MQ=255)
                                                                                         tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                     <  1:1322383/82‑1 (MQ=255)
                                                                                         tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                     >  2:1322383/1‑82 (MQ=255)
                                                                                         tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATttct            <  2:1012822/139‑1 (MQ=255)
                                                                                          atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        >  2:354460/1‑62 (MQ=255)
                                                                                          atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        <  1:354460/62‑1 (MQ=255)
                                                                                          atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaa                                                                      <  1:800459/80‑1 (MQ=255)
                                                                                          atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaa                                                                      >  2:800459/1‑80 (MQ=255)
                                                                                          atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATttctt           >  2:1519045/1‑139 (MQ=255)
                                                                                           tataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCa         >  1:506935/1‑139 (MQ=255)
                                                                                            ataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatata                                                                            >  2:1354698/1‑72 (MQ=255)
                                                                                            ataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatata                                                                            <  1:1354698/72‑1 (MQ=255)
                                                                                            ataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCa         <  1:1519045/139‑1 (MQ=255)
                                                                                              aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAAt                                                                                                       <  2:894239/43‑1 (MQ=255)
                                                                                              aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAAt                                                                                                       >  1:894239/1‑43 (MQ=255)
                                                                                              aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                     >  1:946489/1‑77 (MQ=255)
                                                                                              aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                     <  2:946489/77‑1 (MQ=255)
                                                                                              aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                       <  2:790423/107‑1 (MQ=255)
                                                                                              aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                       >  1:790423/1‑107 (MQ=255)
                                                                                              aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCa         >  2:1593239/1‑137 (MQ=255)
                                                                                              aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCa         <  1:1593239/137‑1 (MQ=255)
                                                                                              aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATg       >  2:1441833/1‑139 (MQ=255)
                                                                                               ccAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                             <  1:1166847/68‑1 (MQ=255)
                                                                                               ccAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTtatat                                                                             >  2:1166847/1‑68 (MQ=255)
                                                                                                cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        <  1:1552973/56‑1 (MQ=255)
                                                                                                cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        >  2:1552973/1‑56 (MQ=255)
                                                                                                cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGat     >  1:1493741/1‑139 (MQ=255)
                                                                                                cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGAAAACACCCATACGTAAAGCTGTTTTGAt                                            >  2:464839/1‑100 (MQ=255)
                                                                                                cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGAAAACACCCATACGTAAAGCTGTTTTGAt                                            <  1:464839/100‑1 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        >  2:375538/1‑53 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                        <  1:375538/53‑1 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttttt                                                                                 >  2:1213844/1‑60 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttttt                                                                                 <  1:1213844/60‑1 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                     <  2:563656/72‑1 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                     >  1:563656/1‑72 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaa                                                                      <  1:87775/71‑1 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaa                                                                      >  2:87775/1‑71 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGtt                                                 >  2:1422972/1‑92 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                  <  1:257386/107‑1 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                  >  2:257386/1‑107 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCAt        >  2:206594/1‑133 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCAt        <  1:206594/133‑1 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt  <  2:456736/139‑1 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt  >  1:851575/1‑139 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt  >  2:1310783/1‑139 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTCAAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt  >  1:751023/1‑139 (MQ=255)
                                                                                                   gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCGCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGtt                                                 <  1:1422972/92‑1 (MQ=255)
                                                                                                       aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATttct            >  2:690453/1‑125 (MQ=255)
                                                                                                       aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATttct            <  1:690453/125‑1 (MQ=255)
                                                                                                        aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACtgtg                          <  2:404060/110‑1 (MQ=255)
                                                                                                        aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACtgtg                          >  1:404060/1‑110 (MQ=255)
                                                                                                             cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaa                                                                      <  2:1287447/61‑1 (MQ=255)
                                                                                                             cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaa                                                                      >  1:1287447/1‑61 (MQ=255)
                                                                                                                |                                                                                                                             
TAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATT  >  NC_000913/2994303‑2994540

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: