Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,121,588:1 +C 9.9% intergenic (‑34/+21) 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,5881.C 9.9% 486.3 / 50.4 121intergenic (‑34/+21)wcaJ/cpsGcolanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base . (57/52);  new base C (7/5);  total (64/57)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 131 total aligned reads displayed.

CTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCG  >  NC_000913/2121440‑2121735
                                                                                                                                                     |                                                                                                                                                   
cTTCGCAAACCAGCCATAGTCCGGCAAGCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTGATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTg                                                                                                                                                      >  3:974512/1‑149 (MQ=255)
    gCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGt                                                                                                                                                 >  4:977287/1‑149 (MQ=255)
          cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGa                                                                                                                                           >  1:1387344/1‑149 (MQ=255)
             ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGa                                                                                                                                        >  2:301097/1‑149 (MQ=255)
             ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGa                                                                                                                                        >  1:348696/1‑149 (MQ=255)
              cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAggg                                                                                                                                            <  3:17313/144‑1 (MQ=255)
              cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAggg                                                                                                                                            >  4:17313/1‑144 (MQ=255)
               aTAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                      >  3:1128364/1‑149 (MQ=255)
               aTAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                      >  3:591701/1‑149 (MQ=255)
                       gCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTaa                                                                                                                                               >  2:1132684/1‑132 (MQ=255)
                       gCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGAGTaa                                                                                                                                               <  1:1132684/132‑1 (MQ=255)
                        cAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTAc                                                                                                                             <  1:716026/149‑1 (MQ=255)
                              tgatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTc                                                                                                                       <  2:1387344/149‑1 (MQ=255)
                               gatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTAc                                                                                                                             <  4:1224511/142‑1 (MQ=255)
                               gatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCa                                                                                                                      <  3:977287/149‑1 (MQ=255)
                                        atCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                    <  4:674466/110‑1 (MQ=255)
                                        atCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                    >  3:674466/1‑110 (MQ=255)
                                        atCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTca                                                                                                             <  1:301097/149‑1 (MQ=255)
                                         tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgggg                                                                                                                                                  <  2:96659/111‑1 (MQ=255)
                                         tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgggg                                                                                                                                                  >  1:96659/1‑111 (MQ=255)
                                           tGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAgggg                                                                                                                                           >  4:392171/1‑115 (MQ=255)
                                           tGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAgggg                                                                                                                                           <  3:392171/116‑2 (MQ=255)
                                             aaaaGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                        <  1:1014811/149‑1 (MQ=255)
                                                aGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGcca                                                                                                                  >  3:198010/1‑134 (MQ=255)
                                                aGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCaa                                                                                                                  <  4:198010/136‑1 (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)
                                                              gagaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTc                                                                                       >  3:513940/1‑149 (MQ=255)
                                                                gaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGa                                                                                                                                        >  3:272691/1‑98 (MQ=255)
                                                                gaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGa                                                                                                                                        <  4:272691/98‑1 (MQ=255)
                                                                 aTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGaa                                                                                                                                          >  4:482968/1‑95 (MQ=255)
                                                                     aCGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGc                                                                                <  2:547131/149‑1 (MQ=255)
                                                                     aCGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGc                                                                                <  4:591701/149‑1 (MQ=255)
                                                                         tGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTcag                                                                                                            <  3:379905/117‑1 (MQ=255)
                                                                          gCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGcac                                                                           <  4:391265/149‑1 (MQ=255)
                                                                          gCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAggg                                                                                                                                            <  1:1435831/84‑1 (MQ=38)
                                                                            aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                      <  1:3988/88‑1 (MQ=255)
                                                                            aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                      >  2:1411557/1‑88 (MQ=255)
                                                                            aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                      >  2:3988/1‑88 (MQ=255)
                                                                            aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgtt                                                                                               <  1:1364116/127‑1 (MQ=255)
                                                                            aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgtt                                                                                               >  2:1364116/1‑127 (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                                                                                                                                                   >  4:848852/1‑53 (MQ=255)
                                                                                                  ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTggg                                                                                                                                                   <  3:848852/53‑1 (MQ=255)
                                                                                                  ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGc                                                                                                                    >  3:1141724/1‑84 (MQ=255)
                                                                                                  ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGc                                                                                                                    <  4:1141724/84‑1 (MQ=255)
                                                                                                  ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAg                                                   <  4:974512/149‑1 (MQ=255)
                                                                                                      ttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                        >  2:653843/1‑92 (MQ=255)
                                                                                                          aTTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcacca                                           >  3:703585/1‑149 (MQ=255)
                                                                                                               tCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgc                                                                                                  <  4:1335360/89‑1 (MQ=255)
                                                                                                               tCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgc                                                                                                  >  3:1335360/1‑89 (MQ=255)
                                                                                                                 aTCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                    >  1:718631/1‑149 (MQ=255)
                                                                                                                  tCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCtt                                                                                        <  3:509076/96‑1 (MQ=255)
                                                                                                                       gttCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTAtt                              <  4:703585/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)
                                                                                                                          ccTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGa                           <  2:718631/149‑1 (MQ=255)
                                                                                                                            tGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcac                                             <  4:148429/129‑1 (MQ=255)
                                                                                                                            tGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcac                                             >  3:148429/1‑129 (MQ=255)
                                                                                                                              ttattaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTg                       >  4:298676/1‑149 (MQ=255)
                                                                                                                                attaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTa                                <  4:1081320/138‑1 (MQ=255)
                                                                                                                                attaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTa                                >  3:1081320/1‑138 (MQ=255)
                                                                                                                                 ttaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTgcgc                    >  1:361580/1‑149 (MQ=255)
                                                                                                                                 ttaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTgcgc                    <  3:298676/149‑1 (MQ=255)
                                                                                                                                 atcttcccTTACCCCACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                    >  3:1352626/6‑133 (MQ=255)
                                                                                                                                 atcttcccTTACCCCACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                    <  4:1352626/128‑1 (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)
                                                                                                                                   aGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTgc                      <  2:808500/145‑1 (MQ=255)
                                                                                                                                   aGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTgc                      >  1:808500/1‑145 (MQ=255)
                                                                                                                                   aGCCACTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTg                             <  4:619521/138‑1 (MQ=255)
                                                                                                                                    ttcccTTACCCCACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg                 >  3:1278939/3‑149 (MQ=255)
                                                                                                                                    gCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttc                                                                                                    >  3:916977/1‑66 (MQ=255)
                                                                                                                                    gCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttc                                                                                                    <  4:916977/66‑1 (MQ=255)
                                                                                                                                     ccccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt                >  1:106976/1‑149 (MQ=255)
                                                                                                                                      cccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGtt               >  3:886001/1‑149 (MQ=255)
                                                                                                                                       ccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa              >  1:1063437/1‑149 (MQ=255)
                                                                                                                                        cTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                                   <  2:311053/127‑1 (MQ=255)
                                                                                                                                        cTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                                   >  1:311053/1‑127 (MQ=255)
                                                                                                                                         ttACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                     <  2:410990/124‑1 (MQ=255)
                                                                                                                                          tACCCCCAGTGGGGTAAGGGAAGATCCGACATTGCTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCAACGGTTCGGTATTg                             >  1:415387/1‑131 (MQ=255)
                                                                                                                                             ccGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCa                   <  2:305708/138‑1 (MQ=255)
                                                                                                                                             ccGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCa                   >  1:305708/1‑138 (MQ=255)
                                                                                                                                             ccGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc        >  1:195409/1‑149 (MQ=255)
                                                                                                                                              cGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGcc       <  2:106976/149‑1 (MQ=255)
                                                                                                                                               gCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTAAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGTGCGATTCCACCTTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGTAGGTTCAAGCGCCa      >  4:902597/1‑149 (MQ=255)
                                                                                                                                                   tGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCg  >  1:155128/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCAGTCG  >  NC_000913/2121440‑2121735

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: