Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,121,583 Δ1 bp 9.3% intergenic (‑29/+26) wcaJ ← / ← cpsG colanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,121,5830G. 9.3% 437.0 / 39.2 129intergenic (‑29/+26)wcaJ/cpsGcolanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (60/56);  new base . (7/5);  total (68/61)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 130 total aligned reads displayed.

CTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCA  >  NC_000913/2121440‑2121731
                                                                                                                                               |                                                                                                                                                     
cTTCGCAAACCAGCCATAGTCCGGCAAGCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTGATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg                                                                                                                                                  >  3:974512/1‑149 (MQ=255)
 ttCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                >  2:716026/1‑149 (MQ=255)
    gCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGt                                                                                                                                             >  4:977287/1‑149 (MQ=255)
          cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGa                                                                                                                                       >  1:1387344/1‑149 (MQ=255)
             ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGa                                                                                                                                    <  1:936220/149‑1 (MQ=255)
             ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGa                                                                                                                                    >  1:348696/1‑149 (MQ=255)
              cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGAt                                                                                                                                   <  2:348696/149‑1 (MQ=255)
              cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAggg                                                                                                                                        >  4:17313/1‑144 (MQ=255)
               aTAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                  >  3:591701/1‑149 (MQ=255)
               aTAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                  >  3:1128364/1‑149 (MQ=255)
                       gCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTaa                                                                                                                                           >  2:1132684/1‑132 (MQ=255)
                       gCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGAGTaa                                                                                                                                           <  1:1132684/132‑1 (MQ=255)
                              tgatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTc                                                                                                                   <  2:1387344/149‑1 (MQ=255)
                               gatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTAc                                                                                                                         <  4:1224511/142‑1 (MQ=255)
                               gatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTAc                                                                                                                         >  3:1224511/1‑142 (MQ=255)
                                        atCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                >  3:674466/1‑110 (MQ=255)
                                        atCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                <  4:674466/110‑1 (MQ=255)
                                         tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgggg                                                                                                                                              <  2:96659/111‑1 (MQ=255)
                                         tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgggg                                                                                                                                              >  1:96659/1‑111 (MQ=255)
                                           tGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAgggg                                                                                                                                       >  4:392171/1‑115 (MQ=255)
                                             aaaaGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                    <  1:1014811/149‑1 (MQ=255)
                                                aGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGcca                                                                                                              >  3:198010/1‑134 (MQ=255)
                                                         ccATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttcgc                                                                                        <  4:1128364/149‑1 (MQ=255)
                                                           aTAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTggg                                                                                                                                               >  1:530306/1‑92 (MQ=255)
                                                           aTAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTggg                                                                                                                                               <  2:530306/92‑1 (MQ=255)
                                                           aTAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGa                                                                                                                               <  3:660778/108‑1 (MQ=255)
                                                           aTAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGa                                                                                                                               >  4:660778/1‑108 (MQ=255)
                                                           aTAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCgcgc                                                                                      >  3:391265/1‑149 (MQ=255)
                                                            tAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTc                                                                                                                   >  2:1204128/1‑119 (MQ=255)
                                                            tAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTc                                                                                                                   <  1:1204128/119‑1 (MQ=255)
                                                            tAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCt                                                                                     >  1:547131/1‑149 (MQ=255)
                                                            tAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCt                                                                                     >  3:338823/1‑149 (MQ=255)
                                                                gaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGa                                                                                                                                    <  4:272691/98‑1 (MQ=255)
                                                                gaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGa                                                                                                                                    >  3:272691/1‑98 (MQ=255)
                                                                 aTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGaa                                                                                                                                      <  3:482968/95‑1 (MQ=255)
                                                                     aCGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGc                                                                            <  4:591701/149‑1 (MQ=255)
                                                                        aTGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGt                                                                                                                                             >  3:206203/1‑81 (MQ=255)
                                                                        aTGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGt                                                                                                                                             <  4:206203/81‑1 (MQ=255)
                                                                         tGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTcag                                                                                                        <  3:379905/117‑1 (MQ=255)
                                                                          gCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGcac                                                                       <  4:391265/149‑1 (MQ=255)
                                                                          gCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAggg                                                                                                                                        >  2:1435831/1‑84 (MQ=38)
                                                                          gCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAggg                                                                                                                                        <  1:1435831/84‑1 (MQ=38)
                                                                            aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                  <  1:3988/88‑1 (MQ=255)
                                                                            aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                  >  2:1411557/1‑88 (MQ=255)
                                                                            aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgtt                                                                                           <  1:1364116/127‑1 (MQ=255)
                                                                             ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACt                                                                                                                        <  3:103838/97‑1 (MQ=255)
                                                                             ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACt                                                                                                                        >  4:103838/1‑97 (MQ=255)
                                                                             ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTCGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGcac                                                                       <  3:930059/146‑1 (MQ=255)
                                                                             ttGGTATTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGCAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGcac                                                                       >  4:930059/1‑146 (MQ=255)
                                                                                                 cTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCa                                                <  4:513940/149‑1 (MQ=255)
                                                                                                  ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTggg                                                                                                                                               <  3:848852/53‑1 (MQ=255)
                                                                                                  ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTggg                                                                                                                                               >  4:848852/1‑53 (MQ=255)
                                                                                                  ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAg                                               <  4:974512/149‑1 (MQ=255)
                                                                                                      ttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                    <  1:653843/92‑1 (MQ=255)
                                                                                                      ttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                    >  2:653843/1‑92 (MQ=255)
                                                                                                          aTTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcacca                                       <  3:752057/149‑1 (MQ=255)
                                                                                                               tCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgc                                                                                              <  4:1335360/89‑1 (MQ=255)
                                                                                                                 aTCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                >  1:718631/1‑149 (MQ=255)
                                                                                                                  tCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCtt                                                                                    <  3:509076/96‑1 (MQ=255)
                                                                                                                  tCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCtt                                                                                    >  4:509076/1‑96 (MQ=255)
                                                                                                                       gttCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTAtt                          <  4:703585/149‑1 (MQ=255)
                                                                                                                         tCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTgg                        >  4:1199538/1‑149 (MQ=255)
                                                                                                                          ccTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGa                       <  2:718631/149‑1 (MQ=255)
                                                                                                                          ccTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGa                       <  3:1199538/149‑1 (MQ=255)
                                                                                                                          ccTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGa                       >  3:293490/1‑149 (MQ=255)
                                                                                                                          ccTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGa                       <  4:338823/149‑1 (MQ=255)
                                                                                                                            tGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcac                                         >  3:148429/1‑129 (MQ=255)
                                                                                                                            tGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcac                                         <  4:148429/129‑1 (MQ=255)
                                                                                                                            tGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                     >  3:929929/1‑149 (MQ=255)
                                                                                                                              ttattaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTg                   >  4:298676/1‑149 (MQ=255)
                                                                                                                              ttattaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTg                   >  4:658698/1‑149 (MQ=255)
                                                                                                                                attaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTa                            <  4:1081320/138‑1 (MQ=255)
                                                                                                                                attaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTa                            >  3:1081320/1‑138 (MQ=255)
                                                                                                                                 ttaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAg                                               >  3:1404496/1‑118 (MQ=255)
                                                                                                                                 ttaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCCCGATTCCACATTCAg                                               <  4:1404496/118‑1 (MQ=255)
                                                                                                                                 atcttcccTTACCCCACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                >  3:1352626/6‑133 (MQ=255)
                                                                                                                                   cttcccTTACCCCACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                 <  2:376712/127‑1 (MQ=255)
                                                                                                                                   cttcccTTACCCCACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                 >  1:376712/4‑130 (MQ=255)
                                                                                                                                   aGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGc                                                                                                                >  3:442858/1‑51 (MQ=255)
                                                                                                                                   aGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGc                                                                                                                <  4:442858/51‑1 (MQ=255)
                                                                                                                                   aGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTg                         >  3:619521/1‑138 (MQ=255)
                                                                                                                                   aGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTgc                  >  1:808500/1‑145 (MQ=255)
                                                                                                                                    ttcccTTACCCCACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg             >  3:1278939/3‑149 (MQ=255)
                                                                                                                                    gCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttc                                                                                                <  4:916977/66‑1 (MQ=255)
                                                                                                                                    gCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg             <  3:990575/149‑1 (MQ=255)
                                                                                                                                     tcccTTACCCCACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt            >  1:325597/2‑149 (MQ=255)
                                                                                                                                     ccccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt            >  1:23112/1‑149 (MQ=255)
                                                                                                                                     ccccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt            >  1:106976/1‑149 (MQ=255)
                                                                                                                                      cccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGtt           >  3:886001/1‑149 (MQ=255)
                                                                                                                                       ccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa          >  1:1063437/1‑149 (MQ=255)
                                                                                                                                        cTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACAt                                                   <  2:1203946/107‑1 (MQ=255)
                                                                                                                                        cTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                               >  1:311053/1‑127 (MQ=255)
                                                                                                                                        cTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                               <  2:311053/127‑1 (MQ=255)
                                                                                                                                         ttACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                 >  1:410990/1‑124 (MQ=255)
                                                                                                                                         ttACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                 <  2:410990/124‑1 (MQ=255)
                                                                                                                                          tACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTg                         <  2:415387/131‑1 (MQ=255)
                                                                                                                                             ccGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCa               <  2:305708/138‑1 (MQ=255)
                                                                                                                                             ccGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCa               >  1:305708/1‑138 (MQ=255)
                                                                                                                                             ccGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc    >  1:195409/1‑149 (MQ=255)
                                                                                                                                               gCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTAAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGTGCGATTCCACCTTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGTAGGTTCAAGCGCCa  >  4:902597/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCA  >  NC_000913/2121440‑2121731

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