Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,088,550 A→G 5.2% S436S (TCT→TCC pgaB ← poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine (PGA) N‑deacetylase outer membrane export lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,088,5500AG 5.2% 241.9 / 2.2 76S436S (TCT→TCCpgaBpoly‑beta‑1,6‑N‑acetyl‑D‑glucosamine (PGA) N‑deacetylase outer membrane export lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base G (2/2);  total (37/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 = 9.87e-01

GATAAGCCGTGATAGCCGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTACACCTGAGCGGGTACGTAATT  >  NC_000913/1088405‑1088695
                                                                                                                                                 |                                                                                                                                                 
gATAAGCCGTGATAGCCGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGga                                                                                                                                                   >  6:15342/1‑146 (MQ=255)
 aTACGCCGTGATAGCCGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGgagaga                                                                                                                                               <  1:96127/149‑1 (MQ=255)
     gCCGTGATAGCCGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACg                                                                                                                                           >  8:99491/1‑149 (MQ=255)
     gCCGTGATAGCCGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACg                                                                                                                                           <  8:36086/149‑1 (MQ=255)
             aGCCGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGgagagaga                                                                                                                                             >  8:333255/1‑139 (MQ=255)
             aGCCGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGgagagaga                                                                                                                                             <  7:333255/139‑1 (MQ=255)
              gCCGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATAtt                                                                                                                                   <  5:15342/148‑1 (MQ=255)
                cGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGGGAGAGACGGTGATATTGtt                                                                                                                                >  1:250076/1‑149 (MQ=255)
                cGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGGGAGAGACGGTGATATTGtt                                                                                                                                >  5:420054/1‑149 (MQ=255)
                  gTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCa                                                                                                                              <  8:428995/149‑1 (MQ=255)
                      aCTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGAt                                                                                                                          >  5:240578/1‑149 (MQ=255)
                        tGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGgagag                                                                                                                                                >  4:16982/1‑125 (MQ=255)
                        tGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGgagag                                                                                                                                                <  3:16982/125‑1 (MQ=255)
                        tGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATg                                                                                                                         >  6:442641/1‑148 (MQ=255)
                           cATCTTCATACTCTGAAAGCAAAGCATCATCGTGGAACAATATGCCCTCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAAtt                                                                                                                     <  4:284212/149‑1 (MQ=255)
                             tCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGga                                                                                                                                                   >  7:152865/1‑117 (MQ=255)
                             tCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGga                                                                                                                                                   <  8:152865/117‑1 (MQ=255)
                                 cATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTtgt                                                                                                                  <  1:234974/146‑1 (MQ=255)
                                 cATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTACAGGATGAATTtgt                                                                                                                  <  3:286015/146‑1 (MQ=255)
                                     aTCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGGGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCtttt                                                                                                            <  6:420054/148‑1 (MQ=255)
                                     aTCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGGGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCtttt                                                                                                            <  2:250076/148‑1 (MQ=255)
                                      tCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCtttt                                                                                                            >  5:335933/1‑147 (MQ=255)
                                      tCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCtttt                                                                                                            <  6:335933/147‑1 (MQ=255)
                                       cTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTtc                                                                                                         >  2:57274/1‑149 (MQ=255)
                                         gAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTtctc                                                                                                       >  2:203324/1‑149 (MQ=255)
                                           aagcaaagcaTCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTccc                                                                                                     <  7:99491/149‑1 (MQ=255)
                                            agcaaagcaTCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTcccc                                                                                                    <  6:255255/149‑1 (MQ=255)
                                               aaagcaTCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGATTTTGTGCTTTTTTCTCCCCTGt                                                                                                 >  7:424881/1‑149 (MQ=255)
                                               aaagcaTAATCGTGGAACAATATGACATAAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTTACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGTATGAATTTGTGCTTTTTTATCACCTGt                                                                                                 >  6:95923/1‑149 (MQ=255)
                                                aagcaTCATAGTGGAACAATATGCCATCAAAAGGATCATGTCCGGCAAGTTCTTCATATTACATGCCAAGTTGTGCTCTGACTCTGTGATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGATTTTGTGCTTTTTTCTCCCCTGtt                                                                                                <  7:169606/149‑1 (MQ=255)
                                                    atcatcGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGt                                                                                             >  2:132001/1‑148 (MQ=255)
                                                       atcGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAg                                                                                          >  8:357436/1‑148 (MQ=255)
                                                          gTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTaa                                                                                           <  4:135246/144‑1 (MQ=255)
                                                          gTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTaa                                                                                           >  3:135246/1‑144 (MQ=255)
                                                          gTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTAtt                                                                                      >  8:78609/1‑149 (MQ=255)
                                                                cAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCTACTTGTGCTCTGACTCTGTCATCGAAAGgagaga                                                                                                                                               <  6:57770/86‑1 (MQ=255)
                                                                cAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGgagaga                                                                                                                                               <  6:64211/86‑1 (MQ=255)
                                                                cAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGgagaga                                                                                                                                               >  5:57770/1‑86 (MQ=255)
                                                                cAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGgagaga                                                                                                                                               >  5:64211/1‑86 (MQ=255)
                                                                  atatGCCATCAAAAGCAGCATGTCCGGCAAGAACTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTATTCTCCCCTGTAGGTAAGTATTTTAATcgc                                                                              >  6:331660/1‑149 (MQ=255)
                                                                      gCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGtt                                                                           >  7:78623/1‑148 (MQ=255)
                                                                        cATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGAATCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAAt                                                                        >  4:405213/1‑149 (MQ=255)
                                                                         aTCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAAtg                                                                       >  2:92208/1‑149 (MQ=255)
                                                                          tCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAAtgt                                                                      >  3:342986/1‑149 (MQ=255)
                                                                          tCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAAtgt                                                                      >  8:196552/1‑149 (MQ=255)
                                                                            aaaaGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAAtgtg                                                                     <  6:240578/148‑1 (MQ=255)
                                                                            aaaaGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTgg                                                                    <  8:78623/149‑1 (MQ=255)
                                                                            aaaaGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTgg                                                                    <  1:92208/149‑1 (MQ=255)
                                                                             aaaGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTgg                                                                    >  7:430496/1‑148 (MQ=255)
                                                                              aaGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGa                                                                  <  5:442641/149‑1 (MQ=255)
                                                                                gcagcaTGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATc                                                                >  7:379310/1‑149 (MQ=255)
                                                                                 cagcaTGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATc                                                                >  8:84605/1‑148 (MQ=255)
                                                                                        tCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATccc                                                        <  1:57274/149‑1 (MQ=255)
                                                                                             cAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCtt                                                   <  6:15101/149‑1 (MQ=255)
                                                                                             cAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCtt                                                   <  4:342986/149‑1 (MQ=255)
                                                                                                    ttCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCttt                                                                                                             >  1:50509/1‑84 (MQ=255)
                                                                                                    ttCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCttt                                                                                                             <  2:50509/84‑1 (MQ=255)
                                                                                                     tCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGt                                                                                                 >  3:420458/1‑95 (MQ=255)
                                                                                                     tCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGt                                                                                                 <  4:420458/95‑1 (MQ=255)
                                                                                                     tCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCATGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATAACgg                                           >  6:251068/1‑149 (MQ=255)
                                                                                                      cATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATAcc                                             <  1:203324/146‑1 (MQ=255)
                                                                                                      cATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATAcc                                             <  5:331660/146‑1 (MQ=255)
                                                                                                       atatAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCa                                         <  7:78609/149‑1 (MQ=255)
                                                                                                        tataACATGCCAACTTGTGCTCTGACTCTGTCATCGAACGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCAt                                        <  6:86484/149‑1 (MQ=255)
                                                                                                          taACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATc                                       >  3:38042/1‑148 (MQ=255)
                                                                                                           aaCATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGc                                                     >  5:408223/1‑133 (MQ=255)
                                                                                                           aaCATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGc                                                     <  6:408223/133‑1 (MQ=255)
                                                                                                                    aCTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTTGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGAt                            <  5:95923/149‑1 (MQ=255)
                                                                                                                     cTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATg                           <  4:38042/149‑1 (MQ=255)
                                                                                                                     cTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATg                           >  2:44218/1‑149 (MQ=255)
                                                                                                                      ttGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGt                          >  8:70609/1‑149 (MQ=255)
                                                                                                                        gtgCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGtt                         <  8:379310/148‑1 (MQ=255)
                                                                                                                        gtgCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGtt                         <  8:430496/148‑1 (MQ=255)
                                                                                                                           ctctGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTaca                     >  7:93483/1‑149 (MQ=255)
                                                                                                                             ctGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGt                          <  5:205461/142‑1 (MQ=255)
                                                                                                                             ctGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGt                          >  6:205461/1‑142 (MQ=255)
                                                                                                                                    tGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTACACCTGAGCgg            >  4:251256/1‑149 (MQ=255)
                                                                                                                                      tCATCGAAAGGAGAGAGACTGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTac                      <  8:219098/137‑1 (MQ=255)
                                                                                                                                      tCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTac                      >  7:219098/1‑137 (MQ=255)
                                                                                                                                      tCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTACACCTGAGCggg           <  7:357436/148‑1 (MQ=255)
                                                                                                                                      tCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTACACCTGAGCggg           >  6:411790/1‑148 (MQ=255)
                                                                                                                                         tCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGt                          >  2:174072/1‑130 (MQ=255)
                                                                                                                                         tCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGt                          <  1:174072/130‑1 (MQ=255)
                                                                                                                                         tCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTac                      >  3:10008/1‑134 (MQ=255)
                                                                                                                                         tCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTac                      <  4:10008/134‑1 (MQ=255)
                                                                                                                                           gAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGCTCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTACACCTGAGCGGGTACGTa     <  1:132001/149‑1 (MQ=255)
                                                                                                                                           gAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCAc                                    <  6:438862/118‑1 (MQ=255)
                                                                                                                                              aGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTACACCTGAGCGGGTACGTAAtt  <  7:70609/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GATAAGCCGTGATAGCCGGTGCACTGGCATCTTCATAATCTGAAAGCAAAGCATCATCGTGGAACAATATGCCATCAAAAGCAGCATGTCCGGCAAGATCTTCATATAACATGCCAACTTGTGCTCTGACTCTGTCATCGAAAGGAGAGAGACGGTGATATTGTTCAGGATGAATTTGTGCTTTTTTCTCCCCTGTTGGTAAGTATTTTACTCGCGTTAATGTGGGATCTAAATCCCAGCTTAATACCGGCATCCACGCATAGATGTTTACACCTGAGCGGGTACGTAATT  >  NC_000913/1088405‑1088695

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