Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,135,866 A→G 5.7% L78P (CTG→CCG)  wzb ← colanic acid production protein‑tyrosine‑phosphatase; Wzc‑P dephosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,135,8660AG 5.7% 194.3 / 2.4 72L78P (CTG→CCG) wzbcolanic acid production protein‑tyrosine‑phosphatase; Wzc‑P dephosphorylase
Reads supporting (aligned to +/- strand):  ref base A (35/30);  new base G (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAGTCCAGCGGACTCCACTTTC  >  NC_000913/2135718‑2136009
                                                                                                                                                    |                                                                                                                                               
cacCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCa                                                                                                                                                 >  1:140325/1‑149 (MQ=255)
  ccGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAgg                                                                                                                                               >  1:1099371/1‑149 (MQ=255)
       gCAAACGTTTCCCGGCTTTTGCGATACGGAGCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGGGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTa                                                                                                                                          <  2:1016955/149‑1 (MQ=255)
        cAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAg                                                                                                                                         >  1:220619/1‑149 (MQ=255)
        cAAACGTTTCACGGCTTTTGCGATACGCATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGATTATCCATGGTCAAAATCAGGTCGTAg                                                                                                                                         >  2:1170262/1‑149 (MQ=255)
          aaCGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGtt                                                                                                                                       <  2:119289/149‑1 (MQ=255)
            cGTTTCCCGGCTGTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGc                                                                                                                                     <  2:439815/149‑1 (MQ=255)
             gTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTg                                                                                                                                      <  1:1160064/147‑1 (MQ=255)
             gTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTg                                                                                                                                      >  2:1160064/1‑147 (MQ=255)
                  ccGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCGGGTCGTAGTTGCGAcaca                                                                                                                               >  2:803844/1‑149 (MQ=255)
                     gCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGAc                                                                                                                            <  2:110723/149‑1 (MQ=255)
                       ttttGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGAcgg                                                                                                                          >  2:533603/1‑149 (MQ=255)
                       ttttGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCATTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGAcgg                                                                                                                          <  2:579865/149‑1 (MQ=255)
                          tGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAACCAGCATCACATTGCCGCCCATCTCCGGTTCCATCTCGCATAAGCGTTCGATATGGCGCTTTTACATGGTCAAAATCAGGTCGTAGTTGCGACACAGAcggcgg                                                                                                                       <  2:285474/149‑1 (MQ=255)
                                    aTCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGGATAAGCGTTCGATATGGCGCTTTTCCATGGGCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGAc                                                                                                             <  1:610541/149‑1 (MQ=255)
                                    aTCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACGTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCGGATTTGAc                                                                                                             >  2:678186/1‑149 (MQ=255)
                                       ggggATTTCACATTCGTTATCCCAGTGACCAAACAGCCTCACTTTGCCGCGCATGTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACggg                                                                                                          <  1:397043/149‑1 (MQ=255)
                                       ggggATTTCACATTCGTAATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCCGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACggg                                                                                                          <  2:971535/149‑1 (MQ=255)
                                        gggATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGACGAGCCTCTACGTTGCCATCTCGCATAAGCGTTCGATATGGAGCTGTTACATGGTCAAAATCAGGTCGTAGTTGCGACACGCACGGCGGGTGCTTTGGCgggg                                                                                                         >  2:72584/1‑148 (MQ=255)
                                          gATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGcac                                                                                                       >  1:847681/1‑149 (MQ=255)
                                           aTTTCACATTAGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCGGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGcaca                                                                                                      <  1:246918/149‑1 (MQ=255)
                                              tCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCTGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCTAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTg                                                                                                   >  1:1104067/1‑149 (MQ=255)
                                              tCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTg                                                                                                   <  1:607410/149‑1 (MQ=255)
                                              tCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTCATTTGACGGGCACAGTg                                                                                                   >  1:868861/1‑149 (MQ=255)
                                               cacaTTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTg                                                                                                   <  1:382410/148‑1 (MQ=255)
                                               cacaTTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTg                                                                                                   >  2:382410/1‑148 (MQ=255)
                                                acaTTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGAc                                                                                                 <  1:905513/149‑1 (MQ=255)
                                                    tCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCCCAGTGACGTTc                                                                                             >  1:29585/1‑149 (MQ=255)
                                                          tCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAgaga                                                                                       >  1:493461/1‑149 (MQ=255)
                                                           cccAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGCCCTTCCAGAGAc                                                                                      >  1:546071/1‑149 (MQ=255)
                                                            ccAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTCGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACa                                                                                     <  2:313656/149‑1 (MQ=255)
                                                             cAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGGTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCGTCCAGCGACAg                                                                                    >  2:1160956/1‑149 (MQ=255)
                                                               gTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGtt                                                                                  <  1:865673/149‑1 (MQ=255)
                                                               gTGACCAAACAGCAGCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATCTGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCTGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGtt                                                                                  <  1:1075763/149‑1 (MQ=255)
                                                                tgcccAAAGAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGAATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCCGGTCGTAGTTGCGAGACAGACGGCGGCTTATTTGACGGGCACAGTGACCTTCCAGAGACAGTTg                                                                                 <  1:917052/146‑1 (MQ=255)
                                                                tGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTg                                                                                 <  1:332868/149‑1 (MQ=255)
                                                                   ccAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATg                                                                              <  2:1079464/149‑1 (MQ=255)
                                                                          gCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTgcg                                                                       <  1:1032649/149‑1 (MQ=255)
                                                                           cATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGACAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCAGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATCTTCTgctg                                                                      >  2:369852/1‑147 (MQ=255)
                                                                                     ccGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGGTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTCCTGCGGCGACGCTGAt                                                            <  1:991877/149‑1 (MQ=255)
                                                                                       gcgcATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAg                                                          >  2:610509/1‑149 (MQ=255)
                                                                                           aTCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGt                                                      >  1:50270/1‑149 (MQ=255)
                                                                                            tctcCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGAGGCTGATTTTACGGGCACAGGGACCTTCCCGAGACAGTTGATGTTCTGCGGCGACGCTGAGAGCGGTa                                                     >  1:769875/1‑149 (MQ=255)
                                                                                            tctcCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCGTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTa                                                     <  1:1063332/149‑1 (MQ=255)
                                                                                             ctcCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGGCGTAGTAGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAg                                                    <  1:194460/149‑1 (MQ=255)
                                                                                                cGGTGGCATCTCTCATACGCGTTCGAAATTGCGCTTATCCATGGTCAAAATCTGGTCGGAGTTGCTACACAGACGTGGGCTGATTTGAGGGGCTCAGTGACCTTCCTGAGACAGTTGATTTTCTGCGGCGACGCTGATAGCGGTAGGAt                                                 <  2:299189/149‑1 (MQ=255)
                                                                                                cGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGAt                                                 <  1:497811/149‑1 (MQ=255)
                                                                                                 ggTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGCCACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATc                                                >  1:725789/1‑149 (MQ=255)
                                                                                                 ggTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATc                                                >  1:976032/1‑149 (MQ=255)
                                                                                                      cATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAgcgc                                           >  1:109271/1‑149 (MQ=255)
                                                                                                            gCATCAGCCTTCGATATCGCGCTTATCCATGCTCAAAATCCGGTCGTCGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGGGACGCTGATAGCGGTAGGATCAGCGCCCTTAc                                     <  1:114084/149‑1 (MQ=255)
                                                                                                            gCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGCCGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTAc                                     >  2:134630/1‑149 (MQ=255)
                                                                                                            gCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGCCGCTGATAGCGGTAGGATCAGCGCCCTTAc                                     >  1:916563/1‑149 (MQ=255)
                                                                                                            gCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTAc                                     <  1:1010965/149‑1 (MQ=255)
                                                                                                             cATATGCGTTCGTTATGGCGCTTTTCCATGGTAAAAATCAGGTCGTATGTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTAcc                                    <  2:1099371/149‑1 (MQ=255)
                                                                                                                aaGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGAc                                 >  2:973650/1‑149 (MQ=255)
                                                                                                                aaGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGGTGAAGTTCTGCGGCGACGCTGATATCGGAAGGATCAGCGCCCTTACCGAc                                 >  1:691833/1‑149 (MQ=255)
                                                                                                                  gCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGCCCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCCGCGCCCTTACCGACCa                               >  2:874917/1‑149 (MQ=255)
                                                                                                                            tGGCGCTTTACCATGGTCAAAATCAGGGCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGTTAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAg                     <  1:62760/149‑1 (MQ=255)
                                                                                                                             ggCGCTTTTCCATGGTCAAAATCGGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAGt                    <  1:803844/149‑1 (MQ=255)
                                                                                                                             ggCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACTGACCAGCGCGCCGAGt                    >  1:298639/1‑149 (MQ=255)
                                                                                                                                gcTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAGTCCa                 >  2:474814/1‑149 (MQ=255)
                                                                                                                                 cTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGCCCAGCGCGCCGAGTCCAg                >  2:431820/1‑149 (MQ=255)
                                                                                                                                 cTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAGTCCAg                <  1:610509/149‑1 (MQ=255)
                                                                                                                                      ccATGGTCAAAATCGGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCACCGAGTCCAGCGGAc           >  2:1023229/1‑149 (MQ=255)
                                                                                                                                        aTGGTCAAAATCAGGTCGTAGTTGCGTCACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGTTAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAGTCAAGCGGACTc         >  1:788806/1‑149 (MQ=255)
                                                                                                                                        aTGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTCCCGACCAGCGCGCCGAGTCCAGCGGACTc         >  1:1114102/1‑149 (MQ=255)
                                                                                                                                        aTGGTCAAAATCAGGTCGTAGTTGCGAAACAGACGGCGGCTTATTTGACGGGCACAGTGACCTTCCAGATACTGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGCCCAGCGCGCCGAGTCCAGCGGACTc         >  1:154371/1‑149 (MQ=255)
                                                                                                                                           gTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCAGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCCGCGCGCCGAGTCCAGCGGACTCCAc      >  2:676104/1‑149 (MQ=255)
                                                                                                                                            tCAAAATCAGGTCGTAGTTGCGCCACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAGTCCAGCGGACTCCACt     >  1:98686/1‑149 (MQ=255)
                                                                                                                                              aaaaTCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTGGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAGTCCAGCGGACTCCACttt   <  2:220619/149‑1 (MQ=255)
                                                                                                                                               aaaTCAGTTCGACGGTGCGACACAGAGGGCGTCTGATTTGAGGGGCTCAGTGACCTTCCTGAGTCAGTTGATGTCTTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAGTCCAGCGGTCTCCACTTTc  <  2:259072/149‑1 (MQ=255)
                                                                                                                                               aaaTCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGAGCAGCGCCCTTACCGACCAGCGCGCCGAGTCCAGCGGACTCCACTTTc  <  2:11968/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCGTAGTTGCGACACAGACGGCGGCTGATTTGACGGGCACAGTGACCTTCCAGAGACAGTTGATGTTCTGCGGCGACGCTGATAGCGGTAGGATCAGCGCCCTTACCGACCAGCGCGCCGAGTCCAGCGGACTCCACTTTC  >  NC_000913/2135718‑2136009

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑