Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,484,897 (C)6→7 coding (524/3594 nt) evgS → hybrid sensory histidine kinase in two‑component regulatory system with EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,484,8911.C96.9% 223.4 / ‑0.5 65coding (518/3594 nt)evgShybrid sensory histidine kinase in two‑component regulatory system with EvgA
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (34/29);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

CGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCAT  >  NC_000913/2484756‑2485017
                                                                                                                                        |                                                                                                                              
cgccgcCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccc                                                                                                                             >  2:8461/1‑139 (MQ=255)
  ccgccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccccc                                                                                                                           >  2:932586/1‑139 (MQ=255)
  ccgccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccccc                                                                                                                           >  2:423740/1‑139 (MQ=255)
     ccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccccccca                                                                                                                        <  1:8461/139‑3 (MQ=255)
       tctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGTCCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAgcccccccag                                                                                                                      >  1:986890/1‑135 (MQ=255)
       tctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                      <  1:932586/139‑5 (MQ=255)
        ctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGTGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccagac                                                                                                                     >  1:320589/1‑134 (MQ=255)
        ctTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAccccccccaga                                                                                                                     <  2:986890/139‑6 (MQ=255)
               atGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTaa                                                                                                               <  1:1535824/139‑1 (MQ=255)
                tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                      >  2:1540527/1‑126 (MQ=255)
                tGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                      <  1:1540527/130‑5 (MQ=255)
                     ttGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATtcatc                                                                                                         >  1:624886/1‑139 (MQ=255)
                           cAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCAt                                                                                                   >  1:880079/1‑139 (MQ=255)
                             accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCACTGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                      <  1:1094150/117‑5 (MQ=255)
                             accaaaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccaga                                                                                                                      >  2:1094150/1‑113 (MQ=255)
                                  accaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCaaa                                                                                            >  1:1287072/1‑139 (MQ=255)
                                  accaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccca                                                                                                                       >  1:1570789/1‑108 (MQ=255)
                                  accaTTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccccccca                                                                                                                       <  2:1570789/111‑4 (MQ=255)
                                     aTTGATTATTACCTTTCCGGCGCTTGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGc                                                                                         <  1:107243/139‑1 (MQ=255)
                                      ttGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGca                                                                                        <  1:423740/139‑1 (MQ=255)
                                               aCCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAtct                                                                              >  1:548043/1‑139 (MQ=255)
                                                    tCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCttt                                                                          >  1:252503/1‑139 (MQ=255)
                                                    tCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGATAATTCATCAATCATTTCCAAAAGCAACAATTATCTCttt                                                                          >  2:1116900/1‑139 (MQ=255)
                                                           cTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                        <  2:717168/134‑1 (MQ=255)
                                                           cTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                        >  1:717168/1‑134 (MQ=255)
                                                           cTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAc                                                                                                                               <  2:38603/78‑1 (MQ=255)
                                                           cTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAc                                                                                                                               >  1:38603/1‑78 (MQ=255)
                                                             ggTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                       <  1:574718/117‑1 (MQ=255)
                                                             ggTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                       >  2:574718/1‑117 (MQ=255)
                                                              gTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAACACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAgg                                                                                                                  >  2:1116891/1‑89 (MQ=255)
                                                              gTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAgg                                                                                                                  <  1:1116891/89‑1 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                        >  1:240989/1‑121 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                        <  2:240989/121‑1 (MQ=255)
                                                                        ttCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAt                                                      <  2:1287072/139‑1 (MQ=255)
                                                                         tCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                       <  1:485152/105‑1 (MQ=255)
                                                                         tCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaa                                                                                       >  2:485152/1‑105 (MQ=255)
                                                                           aCGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAg                                                   >  1:1469920/1‑139 (MQ=255)
                                                                                tCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCatcaatca                                                                                                    <  2:16765/85‑1 (MQ=255)
                                                                                tCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCatcaatca                                                                                                    >  1:16765/1‑85 (MQ=255)
                                                                                tCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                 >  1:1393004/1‑136 (MQ=255)
                                                                                tCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                 >  1:548179/1‑136 (MQ=255)
                                                                                tCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                 <  2:1393004/136‑1 (MQ=255)
                                                                                tCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                 <  2:548179/136‑1 (MQ=255)
                                                                                  aaTGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGt                                            >  1:445929/1‑139 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                        >  1:1460143/1‑109 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                        <  2:1460143/109‑1 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAAt                                                                   >  2:1228713/1‑114 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAAt                                                                   <  1:1228713/114‑1 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                          >  1:860135/1‑139 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                          >  2:809491/1‑139 (MQ=255)
                                                                                    tGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGtct                                          >  1:1179249/1‑139 (MQ=255)
                                                                                       gACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAg                                       <  2:252503/139‑1 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                             <  1:1324933/115‑1 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTatat                                                             >  2:1324933/1‑115 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAtt                                                     >  1:44761/1‑123 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAtt                                                     <  2:44761/123‑1 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCt                                     <  2:320589/139‑1 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAc                                                                                                                               <  2:632759/48‑1 (MQ=255)
                                                                                         ccGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAc                                                                                                                               >  1:632759/1‑48 (MQ=255)
                                                                                            cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                 <  2:1214673/44‑1 (MQ=255)
                                                                                            cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                 >  1:1214673/1‑44 (MQ=255)
                                                                                            cTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTggg                                  >  2:1454592/1‑139 (MQ=255)
                                                                                             ttACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGc                                 <  1:809491/139‑1 (MQ=255)
                                                                                                   cACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                 >  1:236355/1‑117 (MQ=255)
                                                                                                   cACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                 <  2:236355/117‑1 (MQ=255)
                                                                                                    aCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCatcaatca                                                                                                    <  2:1305158/65‑1 (MQ=255)
                                                                                                    aCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCatcaatca                                                                                                    >  1:1305158/1‑65 (MQ=255)
                                                                                                    aCCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGATGGGCACAATGa                          >  1:383296/1‑139 (MQ=255)
                                                                                                         aaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGca                                                                                        >  2:204453/1‑72 (MQ=255)
                                                                                                         aaCCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGca                                                                                        <  1:204453/72‑1 (MQ=255)
                                                                                                           ccAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                 >  1:739594/1‑109 (MQ=255)
                                                                                                           ccAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCa                                                 <  2:739594/109‑1 (MQ=255)
                                                                                                            cAGTAAATATTGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTa                  <  2:1179249/139‑1 (MQ=255)
                                                                                                             aGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAtt                <  2:548043/139‑1 (MQ=255)
                                                                                                                      tGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                       <  2:1583140/76‑1 (MQ=255)
                                                                                                                      tGCTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTAc                                                                       >  1:1583140/1‑76 (MQ=255)
                                                                                                                        cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                            >  2:1611616/1‑85 (MQ=255)
                                                                                                                        cTCGGGTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATc                                                            <  1:1611616/85‑1 (MQ=255)
                                                                                                                            ggTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAAcatcat  >  1:543635/1‑139 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGt                                            <  1:635933/96‑1 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGt                                            >  2:635933/1‑96 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTgg                                   >  2:1152217/1‑105 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTgg                                   <  1:1152217/105‑1 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAtt                        <  1:1407304/116‑1 (MQ=255)
                                                                                                                             gTAGCAAATTACCCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGAtt                        >  2:1407304/1‑116 (MQ=255)
                                                                                                                                        |                                                                                                                              
CGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCAT  >  NC_000913/2484756‑2485017

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