Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 994,634 A→G 6.0% W67R (TGG→CGG) ‡ ssuC ← aliphatic sulfonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913994,6340AG 6.0% 248.4 / 5.7 83W67R (TGG→CGG) ‡ssuCaliphatic sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (50/28);  new base G (2/3);  total (52/31)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

ACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTA  >  NC_000913/994486‑994783
                                                                                                                                                    |                                                                                                                                                     
aCCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAACAGCCAATCAGCGCCCGCCAg                                                                                                                                                      >  2:319883/1‑150 (MQ=255)
  cAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGa                                                                                                                                                    >  2:269266/1‑150 (MQ=255)
      gggATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTg                                                                                                                                                <  1:480368/150‑1 (MQ=255)
         aTCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTGAtt                                                                                                                                             <  1:1188566/150‑2 (MQ=255)
         aTCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATc                                                                                                                                             >  2:979564/1‑150 (MQ=255)
          tCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCg                                                                                                                                            >  2:648551/1‑150 (MQ=255)
           cAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTGAttgg                                                                                                                                           <  1:1042183/150‑4 (MQ=255)
           cAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGc                                                                                                                                           >  2:785249/1‑150 (MQ=255)
             gcgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCa                                                                                                                                         >  1:495400/1‑150 (MQ=255)
              cgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCACCCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCac                                                                                                                                        >  1:453869/1‑149 (MQ=255)
              cgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAg                                                                                                                                        >  2:303823/1‑150 (MQ=255)
               gcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCCACATCGACCCCAGCGATCCGCCAATTGAAAAGCCCATCAGCGACCGCCAGGTGCTGATCGCCAGa                                                                                                                                       >  2:191565/1‑150 (MQ=255)
                      gCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGcc                                                                                                                                <  2:1225683/150‑1 (MQ=255)
                         gCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCcaca                                                                                                                             >  2:534701/1‑150 (MQ=255)
                            cGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCCACCCCAGCGATCCGCCCATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGtt                                                                                                                          >  2:74797/1‑150 (MQ=255)
                               tGCGCAACATCTGAATTGATGTATACAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCACTCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgc                                                                                                                       <  1:138741/150‑1 (MQ=255)
                               tGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgc                                                                                                                       <  1:143640/150‑1 (MQ=255)
                                 cgcAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCACCCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgccc                                                                                                                     >  1:258969/1‑149 (MQ=255)
                                 cgcAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCCGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgccg                                                                                                                     >  2:32158/1‑150 (MQ=255)
                                     aCATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTgg                                                                                                                 >  1:496279/1‑150 (MQ=255)
                                     aCATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTgg                                                                                                                 <  1:1237506/150‑1 (MQ=255)
                                          tGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGgag                                                                                                            >  2:516998/1‑150 (MQ=255)
                                               tGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGt                                                                                                       >  2:735265/1‑150 (MQ=255)
                                                 aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                                     <  1:333679/150‑1 (MQ=255)
                                                 aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                                     >  1:1058682/1‑150 (MQ=255)
                                                 aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                                     <  2:399480/150‑1 (MQ=255)
                                                  ggTATCGAGCAGCAGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCCGCGATCCGCCAATTGAAAAGCCTATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCa                                                                                                    <  2:484813/150‑1 (MQ=255)
                                                  ggTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCa                                                                                                    <  1:1240465/150‑1 (MQ=255)
                                                       ccAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACg                                                                                               >  2:1109676/1‑150 (MQ=255)
                                                         agcagcCGTTCTCCCCAGCGAGATAGCCCGCTACTCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGCAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGcc                                                                                             >  2:623981/1‑150 (MQ=255)
                                                                 gtcaccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCACCGATCCGGCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGAAGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTcaccac                                                                                     <  1:1002276/146‑1 (MQ=255)
                                                                   ctcCCAAGCGAGATAGCCAGCTAATCAGTCCCAAAATCAACCCCACCGATCCGCCAAATGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTGGCCGCTGGCGGAGAGGGTCCAGCCCGCCGTCACCACTc                                                                                   >  2:150666/1‑150 (MQ=255)
                                                                     ccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCt                                                                                 >  1:647256/1‑150 (MQ=255)
                                                                     ccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCt                                                                                 <  1:419529/150‑1 (MQ=255)
                                                                        cAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTcc                                                                              >  2:680587/1‑150 (MQ=255)
                                                                          gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCCCAGTTCCCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCgg                                                                            >  2:1145642/1‑150 (MQ=255)
                                                                          gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCgg                                                                            <  1:1113118/150‑1 (MQ=255)
                                                                          gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCgg                                                                            <  2:1283830/150‑1 (MQ=255)
                                                                            gcgaTAGCCCGCTAATCAGTCCCAAAATCAACCCCACCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCCCAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCATCCGGTg                                                                          <  2:1119119/148‑1 (MQ=255)
                                                                            gagaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTg                                                                          <  1:1203682/150‑1 (MQ=255)
                                                                             agaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAACCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGa                                                                         <  1:1119041/150‑1 (MQ=255)
                                                                              gaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGaa                                                                        <  1:724780/150‑1 (MQ=255)
                                                                               aTAGCCCGCTTATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATGAGCGCCCGCCAGTAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACGACTCCTTCCGGTGAAg                                                                       <  2:508868/150‑1 (MQ=255)
                                                                               aTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAg                                                                       <  2:519425/150‑1 (MQ=255)
                                                                                  gCCCGCTAATCAGTCCCAAAATCAACCCCCGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCGGGGAAGGcc                                                                    >  1:620822/1‑149 (MQ=255)
                                                                                  gCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCa                                                                    >  2:480548/1‑150 (MQ=255)
                                                                                  gCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCa                                                                    >  1:1155344/1‑150 (MQ=255)
                                                                                  gCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCa                                                                    >  1:34213/1‑150 (MQ=255)
                                                                                   cccGCTAATCAGGCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCGGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCaa                                                                   <  1:1316589/150‑1 (MQ=255)
                                                                                       cTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATa                                                               <  1:269266/150‑1 (MQ=255)
                                                                                            cAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGt                                                          >  2:477799/1‑150 (MQ=255)
                                                                                                cccAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCTGGTGAAGGCAAAATACGCGTCGAc                                                      >  1:447409/1‑150 (MQ=255)
                                                                                                 ccAAAATCAACCCCAGCGCTCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGCGCGTCCAGAACGCCGTCACCACTCCTTCCGGGGAAGGCAAAATACGCGTCGAcc                                                     >  2:304202/1‑149 (MQ=255)
                                                                                                 ccAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcc                                                     >  2:277297/1‑149 (MQ=255)
                                                                                                  cAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcag                                                    <  2:647256/150‑1 (MQ=255)
                                                                                                   aaaaTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGCcagc                                                   >  2:961216/1‑150 (MQ=255)
                                                                                                     aaTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagcca                                                 <  2:184904/150‑1 (MQ=255)
                                                                                                      aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGGGCTGATCGCCAGATGCTGCAACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCGTGTGCAGGCAAAATACGCGTCGAcagccag                                                >  2:1167636/1‑150 (MQ=255)
                                                                                                      aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagccag                                                >  1:302837/1‑150 (MQ=255)
                                                                                                      aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagccac                                                >  2:479478/1‑149 (MQ=255)
                                                                                                            cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACCGCCAGCCAAcc                                          >  2:30258/1‑150 (MQ=255)
                                                                                                            cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAAcc                                          >  2:77077/1‑150 (MQ=255)
                                                                                                            cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTGCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCCCTCCTTCCGGTGAAGGCATACTACGCGACGCCAGCCAGCCAAcc                                          >  2:215651/1‑150 (MQ=255)
                                                                                                              cAGCGATCCGCCAAGTGAAACGCCAATCAGCGCCCGCCAGGAGCGGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGCGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCga                                        <  1:1175630/150‑1 (MQ=255)
                                                                                                                     ccGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCCACCGAGGAGGcc                                 >  2:730798/1‑150 (MQ=255)
                                                                                                                        ccAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGt                              >  1:1214279/1‑150 (MQ=255)
                                                                                                                         cAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGtt                             >  2:1126936/1‑150 (MQ=255)
                                                                                                                           aTTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCCGCCAACCGAGGAGGCCAGTTGc                           >  1:972322/1‑150 (MQ=255)
                                                                                                                            ttGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGcc                          >  1:497874/1‑150 (MQ=255)
                                                                                                                            ttGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGcc                          <  2:729963/150‑1 (MQ=255)
                                                                                                                             tGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACCGCCAGCCCACCGAGGAGGCCAGTTGccc                         >  1:1435/1‑149 (MQ=255)
                                                                                                                               aaaaGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCcaca                       >  1:387874/1‑150 (MQ=255)
                                                                                                                                aaaGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCacac                      <  1:648551/150‑1 (MQ=255)
                                                                                                                                 aaGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACAcc                     <  2:495400/150‑1 (MQ=255)
                                                                                                                                  aGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCg                    >  1:1047991/1‑150 (MQ=255)
                                                                                                                                   ggcgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCACCCGAGGAGGCCAGTTGCCACACCGc                   >  2:647886/5‑150 (MQ=255)
                                                                                                                                   gCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACCCCGc                   >  2:36688/1‑150 (MQ=255)
                                                                                                                                     cgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCACCGCTGGCGGAGCGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGccc                 >  2:1135920/3‑149 (MQ=255)
                                                                                                                                          aGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGCCAAAAGACGCGTCGACAGCCCGCCAACCGAGGAGGCCCGGTGCCACCCCGCCACGATg            >  2:784463/1‑150 (MQ=255)
                                                                                                                                                cgccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCAcc      <  2:557747/150‑1 (MQ=255)
                                                                                                                                                  ccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCgg    >  1:969670/1‑150 (MQ=255)
                                                                                                                                                  ccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCgg    <  1:516998/150‑1 (MQ=255)
                                                                                                                                                  ccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCgg    >  2:1273107/1‑150 (MQ=255)
                                                                                                                                                    aGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTa  <  2:118877/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTA  >  NC_000913/994486‑994783

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