Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 722,654 C→T 5.5% A587A (GCG→GCA kdpD ← sensory histidine kinase KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913722,6540CT 5.5% 230.6 / 2.4 73A587A (GCG→GCAkdpDsensory histidine kinase KdpD
Reads supporting (aligned to +/- strand):  ref base C (30/39);  new base T (2/2);  total (32/41)
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.75e-01

AACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAACGTGGAGGCAATAAA  >  NC_000913/722505‑722803
                                                                                                                                                     |                                                                                                                                                     
aaCAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTgccc                                                                                                                                                       >  2:699721/1‑150 (MQ=255)
  cAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTgcccgc                                                                                                                                                     >  2:1762262/1‑150 (MQ=255)
  cAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTgcccgc                                                                                                                                                     >  2:65994/1‑150 (MQ=255)
    gCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTgcccgcgc                                                                                                                                                   <  1:466334/150‑1 (MQ=255)
     cGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTgcccgcgcc                                                                                                                                                  <  1:1237973/150‑1 (MQ=255)
            gTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAATGTGTCGGTGCCCGCGCCCgcaggc                                                                                                                                           >  2:329284/1‑150 (MQ=255)
            gTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCgcaggc                                                                                                                                           >  1:795248/1‑150 (MQ=255)
              ctcCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCgcaggcag                                                                                                                                         <  2:648024/150‑1 (MQ=255)
               tcCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGcaggcagg                                                                                                                                        <  1:982353/150‑1 (MQ=255)
                ccAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCaggcaggc                                                                                                                                       >  1:518485/1‑150 (MQ=255)
                  agcagGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCt                                                                                                                                     >  2:715123/1‑150 (MQ=255)
                     agGCGCTGCTGTTCCGGGACCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTa                                                                                                                                  <  1:369639/150‑1 (MQ=255)
                       gcgcTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATc                                                                                                                                <  2:1221549/150‑1 (MQ=255)
                        cgcTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCa                                                                                                                               >  1:1437690/1‑150 (MQ=255)
                          ctgctgTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCaaa                                                                                                                             <  2:573364/150‑1 (MQ=255)
                           tgctgTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCaaaa                                                                                                                            >  2:1975099/1‑150 (MQ=255)
                            gctgTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACTGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAAc                                                                                                                           <  2:587813/150‑1 (MQ=255)
                            gctgTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAAc                                                                                                                           <  1:1762245/150‑1 (MQ=255)
                             ctgTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACt                                                                                                                          >  2:1724984/1‑150 (MQ=255)
                                  ccGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGATACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACt                                                                                                                     <  1:11460/150‑1 (MQ=255)
                                   cGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCgcaggcag                                                                                                                                         <  2:496482/129‑1 (MQ=255)
                                   cGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCgcaggcag                                                                                                                                         >  1:496482/1‑129 (MQ=255)
                                     ggATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTgcg                                                                                                                  <  2:1476193/150‑1 (MQ=255)
                                       atcatcAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGcgcg                                                                                                                >  2:1731820/1‑150 (MQ=255)
                                        tcatcaACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGa                                                                                                               >  1:73423/1‑150 (MQ=255)
                                          atcaACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCTgcgggcgcgggc                                                                                                                                           >  1:1082664/1‑120 (MQ=255)
                                          atcaACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCTgcgggcgcgggc                                                                                                                                           <  2:1082664/120‑1 (MQ=255)
                                          atcaACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATg                                                                                                             >  2:348526/1‑150 (MQ=255)
                                          atcaACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATg                                                                                                             <  1:699721/150‑1 (MQ=255)
                                          atcaACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATg                                                                                                             >  2:67057/1‑150 (MQ=255)
                                           tcaACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATgg                                                                                                            <  1:1762262/150‑1 (MQ=255)
                                           tcaACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATgg                                                                                                            <  1:65994/150‑1 (MQ=255)
                                           tcaACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATgg                                                                                                            <  1:1519015/150‑1 (MQ=255)
                                                 ggCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGt                                                                                                      >  2:1769447/1‑150 (MQ=255)
                                                  gcgcAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTc                                                                                                     <  1:715123/150‑1 (MQ=255)
                                                  gcgcAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTc                                                                                                     <  2:518485/150‑1 (MQ=255)
                                                   ggcAGATTCCGCGTCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGTTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTcc                                                                                                    <  2:925936/149‑1 (MQ=255)
                                                          tCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTAACAAAACTCCACTGCGCGATGGCATCGTCCCACGGCg                                                                                             >  2:82787/1‑150 (MQ=255)
                                                            cccGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTc                                                                                           <  2:576054/150‑1 (MQ=255)
                                                                gCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTc                                                                                       >  1:1772067/1‑150 (MQ=255)
                                                                 cTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAAGGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTGCCACGGCGTCATTcg                                                                                      <  2:968502/150‑2 (MQ=255)
                                                                   ccaccaccaccaGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACGCCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCtt                                                                                    >  2:1388137/1‑150 (MQ=255)
                                                                     accaccaccaGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGc                                                                                  >  2:605902/1‑150 (MQ=255)
                                                                             caGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATgtgt                                                                          <  1:605902/150‑1 (MQ=255)
                                                                             caGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATgtgt                                                                          <  2:941893/150‑1 (MQ=255)
                                                                                    tAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGc                                                                   <  1:941936/150‑1 (MQ=255)
                                                                                        ttttCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCa                                                               >  2:124545/1‑150 (MQ=255)
                                                                                        ttttCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCa                                                               >  2:1292162/1‑150 (MQ=255)
                                                                                             tcGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTaa                                                                       <  2:1158434/137‑1 (MQ=255)
                                                                                             tcGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTaa                                                                       >  1:1158434/1‑137 (MQ=255)
                                                                                              cgccgcTTTTTAGCGGCAAAATCTGGTACGGTACACCTGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTAc                                                         <  1:1724984/150‑1 (MQ=255)
                                                                                              cgccgcTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTAc                                                         <  1:1467701/150‑1 (MQ=255)
                                                                                                  gcTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGtt                                                     >  1:126672/1‑150 (MQ=255)
                                                                                                   cTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTg                                                    <  2:73423/150‑1 (MQ=255)
                                                                                                       ttAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCAt                                                <  1:1292162/150‑1 (MQ=255)
                                                                                                            ggCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGc                                           >  1:1100335/1‑150 (MQ=255)
                                                                                                            ggCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGc                                           <  1:1769447/150‑1 (MQ=255)
                                                                                                               aaaaTCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGcaac                                        <  1:1023810/150‑1 (MQ=255)
                                                                                                                aaaTCTGGTACGGTACGCCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGcaaca                                       <  1:1388137/150‑1 (MQ=255)
                                                                                                                aaaTCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGcaaca                                       <  1:884388/150‑1 (MQ=255)
                                                                                                                     tGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAacaacatc                                  <  1:16167/150‑3 (MQ=255)
                                                                                                                     tGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACAcc                                  <  2:641287/150‑1 (MQ=255)
                                                                                                                     tGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACAcc                                  <  2:180536/150‑1 (MQ=255)
                                                                                                                      ggtacggtacACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCt                                 >  2:618057/1‑150 (MQ=255)
                                                                                                                       gtacggtacACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCTCGATGGCATCGTCCCACGGCGTCATTCCTTGCTCAAGTCTTAACGGCAGCAATTTACCGTTGTCATCTGGCAACAACACCtc                                >  2:1975084/1‑149 (MQ=255)
                                                                                                                         acggtacACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGAc                              <  1:1731820/150‑1 (MQ=255)
                                                                                                                          cggtacACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACt                             <  2:1772067/150‑1 (MQ=255)
                                                                                                                           ggtacACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTTCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTg                            >  1:569235/1‑150 (MQ=255)
                                                                                                                           ggtacACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTg                            >  2:590737/1‑150 (MQ=255)
                                                                                                                            gtacACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGc                           >  1:1953872/1‑150 (MQ=255)
                                                                                                                               cacCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCggg                        <  1:1377673/150‑1 (MQ=255)
                                                                                                                                acCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGc                       >  1:1136405/1‑150 (MQ=255)
                                                                                                                                 ccGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCa                      >  1:20707/1‑150 (MQ=255)
                                                                                                                                      tAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGaaa                 >  1:1238267/1‑150 (MQ=255)
                                                                                                                                           tgtCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAACGTgg            <  2:733234/150‑1 (MQ=255)
                                                                                                                                            gtCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAACGTGGa           <  2:1437690/150‑1 (MQ=255)
                                                                                                                                              gcctGCCTGCGGGCGCGGGCAGGCCTTTATCAAAACTCCACTGCGTGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGtt                                                                         <  2:366761/83‑1 (MQ=255)
                                                                                                                                              gcctGCCTGCGGGCGCGGGCAGGCCTTTATCAAAACTCCACTGCGTGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGtt                                                                         >  1:366761/4‑86 (MQ=255)
                                                                                                                                                  gcccgcgcccgcAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAACGTGGAGGCAAt     <  1:348526/150‑1 (MQ=255)
                                                                                                                                                    ccgcgcccgcAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAACGTGGAGGCAATaa   >  1:250691/1‑150 (MQ=255)
                                                                                                                                                     cgcgcccgcAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAACGTGGAGGCAATaaa  >  1:1076879/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACAGCGTAAACGTCTCCAGCAGGCGCTGCTGTTCCGGGATCATCAACTGGCGCAGATTCCCCGGCTCCACCACCACCAGCCCGTAGGTTTTCTCGCCGCTTTTTAGCGGCAAAATCTGGTACGGTACACCGGGTAACGTGTCGGTGCCCGCGCCCGCAGGCAGGCCTTTATCAAAACTCCACTGCGCGATGGCATCGTCCCACGGCGTCATTCCTTGCGGATGTGTTAACGGCTGCAATTTACCGTTGTCATCGGGCAACAACACCTGACTGCGGGCATGAAACGTGGAGGCAATAAA  >  NC_000913/722505‑722803

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