Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 478,533 C→T A29T (GCC→ACC)  ylaB ← putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,5330CT100.0% 297.9 / NA 90A29T (GCC→ACC) ylaBputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (49/41);  total (49/41)

CTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGCCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAA  >  NC_000913/478407‑478667
                                                                                                                              |                                                                                                                                      
cTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGACAGATGCt                                                                                                                            <  1:1127562/139‑1 (MQ=255)
  tGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTa                                                                                                                          >  1:516914/1‑139 (MQ=255)
   gCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTaa                                                                                                                         >  2:1332096/1‑139 (MQ=255)
               gCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                             <  1:1462492/139‑1 (MQ=255)
                cTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                                  <  1:195640/117‑1 (MQ=255)
                cTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                                  >  2:195640/1‑117 (MQ=255)
                     gCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATAc                                                                                                       >  1:669013/1‑139 (MQ=255)
                     gCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATAc                                                                                                       >  1:317493/1‑139 (MQ=255)
                     gCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATAc                                                                                                       >  2:138953/1‑139 (MQ=255)
                           cATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAg                                                                                                 >  2:640695/1‑139 (MQ=255)
                           cATCTTTCCATTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGACTCCAGCTCTTCAATAAACGATGTATCTACCTGCTGATGGGTCTGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAg                                                                                                 >  2:1540732/1‑139 (MQ=255)
                                 tCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAg                                                                                           >  1:1094841/1‑139 (MQ=255)
                                   ccTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                         >  1:902976/1‑139 (MQ=255)
                                   ccTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                         <  1:640695/139‑1 (MQ=255)
                                        tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                                   >  2:253828/1‑124 (MQ=255)
                                        tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                                   <  1:253828/124‑1 (MQ=255)
                                        tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTccc                                                                                    >  1:221100/1‑139 (MQ=255)
                                          tCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGa                                                                                  >  2:1378581/1‑139 (MQ=255)
                                          tCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGa                                                                                  >  1:1605095/1‑139 (MQ=255)
                                           cGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTg                                                                                                     <  1:1498039/119‑1 (MQ=255)
                                           cGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTg                                                                                                     >  2:1498039/1‑119 (MQ=255)
                                             ccACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAAt                                                                               >  2:1099296/1‑139 (MQ=255)
                                                  ttATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGc                                                                          <  1:428841/139‑1 (MQ=255)
                                                      tGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGAc                                                                      <  2:317493/139‑1 (MQ=255)
                                                        gCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCa                                                                    >  1:1011062/1‑139 (MQ=255)
                                                           cGAATAGCGACGCGGGCGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                         <  2:1092663/115‑1 (MQ=255)
                                                           cGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                         >  1:1092663/1‑115 (MQ=255)
                                                           cGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGAt                                                                 >  2:106073/1‑139 (MQ=255)
                                                            gAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCaa                                                                                                          <  1:1654877/97‑1 (MQ=255)
                                                            gAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCaa                                                                                                          >  2:1654877/1‑97 (MQ=255)
                                                            gAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATg                                                                >  1:1648323/1‑139 (MQ=255)
                                                             aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                             >  1:144192/1‑93 (MQ=255)
                                                             aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                             <  2:144192/93‑1 (MQ=255)
                                                             aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTcc                                                                                     <  1:1339306/117‑1 (MQ=255)
                                                             aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTcc                                                                                     >  2:1339306/1‑117 (MQ=255)
                                                             aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGATCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGt                                                               >  2:223680/1‑139 (MQ=255)
                                                                 gCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCgtg                                                           <  2:669013/139‑1 (MQ=255)
                                                                  cGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                         <  2:1200875/108‑1 (MQ=255)
                                                                  cGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                         >  1:1200875/1‑108 (MQ=255)
                                                                   gACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGtt                                                         >  1:319083/1‑139 (MQ=255)
                                                                   gACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGtt                                                         <  2:1011062/139‑1 (MQ=255)
                                                                   gACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGtt                                                         <  1:1378581/139‑1 (MQ=255)
                                                                     cgcgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTtct                                                       <  2:516914/139‑1 (MQ=255)
                                                                        gggAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCAc                                                    <  1:510199/139‑1 (MQ=255)
                                                                                 aGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                                           <  1:1332096/139‑1 (MQ=255)
                                                                                 aGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                                           <  1:138953/139‑1 (MQ=255)
                                                                                  ggTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAg                                                                   <  2:1043774/114‑1 (MQ=255)
                                                                                  ggTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAACACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAg                                                                   >  1:1043774/1‑114 (MQ=255)
                                                                                     aTCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCATGTTCTCACCAGCATATCCtta                                       >  2:581449/1‑139 (MQ=255)
                                                                                      tCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGAcc                                                                     <  2:1676729/108‑1 (MQ=255)
                                                                                      tCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGAcc                                                                     >  1:1676729/1‑108 (MQ=255)
                                                                                       ccAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCAAGATGCTTAAGCCGACAGgcagc                                                                                                            <  1:1523732/68‑1 (MQ=255)
                                                                                       ccAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCAAGATGCTTAAGCCGACAGgcagc                                                                                                            >  2:1523732/1‑68 (MQ=255)
                                                                                        cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTcc                                                                                     <  1:119186/90‑1 (MQ=255)
                                                                                        cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTcc                                                                                     >  2:119186/1‑90 (MQ=255)
                                                                                        cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCttatta                                    >  2:927842/1‑139 (MQ=255)
                                                                                        cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCttatta                                    >  2:207352/1‑139 (MQ=255)
                                                                                         aGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTGAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTaa                                   <  1:1099296/139‑1 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAgg                                                                                                                >  1:1093885/1‑60 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAgg                                                                                                                <  2:1093885/60‑1 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAg                                                                   >  1:593135/1‑105 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAg                                                                   <  2:593135/105‑1 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCata                                             >  1:521508/1‑127 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCata                                             <  2:521508/127‑1 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCa                                               >  1:427355/1‑125 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCa                                               <  2:427355/125‑1 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 <  2:902976/139‑1 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 >  2:579657/1‑139 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 <  2:319083/139‑1 (MQ=255)
                                                                                           ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 <  1:653879/139‑1 (MQ=255)
                                                                                            tctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAacac                                >  2:1610446/1‑139 (MQ=255)
                                                                                            tctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAacac                                <  2:1648323/139‑1 (MQ=255)
                                                                                             ctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTcc                                                                                     >  1:1539532/1‑85 (MQ=255)
                                                                                             ctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTcc                                                                                     <  2:1539532/85‑1 (MQ=255)
                                                                                                          tGTTTGTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGc                                                                          <  2:1099194/83‑1 (MQ=255)
                                                                                                          tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGc                                                                          >  1:1099194/1‑83 (MQ=255)
                                                                                                             ttCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 <  2:320572/121‑1 (MQ=255)
                                                                                                             ttCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 >  1:320572/1‑121 (MQ=255)
                                                                                                                  ccTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 >  1:1176519/1‑116 (MQ=255)
                                                                                                                  ccTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 <  2:1176519/116‑1 (MQ=255)
                                                                                                                     gctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATtg       <  1:207352/139‑1 (MQ=255)
                                                                                                                      ctgATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                             <  1:487527/36‑1 (MQ=255)
                                                                                                                      ctgATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                             >  2:487527/1‑36 (MQ=255)
                                                                                                                      ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAg                                                                   >  1:1051518/1‑78 (MQ=255)
                                                                                                                      ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAg                                                                   <  2:1051518/78‑1 (MQ=255)
                                                                                                                      ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 >  2:435944/1‑112 (MQ=255)
                                                                                                                      ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                                 <  1:435944/112‑1 (MQ=255)
                                                                                                                         aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGcc                                                                         >  1:730489/1‑69 (MQ=255)
                                                                                                                         aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGcc                                                                         <  2:730489/69‑1 (MQ=255)
                                                                                                                         aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGa   >  2:1147015/1‑139 (MQ=255)
                                                                                                                         aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGa   >  2:1450459/1‑139 (MQ=255)
                                                                                                                          tGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGaa  >  2:604689/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                      
CTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGCCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGTGGAA  >  NC_000913/478407‑478667

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