Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,584 A→G 5.2% D108G (GAC→GGC) ‡ mcbR → colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,5840AG 5.2% 250.5 / 3.7 77D108G (GAC→GGC) ‡mcbRcolanic acid and biofilm gene transcriptional regulator, MqsR‑controlled
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-02

GTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAG  >  NC_000913/1520436‑1520729
                                                                                                                                                    |                                                                                                                                                 
gTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGa                                                                                                                                                   <  3:239805/149‑1 (MQ=255)
   tCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCt                                                                                                                                                <  1:503844/149‑1 (MQ=255)
       tGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCg                                                                                                                                            >  4:514027/1‑149 (MQ=255)
             cgcTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTc                                                                                                                                      <  2:181868/149‑1 (MQ=255)
               cTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCa                                                                                                                                    <  2:776352/149‑1 (MQ=255)
                        gCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCt                                                                                                                           >  3:443117/1‑149 (MQ=255)
                        gCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCt                                                                                                                           <  1:582335/149‑1 (MQ=255)
                         cACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCtg                                                                                                                          <  4:22086/149‑2 (MQ=255)
                          aCCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCg                                                                                                                         >  3:532645/1‑149 (MQ=255)
                                acaAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAAt                                                                                                                   <  1:230277/149‑1 (MQ=255)
                                 caAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGaa                                                                                                                                         <  4:783595/126‑1 (MQ=255)
                                 caAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGaa                                                                                                                                         >  3:783595/1‑126 (MQ=255)
                                  aaGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTa                                                                                                                 >  3:787566/1‑149 (MQ=255)
                                        ttACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAcagcag                                                                                                           >  4:586560/1‑149 (MQ=255)
                                         tACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAgg                                                                                                          >  4:418163/1‑149 (MQ=255)
                                            aGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGgcgc                                                                                                       >  4:236834/1‑149 (MQ=255)
                                            aGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGgcgc                                                                                                       >  1:486265/1‑149 (MQ=255)
                                                  ggAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGaaa                                                                                                 >  4:142813/1‑149 (MQ=255)
                                                   gAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGaaaa                                                                                                <  4:106019/149‑1 (MQ=255)
                                                   gAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGaaaa                                                                                                <  4:443117/149‑1 (MQ=255)
                                                   gAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGaaaa                                                                                                <  2:403835/149‑1 (MQ=255)
                                                   gAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGaaaa                                                                                                >  2:640668/1‑149 (MQ=255)
                                                     aGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCCTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAgg                                                                                              <  4:537893/149‑1 (MQ=255)
                                                       tGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAca                                                                                                               <  3:315004/130‑1 (MQ=255)
                                                       tGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTAca                                                                                                               >  4:315004/1‑130 (MQ=255)
                                                        ggggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCaa                                                                                                     <  4:647401/139‑1 (MQ=255)
                                                        ggggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCaa                                                                                                     >  3:647401/1‑139 (MQ=255)
                                                         gggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGa                                                                                          >  2:521429/1‑149 (MQ=255)
                                                         gggAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTCGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGTATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGa                                                                                          >  4:675354/1‑149 (MQ=255)
                                                           gAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACa                                                                                        >  2:735918/1‑149 (MQ=255)
                                                            aaaCGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGGTCTTGCggg                                                                                                                                          <  3:404006/98‑2 (MQ=255)
                                                            aaaCGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGGTCTTGCggg                                                                                                                                          >  4:404006/1‑97 (MQ=255)
                                                             aaCGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATg                                                                                      <  1:563474/149‑1 (MQ=255)
                                                                 tCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAAc                                                                                  >  2:525940/1‑149 (MQ=255)
                                                                 tCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAAc                                                                                  >  4:649188/1‑149 (MQ=255)
                                                                   aaTTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACaa                                                                                <  4:787566/149‑1 (MQ=255)
                                                                    aTTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACaaa                                                                               <  3:611987/149‑1 (MQ=255)
                                                                       aGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCa                                                                            <  2:288036/149‑1 (MQ=255)
                                                                         aCGAAATCAATCGGAACCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCAtt                                                                          <  1:409439/149‑1 (MQ=255)
                                                                            aaaTCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAAt                                                                       <  3:456/149‑1 (MQ=255)
                                                                                 aatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                        >  2:255161/1‑143 (MQ=255)
                                                                                 aatcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTaa                                                                        <  1:255161/143‑1 (MQ=255)
                                                                                  atcGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAAc                                                                 >  3:700927/1‑149 (MQ=255)
                                                                                     ggATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAgg                                                              >  2:356648/1‑149 (MQ=255)
                                                                                     ggATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAgg                                                              <  1:735918/149‑1 (MQ=255)
                                                                                      gATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGc                                                             >  3:383589/1‑149 (MQ=255)
                                                                                         ccGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTAt                                                          >  2:183431/1‑149 (MQ=255)
                                                                                             tACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCg                                                      >  1:369528/1‑149 (MQ=255)
                                                                                               cGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCt                                                    <  3:512145/149‑1 (MQ=255)
                                                                                                  aCTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAg                                                 >  3:578736/1‑149 (MQ=255)
                                                                                                   cTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGc                                                <  2:91546/149‑1 (MQ=255)
                                                                                                    tGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCg                                               >  1:97625/1‑149 (MQ=255)
                                                                                                       aaTTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAtt                                            <  4:383589/149‑1 (MQ=255)
                                                                                                         ttAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTa                                          <  2:97625/149‑1 (MQ=255)
                                                                                                           aaTGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACATGCTATTTCGCTTAGCGATTTatc                                        <  4:573858/149‑1 (MQ=255)
                                                                                                           aaTGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatc                                        >  3:148665/1‑149 (MQ=255)
                                                                                                            aTGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatca                                       <  3:236834/149‑1 (MQ=255)
                                                                                                            aTGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatca                                       <  3:514027/149‑1 (MQ=255)
                                                                                                                     gCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCaaa                              <  1:809805/149‑1 (MQ=255)
                                                                                                                      ctctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAat                             >  1:766171/1‑149 (MQ=255)
                                                                                                                       tctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatcatc                                     <  4:158851/140‑1 (MQ=255)
                                                                                                                       tctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTatcatc                                     >  3:158851/1‑140 (MQ=255)
                                                                                                                       tctGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAata                            >  3:197104/1‑149 (MQ=255)
                                                                                                                          ggCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGc                         <  1:183431/149‑1 (MQ=255)
                                                                                                                             tGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCa                      >  2:184082/1‑149 (MQ=255)
                                                                                                                             tGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCa                      >  4:341205/1‑149 (MQ=255)
                                                                                                                                    aaCCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCtgtg               <  4:532645/149‑1 (MQ=255)
                                                                                                                                        tCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTgag           >  1:6085/1‑149 (MQ=255)
                                                                                                                                         cACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTgaga          <  4:290526/149‑1 (MQ=255)
                                                                                                                                         cACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTgaga          >  3:249444/1‑149 (MQ=255)
                                                                                                                                          aCCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgat         >  3:391520/1‑149 (MQ=255)
                                                                                                                                           ttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAtt                                            <  4:139193/111‑1 (MQ=255)
                                                                                                                                           ttccGCAAGGTCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGAtt                                            >  3:139193/3‑113 (MQ=255)
                                                                                                                                            cccGCAATACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCa                                                                            <  4:481467/80‑1 (MQ=255)
                                                                                                                                            cccGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCa                                                                            >  3:481467/1‑80 (MQ=255)
                                                                                                                                             ccGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAgatgat      <  4:148665/149‑1 (MQ=255)
                                                                                                                                                cAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGa   <  3:418163/149‑1 (MQ=255)
                                                                                                                                                 aaGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAGGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAg  <  3:47755/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTTTCGGTGAATGCGCTTTCTGTCGCACCTGCACAAGCATTTACAGTTCCGGAAGTGGGGAAACGTCAATTAGACGAAATCAATCGGATCCGCTACGAACTGGAATTAATGGCAGTTGCTCTGGCTGTTGAAAACCTCACCCCGCAAGACCTTGCGGAACTCCAGGAATTGCTCGAGAAATTACAGCAGGCGCAAGAAAAGGGCGACATGGAACAAATCATTAATGTAAACAGGCTATTTCGCTTAGCGATTTATCATCGCTCAAATATGCCCATCCTGTGTGAGATGATTGAG  >  NC_000913/1520436‑1520729

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